[区间GCD预处理 树状数组 离线] HDU 5869 Different GCD Subarray Query

首先确定一个右端点 向左做后缀gcd的值是不超过log的 因为gcd必然递减 每次至少除以2
那么可以对每个右端点预处理出来
然后要求区间不同的gcd个数 这里要用到1878: [SDOI2009]HH的项链的技巧
用树状数组离线处理
按右端点的顺序处理询问
令pos[x]等于x这个值出现的最靠右的位置 显然这里是对计算贡献最有利的 然后用树状数组维护下 对于某个位置 有多少个pos的取值是在这里的

#include<cstdio>#include<cstdlib>#include<algorithm>#include<cstring>#define cl(x) memset(x,0,sizeof(x))using namespace std;typedef long long ll;inline char nc(){  static char buf[100000],*p1=buf,*p2=buf;  if (p1==p2) { p2=(p1=buf)+fread(buf,1,100000,stdin); if (p1==p2) return EOF; }  return *p1++;}inline int read(int &x){  char c=nc(),b=1;  for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1; else if (c==EOF) return 0;  for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b; return 1;}const int N=100005;int n,m,a[N];struct abcd{  int l,r,d;  abcd(int l=0,int r=0,int d=0):l(l),r(r),d(d) {}}G[N][35];int cnt[N];inline void Init(){  for (int i=1;i<=n;i++){    for (int j=1;j<=cnt[i-1];j++)      G[i][j]=G[i-1][j],G[i][j].d=__gcd(G[i][j].d,a[i]);    cnt[i]=cnt[i-1]+1;    G[i][cnt[i]]=abcd(i,i,a[i]);    int pnt=0;    for (int j=1;j<=cnt[i];j++){      if (G[i][j].d!=G[i][pnt].d)    G[i][++pnt]=G[i][j];      else    G[i][pnt].r=G[i][j].r;    }    cnt[i]=pnt;  }}namespace BIT{  int maxn; ll c[N];  inline void init(int n){ maxn=n; for (int i=1;i<=n;i++) c[i]=0; }  inline void add(int x,int r){ for (int i=x;i<=maxn;i+=i&-i) c[i]+=r; }  inline ll sum(int x) { int ret=0; for (int i=x;i;i-=i&-i) ret+=c[i]; return ret; }  inline ll sum(int l,int r) { return sum(r)-sum(l-1); }}struct event{  int l,r,idx;  bool operator < (const event &B) const{ return r<B.r; }}eve[N];int tot;ll ans[N];int pos[1000005];int main(){  int T,l,r,Case=0;  freopen("t.in","r",stdin);  freopen("t.out","w",stdout);  while (read(n)){    read(tot);     for (int i=1;i<=n;i++) read(a[i]);    Init(); BIT::init(n);    for (int i=1;i<=tot;i++) read(eve[i].l),read(eve[i].r),eve[i].idx=i;    sort(eve+1,eve+tot+1);    int pnt=1;    for (int i=1;i<=n;i++){      for (int j=1;j<=cnt[i];j++){    int d=G[i][j].d,r=G[i][j].r;    if (pos[d]) BIT::add(pos[d],-1);    pos[d]=max(r,pos[d]); BIT::add(pos[d],1);      }      while (pnt<=tot && eve[pnt].r==i)    ans[eve[pnt].idx]=BIT::sum(eve[pnt].l,eve[pnt].r),pnt++;    }    for (int i=1;i<=tot;i++)      printf("%I64d\n",ans[i]);    cl(pos);  }  return 0;}