HDU 5726 GCD

题目分析

这道题时间给的比较长，但是处理不好很容易超时，自己刚拿到这道题只觉得求一段区间的gcd好求，但是问有多少个跟该区间的gcd一样，这个我感觉不怎么好处理，然后看了网上的做法，打ST表，这个很简单，就是一个动态规划，不理解的可以学一下，很简单的，第一步处理完之后很多需要处理，刚开始我是这样想的，处理出给定的区间的gcd之后，然后利用gcd的性质（选定一个点作为左端点，再不断枚举右端点会发现区间的gcd值单调递减），这样就可以二分了，但是对于通过枚举左端点，然后二分右端点，很明显这样可以得到结果，但是有Q次询问，Q的数量达到了100000，因此毫无疑问此做法超时。因此做法是把所有区间的gcd值全部求一遍，当然不是暴力直接求一遍，而是通过二分求一遍，这样一来，时间复杂度降为真正的nlogn。当然这一步处理需要你对二分掌握比较好。

#include <map>
#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
const int maxn = 1e5+5;
typedef long long LL;
int dp[maxn][20],ans[maxn];
map <int, LL> MP;

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

void ST_build(int n){  // ST表
    for(int j = 1; (1<<j) <= n; j++)
        for(int i = 1; i+(1<<j)-1 <= n; i++)
            dp[i][j] = gcd(dp[i][j-1], dp[i+(1<<(j-1))][j-1]);
}

int query(int l,int r){ //查询
    int cnt = log2(r-l+1), len = 1<<cnt;
    return gcd(dp[l][cnt], dp[r-len+1][cnt]);
}

int main(){
    int T,N,Q;
    scanf("%d", &T);
    for(int kase = 1; kase <= T; kase++){
        scanf("%d", &N);
        for(int i = 1; i <= N; i++) scanf("%d", &dp[i][0]);
        MP.clear();
        ST_build(N); 
        int l,r;
        scanf("%d", &Q);
        for(int i = 1; i <= Q; i++){
            scanf("%d%d", &l, &r);
            MP.insert(make_pair(ans[i] = query(l,r), 0));  //初始化map
        }
        for(int i = 1; i <= N; i++){
            int beg = i;
            for(int j = beg; j <= N; ){
                int val = query(beg, j);
                l = j;
                r = N;
                while(l < r){  //二分法找右端点
                    int mid = (l+r+1)>>1;
                    if(query(beg, mid) < val)
                        r = mid-1;
                    else
                        l = mid;
                }
                if(MP.find(val) != MP.end()) MP[val] += (l-j+1);
                j = r+1;
            }
        }
        printf("Case #%d:\n", kase);
        for(int i = 1; i <= Q; i++)
            printf("%d %I64d\n", ans[i], MP[ans[i]]);
    }
    return 0;
}