2016多校联合训练赛 第一场1004 GCD hdu 5726

题目连接： http://acm.hdu.edu.cn/showproblem.php?pid=5726
题意：给出一组序列，和一些区间[l, r]询问，要求求出[l, r]区间所有数的最大公约数（GCD），并给出这个序列中，有多少个区间的GCD也是这个值。
思路：比赛时一直在琢磨莫队算法和线段树的做法，但是想了很久觉得莫队可能会TLE，所以换了线段树的做法，但是可能是线段树写的太挫了，就继续TLE，无奈改成了RMQ的ST做法。
具体是这样的，我们用ST的做法在$O(nlog(n))$可以预处理出dp[i][j]，表示以第i个数为起点的，长度为$2^j$的区间的GCD，然后我们在以后的询问里就都可以$O(1)$的查询出给定区间的GCD了。然后由于数据范围是$0<=a_i <= 10^9$那么我们可以得出区间的GCD最多有$log(10^9)$种，因为它最多有$log(10^9)$种质因数。然后我们枚举起点，二分的找出其GCD对于答案的贡献，每个GCD的值需存在一个map里。

具体代码（写的不是一般的搓）：

#include <cstdio>
#include <map>
#include <cmath>
#include <iostream>
#include <algorithm>

typedef long long ll;
const int MAXN = 100000 + 10;
int dp[MAXN][70];
int a[MAXN];
int n;
std::map< int, ll > m;

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

inline void init(){
    for(int i = 1; i <= n; ++i){
        dp[i][0] = a[i];
    }
    for(int i = 1, x = 1; x <= n; ++i, x <<= 1){
        for(int j = 0; j + x <= n; ++j){
            dp[j][i] = gcd(dp[j][i - 1], dp[j + x][i - 1]);
        }
    }
    return;
}

inline int query(int l, int r){
    int k = log(r - l + 1) / log(2);
    return gcd(dp[l][k], dp[r - (1 << k) + 1][k]);
}

int main(){
    int cas, t = 0;
    scanf("%d", &cas);
    while(cas--){
        scanf("%d", &n);
        for(int i = 1; i <= n; ++i){
            scanf("%d", a + i);
        }
        init();
        m.clear();
        for(int i = 1; i <= n; ++i){
            for(int cur = i, d = a[i]; cur <= n; ){
                int l = cur, r = n;
                while(l < r){
                    int mid = (l + r + 1) >> 1;
                    if(query(l, r) == d){
                        l = mid;
                    }else{
                        r = mid - 1;
                    }
                }
                m[d] += l - cur + 1;
                cur = l + 1;
                d = gcd(d, a[l + 1]);
            }    
        }
        printf("Case #%d:\n", ++t);
        int q;
        scanf("%d", &q);
        while(q--){
            int l, r;
            scanf("%d %d", &l, &r);    
            int res = query(l, r);
            printf("%d %lld\n", res, m[res]);    
        }
    }
    return 0;
}