RMQ预处理gcd

#include
using namespace std;
#define LL long long
const int MAXN=100009;
int st[MAXN][20];
mapans;
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<n)break;    //退出条件：从j位起长为(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<>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