[CF 475D] CGCDSSQ (RMQ)

题目链接：http://codeforces.com/contest/475/problem/D

是昨天晚上的CF题目，题意是给定你n个数，问你所有子区间内的最小公约数是x的个数是多少

 

问的康神，了解了ST表。

其实我也没太明白，看了看博文觉得还是不错的：http://hplonline20100103.blog.163.com/blog/static/1361364342010040044244/

然后就是logn去扫了

 

代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;

int ST[250100][19];

int getGCD(int l,int r){
    int u = 31-__builtin_clz(r-l+1);
    return __gcd(ST[l][u],ST[r+1-(1<<u)][u]);
}

map<int,LL> mp;

int main(){
    int n;
    scanf("%d",&n);
    for(int i=0;i<n;i++) scanf("%d",&ST[i][0]);
    for(int j=0;(1<<j)<=n;j++) for(int i=0;i<n;i++)
        ST[i][j+1]=__gcd(ST[i][j],ST[i+(1<<j)][j]);
    for(int i=0;i<n;i++){
        int ed = i;
        while( ed<n ){
            int g = getGCD(i,ed);
            int l=ed,r=n-1,ge=-1;
            while( l<=r ){
                if( r-l<=1 ){
                    if(getGCD(i,r)==g){ ge=r;break; }
                    if( getGCD(i,l)==g )ge=l;
                    break;
                }
                int m = l+r>>1;
                if( getGCD(i,m)==g ) l = m;
                else r = m;
            }
            mp[g]+=ge-ed+1;
            ed = ge+1;
        }
    }
    int Q;
    scanf("%d",&Q);
    while(Q--){
        int x;
        scanf("%d",&x);
        printf("%I64d\n",mp[x]);
    }
    return 0;
}