hdu5726-GCD-ST表+二分

先用st表处理出所有l-r的GCD值，然后二分求得这些值一共出现了多少次。

#include<bits/stdc++.h>

#define inf 0x3f3f3f3f

const int maxn=100000;

using namespace std;

typedef pair<int,int> P;

int l,r;

int icase;

int t,n,q;

int a[maxn+1];

int dp[maxn+1][20];

map<int,long long> m;

int gcd(int a,int b){
   return b==0?a:gcd(b,a%b);
}

void RMQ(){
   for(int i=1;i<=n;i++){
        dp[i][0]=a[i];
   }
   for(int j=1;j<20;j++){
        for(int i=1;i<=n;i++){
                if(i+((1<<(j-1))-1)<=n)
                dp[i][j]=gcd(dp[i][j-1],dp[i+(1<<(j-1))][j-1]);
                //printf("%d
",dp[i][j]);
        }
   }
}

int query(int l,int r){
    int k=(int)log2(double(r-l+1));
    return gcd(dp[l][k],dp[r-(1<<k)+1][k]);
}

void Cal(){
   for(int i=1;i<=n;i++){
        int GCD=a[i],tail=i;
        while(tail<=n){
             int lb=tail,ub=n,mid=(lb+ub)>>1;
             while(ub-lb>1){
                mid=(lb+ub)>>1;
                if(query(i,mid)==GCD) lb=mid;
                else ub=mid;
             }
             if(query(i,ub)==GCD){
                m[GCD]+=(ub-tail+1);
                tail=ub+1;
                GCD=query(i,tail);
             } else if(query(i,mid)==GCD){
                m[GCD]+=(mid-tail+1);
                tail=mid+1;
                GCD=query(i,tail);
             }
               else{
                m[GCD]+=(lb-tail+1);
                tail=lb+1;
                GCD=query(i,tail);
             }
        }
   }
}

void solve(){
  scanf("%d",&n);
  for(int i=1;i<=n;i++){
        scanf("%d",&a[i]);
  }
  RMQ();
  Cal();
  scanf("%d",&q);
  printf("Case #%d:
",++icase);
  while(q--){
        scanf("%d%d",&l,&r);
        int res=query(l,r);
        printf("%d %I64d
",res,m[res]);
  }
}

void init(){
  m.clear();
  memset(dp,0,sizeof(dp));
}

int main()
{
    scanf("%d",&t);
    while(t--){
        init();
        solve();
    }
    return 0;
}