HDU4630 线段树离线操作 && GCD性质

题意：

n个数，m个询问，每次查询一个区间的中任意两个数的最大GCD

思路：

离线操作，按照查询的右端点将查询从小到大排序。预处理出所有数的因子用二维vector保存。遍历整个这n个数，每个数遍历它的因子，如果这个因子在前面出现过，则将这个因子维护到线段树中，并将因子最新出现位更新，如果没有出现过则将当前此因子最新出现的位置标记为当前。线段树用来维护区间的最大因子。因为区间两次以上出现过的最大因子即为最大GCD，所以查询区间最大值即可。

#include<bits/stdc++.h>
using namespace std;
const int maxn=50005;
struct nn{
    int l,r,id;
}po[maxn];
struct node{
    int l,r,maxx;
}tree[maxn<<2];
vector<int>ve[maxn];
int a[maxn],pre[maxn],ans[maxn];
bool cmp(nn a,nn b){
    return a.r<b.r;
}
void init(){
    for(int i=1;i<maxn;i++){
        for(int j=i;j<maxn;j+=i){
            ve[j].push_back(i);
        }
    }
}
void build(int l,int r,int i){
//    cout<<l<<" "<<r<<endl;
    tree[i].l=l;tree[i].r=r;
    tree[i].maxx=0;
    if(l==r)return;
    int mid=l+r>>1;
    build(l,mid,i<<1);
    build(mid+1,r,i<<1|1);
}
void update(int pos,int val,int i){
    tree[i].maxx=max(val,tree[i].maxx);
    if(tree[i].l==tree[i].r)return;
    int mid=tree[i].l+tree[i].r>>1;
    if(mid>=pos)update(pos,val,i<<1);
    else update(pos,val,i<<1|1);
}
int query(int l,int r,int i){
    if(tree[i].l==l&&tree[i].r==r)
        return tree[i].maxx;
    int mid=tree[i].l+tree[i].r>>1;
    if(mid>=r) query(l,r,i<<1);
    else if(mid+1<=l) query(l,r,i<<1|1);
    else return max(query(l,mid,i<<1),query(mid+1,r,i<<1|1));
}
int main(){
//    ios::sync_with_stdio(0);
//    cin.tie(0);cout.tie(0);
    init();
    int t;scanf("%d",&t);
    while(t--){
        int n;scanf("%d",&n);
        build(1,n,1);
        for(int i=1;i<=n;i++){
            scanf("%d",&a[i]);
        }
        int m;scanf("%d",&m);
        for(int i=0;i<m;i++){
                scanf("%d%d",&po[i].l,&po[i].r);
//            cin>>po[i].l>>po[i].r;
            po[i].id=i;
        }
        sort(po,po+m,cmp);
//        for(int i=0;i<m;i++){
//            cout<<po[i].l<<" "<<po[i].r<<" "<<po[i].id<<endl;
//        }
        memset(pre,-1,sizeof(pre));
        for(int i=1,j=0;i<=n&&j<m;i++){
            for(int k=0;k<ve[a[i]].size();k++){
                int tmp=ve[a[i]][k];
//                cout<<tmp<<endl;
                if(pre[tmp]!=-1){
                    update(pre[tmp],tmp,1);
                }
                pre[tmp]=i;//一次都没出现过，第一次出现不用往线段树里加
            }
            while(j<m&&po[j].r==i){
                ans[po[j].id]=query(po[j].l,po[j].r,1);
                j++;
            }
//            cout<<i<<endl;
        }
        for(int i=0;i<m;i++){
                printf("%d\n",ans[i]);
//            cout<<ans[i]<<endl;
        }
    }
    return 0;
}