http://acm.hdu.edu.cn/showproblem.php?pid=5875
因为取膜一个比自己大的数不会改变,所以每次都找第一个比当前值大的数的位置,用线段树来实现,据说使用RMQ加二分会更优雅,常数小
线段树查询函数,线段树记录区间最小值,如果比查询元素小,则继续向下进行,每次优先递归左半部,找到位置就可以全部return掉了
#include<bits/stdc++.h>
using namespace std;
const int maxn=1e5+7;
int seg[maxn*4];
int a[maxn];
void pushup(int node)
{
seg[node]=min(seg[node<<1],seg[node<<1|1]);
return ;
}
void build(int l,int r,int node)
{
if(l==r)
{
seg[node]=a[l];
return ;
}
int m=(l+r)>>1;
build(l,m,node<<1);
build(m+1,r,node<<1|1);
pushup(node);
}
int flog=0,index,v;
void query(int x,int le,int ri,int l,int r,int node)
{
if(flog)
return ;
if(l==r)
{
index=l;
flog=1;
return ;
}
// cout<<l<<" "<<r<<endl;
int m=(l+r)>>1;
if(m>=le&&seg[node<<1]<=x)
query(x,le,ri,l,m,node<<1);
if(m<ri&&seg[node<<1|1]<=x)
query(x,le,ri,m+1,r,node<<1|1);
}
int main()
{
int T,n,i;
cin>>T;
while(T--)
{
cin>>n;
for(i=1;i<=n;i++)
{
scanf("%d",&a[i]);
}
build(1,n,1);
int m,l,r;
cin>>m;
while(m--)
{
scanf("%d %d",&l,&r);
int now=a[l];
while(l!=r)
{
// cout<<" "<<now<<" "<<l<<endl;
flog=0;
query(now,l+1,r,1,n,1);
if(flog)
l=index;
else
break;
now%=a[l];
}
printf("%d\n",now);
}
}
return 0;
}