HDU 4630-No Pain No Game（线段树+离线处理）

题意：

给你n个数的序列a，q个询问，每个询问给l,r，求在下标i在[l，r]的区间任意两个数的最大公约数中的最大值

分析：

有了hdu3333经验，我们从左向右扫序列，如果当前数的约数在前面出现过，那这个约数可能就是最大的答案。所以我们枚举当前数的所有约数，用线段树维护区间最大值，查询序列离线处理保证查询的正确。

#include <map>
#include <set>
#include <list>
#include <cmath>
#include <queue>
#include <stack>
#include <cstdio>
#include <vector>
#include <string>
#include <cctype>
#include <complex>
#include <cassert>
#include <utility>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
using namespace std;
typedef pair<int,int> PII;
typedef long long ll;
#define lson l,m,rt<<1
#define pi acos(-1.0)
#define rson m+1,r,rt<<1|1
#define All 1,N,1
#define N 50010
#define read freopen("in.txt", "r", stdin)
const ll  INFll = 0x3f3f3f3f3f3f3f3fLL;
const int INF= 0x7ffffff;
const int mod =  1000000007;
//bef当前数的约数前一次出现的位置
int maxv[N*4],n,bef
,e,a
,maxn
;
struct que{
int l,r,id;
}q
;
bool cmp(que x,que y){
return x.r<y.r;
}
void pushup(int rt){
maxv[rt]=max(maxv[rt<<1],maxv[rt<<1|1]);
}
void build(int l,int r,int rt){
maxv[rt]=1;
if(l==r)return;
int m=(l+r)>>1;
build(lson);
build(rson);
}
void update(int p,int v,int l,int r,int rt){
if(l==r){
maxv[rt]=max(maxv[rt],v);
return;
}
int m=(l+r)>>1;
if(p<=m)update(p,v,lson);
else update(p,v,rson);
pushup(rt);
}
int query(int L,int R,int l,int r,int rt){
if(L<=l&&R>=r){
return maxv[rt];
}
int m=(l+r)>>1;
int num=0;
if(L<=m)num=max(num,query(L,R,lson));
if(R>m)num=max(num,query(L,R,rson));
return num;
}
void solve(){
sort(q+1,q+e+1,cmp);
memset(bef,0,sizeof(bef));
build(1,n,1);
int b=1;
for(int i=1;i<=n;++i){
for(int j=1;j*j<=a[i];j++){
if(a[i]%j==0){
if(bef[j]){
update(bef[j],j,1,n,1);
}
int tmp=a[i]/j;
if(bef[tmp]&&tmp!=j)
update(bef[tmp],tmp,1,n,1);
bef[j]=i;
bef[tmp]=i;
}
}
while(b<=e&&q[b].r==i){
if(q[b].l==q[b].r)
maxn[q[b].id]=0;
else{
maxn[q[b].id]=query(q[b].l,q[b].r,1,n,1);
}
b++;
}
}
for(int i=1;i<=e;++i)
printf("%d\n",maxn[i]);
}
int main()
{
int t;
scanf("%d",&t);
while(t--){
scanf("%d",&n);
for(int i=1;i<=n;++i)
scanf("%d",&a[i]);
scanf("%d",&e);
for(int i=1;i<=e;++i){
scanf("%d%d",&q[i].l,&q[i].r);
q[i].id=i;
}
solve();
}
return 0;
}