题意:给你n个数,每次查询[l,r]区间里所有子区间不同gcd的个数。
我们先来看如何查询区间不同数的个数: spoj D-query 区间不同数个数 主席树||离线+树状数组
上面这道题讨论的是单个点,既对答案有贡献的是单个点,而这里是子区间要怎么处理。
同理,对每个gcd有多个子区间,定义子区间左端点为子区间大小,那么只有最大的那个子区间才是有效的。
这里还要用到gcd性质,离线后对当前点它只能和前i-1个后缀取gcd,以及它自己作为gcd,这里只会产生log(a[i])个gcd。
所以保存下来暴力枚举即可。
ps:就统计区间不同数来说,还有在线做法,既对每个前缀建一棵线段树(主席树)hdu 5790 prefix 主席树
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <map>
#include <algorithm>
using namespace std;
const int maxn=1e5+7;
struct node{
int l,r,id;
bool operator<(node t)const {
return r<t.r;
}
}q[maxn];
map<int,int> mp[3];
map<int,int> ::iterator it;
int a[maxn];
int ans[maxn];
int data[maxn];
int sum(int i){
int s=0;
while(i>0){
s+=data[i];
i-=i&-i;
}
return s;
}
void add(int i,int x){
while(i<maxn){
data[i]+=x;
i+=i&-i;
}
}
int main()
{
int n,m;
while(~scanf("%d%d",&n,&m)){
fill(data,data+n+5,0);
for(int i=0;i<3;i++)
mp[i].clear();
for(int i=1;i<=n;i++){
scanf("%d",&a[i]);
}
for(int i=1;i<=m;i++){
scanf("%d%d",&q[i].l,&q[i].r);
q[i].id=i;
}
sort(q+1,q+m+1);for(int i=1;i<=m;i++)
q[0].r=0;
for(int i=1;i<=m;i++){
for(int j=q[i-1].r+1;j<=q[i].r;j++){
for(it=mp[(j&1)^1].begin();it!=mp[(j&1)^1].end();it++){
int tmp=__gcd(a[j],it->first);
mp[j&1][tmp]=1;
if(!mp[2][tmp])mp[2][tmp]=mp[2][it->first],add(mp[2][it->first],1);
else if(mp[2][it->first]>mp[2][tmp]){
add(mp[2][tmp],-1);
add(mp[2][it->first],1);
mp[2][tmp]=mp[2][it->first];
}
}
mp[j&1][a[j]]=1;
if(mp[2][a[j]]){
add(mp[2][a[j]],-1);
add(j,1);
}
else add(j,1);
mp[2][a[j]]=j;
mp[(j&1)^1].clear();
}
ans[q[i].id]=sum(q[i].r)-sum(q[i].l-1);
}
for(int i=1;i<=m;i++){
printf("%d\n",ans[i]);
}
}
return 0;
}