HDU 5869 Different GCD Subarray Query (离线处理 树状数组)

题意：给你一个长度为n的数组，q次询问区间[L, R], 求区间内的所有子段的不同gcd值有多少种.

思路:有个知识：固定一个数为右端点，区间所能形成的gcd最多有logn个。所以可以预处理出每个数为右端点所能形

成的gcd，相同gcd取左端点靠右的。这样对询问的R排序就能离线处理，用树状数组维护即可。

代码：

#include<bits/stdc++.h>
using namespace std;
const int maxn = 1e6+15;
int n, q, tree[maxn];
int a[maxn], ans[maxn], vis[maxn];
typedef pair<int, int> pp;
vector<pp> g[maxn];
struct node
{
    int l, r, id;
    bool operator < (const node &a) const
    {
        if(r == a.r) return l < a.l;
        return r < a.r;
    }
}qq[maxn];

int lowbit(int x)
{
    return x&(-x);
}

void update(int pos, int val)
{
    while(pos < maxn)
    {
        tree[pos] += val;
        pos += lowbit(pos);
    }
}

int query(int pos)
{
    int res = 0;
    while(pos)
    {
        res += tree[pos];
        pos -= lowbit(pos);
    }
    return res;
}

void solve()
{
    for(int i = 1; i <= n; i++)
    {
        g[i].push_back(pp(a[i], i));
        int pre = a[i];
        for(int j = 0; j < g[i-1].size(); j++)
        {
            int tmp = __gcd(a[i], g[i-1][j].first);
            if(tmp != pre)
                g[i].push_back(pp(tmp, g[i-1][j].second)), pre = tmp;
        }
    }
    int cur = 0;
    for(int i = 1; i <= q; i++)
    {
        while(cur < qq[i].r)
        {
            cur++;
            for(int j = 0; j < g[cur].size(); j++)
            {
                int tmp = g[cur][j].first;
                int pos = g[cur][j].second;
                if(vis[tmp]) update(vis[tmp], -1);
                vis[tmp] = pos;
                update(vis[tmp], 1);
            }
        }
        ans[qq[i].id] = query(qq[i].r)-query(qq[i].l-1);
    }
    for(int i = 1; i <= q; i++)
        printf("%d\n", ans[i]);
}

int main(void)
{
    while(cin >> n >> q)
    {
        memset(vis, 0, sizeof(vis));
        memset(tree, 0, sizeof(tree));
        for(int i = 0; i <= n; i++)
            g[i].clear();
        for(int i = 1; i <= n; i++)
            scanf("%d", &a[i]);
        for(int i = 1; i <= q; i++)
            scanf("%d%d", &qq[i].l, &qq[i].r), qq[i].id = i;
        sort(qq+1, qq+1+q);
        solve();
    }
    return 0;
}