hdu 5726 GCD 倍增+ 二分

题目链接

给n个数， 定义一个运算f[l,r] = gcd(al, al+1,....ar)。 然后给你m个询问， 每次询问给出l, r。 求出f[l, r]的值以及有多少对l', r' 使得f[l, r] = f[l', r']。

第一个很简单， 用倍增的思想就可以了。 

然后是第二个， 我们枚举每一个左端点i， 显然f[i, j]是只降不增的。 那么我们可以二分找到所有使得f[i, j]下降的值j。 因为gcd每次至少变为原来的二分之一， 而ai最大为1e9. 所以最多只有log2(1e9)个这样的点。 我们都找出来然后放到map里就可以了。 具体看代码

#include <bits/stdc++.h>

using namespace std;
#define ll long longint n;
const int maxn = 1e5+5;
int a[maxn], f[maxn][20], mm[maxn];
map <int, ll> mp;
int gcd(int a, int b)
{
    return b?gcd(b, a%b):a;
}
void initrmq()
{
    mm[0] = -1;
    for(int i = 1; i <= n; i++) {
        mm[i] = ((i&(i-1))==0)?mm[i-1]+1:mm[i-1];
    }
    for(int j = 1; j < 17; j++) {
        for(int i = 1; i + (1<<j)-1 <= n; i++) {
            f[i][j] = gcd(f[i][j-1], f[i+(1<<(j-1))][j-1]);
        }
    }
}
int query(int l, int r)
{
    int k = mm[r-l+1];
    return gcd(f[l][k], f[r-(1<<k)+1][k]);
}
void pre()
{
    for(int i = 1; i <= n; i++) {
        int g = f[i][0];
        int L = i, tmp;
        while(L <= n) {
            int l = L, r = n;
            while(l <= r) {
                int mid = l+r>>1;
                if(query(i, mid) == g) {
                    tmp = mid;
                    l = mid+1;
                } else {
                    r = mid-1;
                }
            }
            mp[g] += (tmp-L+1);
            L = tmp+1;
            g = gcd(g, f[L][0]);
        }
    }
}
int main()
{
    int t, m, l, r;
    cin>>t;
    for(int casee = 1; casee <= t; casee++) {
        cin>>n;
        for(int i = 1; i <= n; i++) {
            scanf("%d", &a[i]);
            f[i][0] = a[i];
        }
        mp.clear();
        initrmq();
        pre();
        cin>>m;
        printf("Case #%d:
", casee);
        for(int i = 0; i < m; i++) {
            scanf("%d%d", &l, &r);
            int ans = query(l, r);
            printf("%d %lld
", ans, mp[ans]);
        }
    }
}