hdu 4630 No Pain No Game

树状数组+离线处理，输入结束后按照L从到小排序，然后在计算约数的时候，把约束出现次数大于1的存入树状数组，然后根据R进行查询，最后全部输出

#include <cstdio>
#include <cstring>
#include <cstdlib>
#define N 50005

struct node
{
    int l,r,pos;
};

node v[N];
int n,m;
int num[N],arr[N],p[N],q[N];

int max(int x,int y)
{
    return x>y?x:y;
}

int cmp(const void *p1,const void *p2)
{
    return ((node*)p2)->l - ((node*)p1)->l;
}

void add(int cur,int k)
{
    for(;cur<=n;cur+=(cur&(-cur))) arr[cur]=max(arr[cur],k);
}

int sch(int cur)
{
    int s=0;
    for(;cur>0;cur-=(cur&(-cur))) s=max(s,arr[cur]);
    return s;
}

int main()
{
  //  freopen("in.txt","r",stdin);
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&n);
        for(int i=1; i<=n; ++i) scanf("%d",&num[i]);
        scanf("%d",&m);
        for(int i=0; i<m; ++i) scanf("%d%d",&v[i].l,&v[i].r),v[i].pos=i;

        for(int i=0; i<N; ++i) arr[i]=p[i]=0;
        qsort(v,n,sizeof(node),cmp);

        for(int i=0,j=n; i<m; )
        {
            for(;j>0&&j>=v[i].l;--j)
                for(int k=1; k*k<=num[j]; ++k)
                    if(num[j]%k==0)
                    {
                        if(p[k]) add(p[k],k);
                        p[k]=j;
                        if(k!=num[j]/k&&p[num[j]/k]) add(p[num[j]/k],num[j]/k);
                        p[num[j]/k]=j;
                    }
            for(;i<m&&v[i].l>j;++i) q[v[i].pos]=sch(v[i].r);
        }

        for(int i=0;i<m;++i) printf("%d\n",q[i]);
    }
    return 0;
}