HDU-4630-No Pain No Game

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int maxn=50100;
struct node
{
int l;
int r;
int pos;
bool operator <(const node& a)const
{
return l<a.l;
}
}q[maxn];
int n,m,a[maxn],pre[maxn],c[maxn],ans[maxn];
int lowbit(int x)
{
return x&(-x);
}
void Update(int i,int val)
{
for(;i<=n;c[i]=max(c[i],val),i+=lowbit(i));
}
int Sum(int i)
{
int ans=0;
for(;i>0;ans=max(ans,c[i]),i-=lowbit(i));
return ans;
}
int main()
{
int T;
scanf("%d",&T);
while(T--)
{
memset(pre,0,sizeof(pre));
memset(c,0,sizeof(c));
scanf("%d",&n);
for(int i=1;i<=n;i++)
scanf("%d",&a[i]);
scanf("%d",&m);
for(int i=1;i<=m;i++)
{
scanf("%d%d",&q[i].l,&q[i].r);
q[i].pos=i;
}
sort(q+1,q+m+1);
for(int i=m,j=n;i>0;)
{
while(j>=1&&j>=q[i].l)
{
for(int k=1;k*k<=a[j];k++)
if(a[j]%k==0)
{
if(pre[k])
Update(pre[k],k);
pre[k]=j;
if(k*k==a[j])
break;
if(pre[a[j]/k])
Update(pre[a[j]/k],a[j]/k);
pre[a[j]/k]=j;
}
j--;
}
while(q[i].l>j)
{
ans[q[i].pos]=Sum(q[i].r);
i--;
}
}
for(int i=1;i<=m;i++)
printf("%d\n",ans[i]);
}
return 0;
}