题意:
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;
}