RMQ预处理gcd

#include<bits/stdc++.h>
using namespace std;
#define LL long long
const int MAXN=100009;
int st[MAXN][20];
map<int,LL>ans;
void RMQ(int n){                    //RMQ一般就是指ST表
    int Len=log2(n);                //得到了要维护的深度
    for(int i=1;i<=Len;i++){          //i深度，在输入的时候st[j][0]已经存好，这里i直接从1开始
        for(int j=1;j<=n;j++){        //j位置，注意输入时存的下标是从1~N
            if(j-1+(1<<i)>n)break;    //退出条件：从j位起长为(1<<i)区间最后一位置j+(1<<i)-1，故>n跳出
            st[j][i]=__gcd(st[j][i-1],st[j+(1<<(i-1))][i-1]);//st[j][i]指第j位起长为2^i长度区间GCD
        }
    }
}
int gcd(int l,int r){            //传入查找区间
    int Len=log2(r-l+1);        //在第几深度查找
    return __gcd(st[l][Len],st[r-(1<<Len)+1][Len]);//左右区间GCD的GCD即为所求（区间可部分覆盖）
}
//gcd从左端点开始要么是一长段不变 , 要么是直接减少一半以上
void init(int n){                 //就是用来更新ANS
    for(int i=1;i<=n;i++){        //枚举左界，二分求出每次变化的点
        int now=i;                //在以i为左界下，now是当前判到的变化点，即由i开始往右扫区间GCD有变化的点，左界i初始化为第一个
        while(now<=n){            //变化点未去到尽头就要继续，因为要把以I为左界的右边的区间全部扫完
            int l=now,r=n,x=now;    //解释x:在以i为左界下，以当前变化点now为右界的GCD记为G1
                                    //则x是在以i为左界下，能使GCD也是G1的最大右界
            while(l<=r){                    //二分未结束
                int mid=l+r>>1;             //二分区间中点
                if(gcd(i,mid)==gcd(i,now))  //如果i,mid区间与i,now区间GCD下相同
                    x=mid,l=mid+1;          //表明下一变化点在右半区间，则最大右界移到中点，然后左界右移
                else r=mid-1;               //否则下一变化点在左半区间，则右界左移，最大右界不动
            }
            ans[gcd(i,now)]+=x-now+1;       //一整段gcd都相同，最大右界减当前变化点的长度加进对应答案里计数
            now=x+1;        //下一次变化点就是当前变化点下GCD不变的最大右界
        }
    }
}
int main(){
    int t,ca=0;cin>>t;
    while(t--){
        printf("Case #%d:\n",++ca);ans.clear();
        int n;cin>>n;for(int i=1;i<=n;i++)cin>>st[i][0];//输入N个数，注意是从1开始存的
        RMQ(n);              //处理ST表
        init(n);             //处理GCD映射MAP（GCD值，GCD数量）
        int q;cin>>q;        //询问数
        while(q--){
            int l,r;cin>>l>>r;int g=gcd(l,r);    //找出区间GCD
            printf("%d %lld\n",g,ans[g]);        //输出区间找到的GCD及一样GCD的区间数
        }
    }
} 
RMQ预处理gcd（hdu 5726 GCD）
题意：n的数组 , q次询问 , 每次询问(l,r) , 求区间内所有数的gcd , 以及有多少个区间的gcd和这个区间相同(包括自己)
Sample Input
1
5
1 2 4 6 7
4
1 5
2 4
3 4
4 4
Sample Output
Case #1:
1 8
2 4
2 4
6 1