hdu5726-GCD-ST表+二分

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

 #include<bits/stdc++.h>

 #define inf 0x3f3f3f3f

 const int maxn=;

 using namespace std;

 typedef pair<int,int> P;

 int l,r;

 int icase;

 int t,n,q;

 int a[maxn+];

 int dp[maxn+][];

 map<int,long long> m;

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

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

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

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

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

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

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