HDU 5726 GCD（RMQ+二分）

题目链接：点击打开链接

题目大意：n个数，q次查询，每次查询输入l和r，输出区间GCD等于区间[l，r]的GCD的区间数。

解题思路：n的范围为1e5，可以先用ST预处理出所有区间的GCD，时间复杂度为O（nlogn）处理过之后，对于任意区间都可以实现O（1）的查询了。

然后就是计数的问题了，我们先定义一个map，map[key]表示GCD为key的区间个数。对于任意左端点固定的区间，随着右端点的增大，区间的GCD是单调非增的。

所以枚举左端点i，我们先设定一个j=i，然后左边界为j+1，右边界为n，开始二分查找GCD（i，l）==GCD（i，j）的最小下标l。然后向map[GCD[i，j]]添加区间[j，l]的长度，也就是j-l+1。然后将j更新为l+1，再重新二分求长度进行计数。

最后就可以进行对于每个GCD，就可以进行O（1）的查询了。

时间复杂度O（nlognlogn）

#include <set>
#include <map>
#include <cmath>
#include <stack>
#include <queue>
#include <vector>
#include <string>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long LL;
#define FIN freopen("in.txt", "r", stdin);
#define FOUT freopen("out.txt", "w", stdout);
#define lson l, mid, cur << 1
#define rson mid + 1, r, cur << 1 | 1
const int INF = 0x3f3f3f3f;
const int MAXN = 1e5 + 50;
const int MOD = 1e9 + 7;

int n, q;
LL num[MAXN], F[MAXN][35];
map<LL,LL> m;

void RMQ_init(int n)
{
    for (int i = 1; i <= n; i++)
        F[i][0] = num[i];
    for (int j = 1; (1 << j) <= n; 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]);
}

LL RMQ_query(int l, int r)
{
    int k = (int)log2(r - l + 1);
    return __gcd(F[l][k], F[r - (1 << k) + 1][k]);
}

int main()
{
#ifndef ONLINE_JUDGE
    FIN;
#endif // ONLINE_JUDGE
    int tcase;
    scanf("%d", &tcase);
    for (int c = 1; c <= tcase; c++)
    {
        scanf("%d", &n);
        for (int i = 1; i <= n; i++)
            scanf("%lld", &num[i]);
        RMQ_init(n);
        m.clear();
        for (int i = 1; i <= n; i++)
        {
            int j = i;
            LL tmpgcd = num[i];
            while (j <= n)
            {
                int l = j, r = n;
                while (l <= r)
                {
                    int mid = (l + r) / 2;
                    if (RMQ_query(i, mid) >= tmpgcd)
                        l = mid + 1;
                    else
                        r = mid - 1;
                }
                //二分结束时，l = h + 1
                //h是满足GCD[i，j] == GCD[i，h]的最下下标，所以区间长度就是h - j + 1，也就是l - j。
                //然后把j更新为l
                m[tmpgcd] += l - j; 
                j = l;
                if (j > n)
                    break;
                tmpgcd = RMQ_query(i, j);
            }
        }
        scanf("%d", &q);
        printf("Case #%d:\n", c);
        while (q--)
        {
            int a, b;
            scanf("%d%d", &a, &b);
            LL gcd = RMQ_query(a, b);
            printf("%lld %lld\n", gcd, m[gcd]);
        }
    }
    return 0;
}