HDU5869 Different GCD Subarray Query(线段树&&离线处理) 2016大连网络赛

题意：给定一个a数组，每次询问一个区间[l,r]求这个区间内所有子区间的gcd的种类数目。

题解 ： 对于这个题目我们发现是求一段连续区间不同gcd的个数，一定注意是连续，我们注意到对于连续的区间 gcd 下降的非常快。每次这个数至少除以2，这样的话，这个级别就变成了 logn 级别的处理了。我们就可以预处理出这个数向左有多少个不同的gcd 。这个怎么处理呢 ？ 注意到 gcd (a,b,c) ＝gcd( gcd(a,b),c) 所以我们可以借助上一次的gcd 的值进行处理 如果两个gcd 的值相同我们保留最靠右的那个。这是为什么呢 ？这就和我们怎么处理查询有关系了。

怎么处理这个查询呢 ？ 我们先把所有的询问全部存下来并且按照右端点排序，排序完之后我们从1 到 n 扫描 每次维护 一个新 gcd 最靠右的位置让这个位置 ＋ 1 ，让这个gcd 上次出现的位置 －1 就可以了，每次查询一个 l 到 r 的区间和然后最后一起输出就可以了。

对于这种一个区间内多少个不同的数这种问题而言我们经常考虑离线处理这些数。

#include <iostream>
#include <algorithm>
#include <cstring>
#include <vector>
using namespace std;
const int maxn = 100010;
int a[maxn] = {0};
int n,q;
struct node {
    int pos;
    int gcd;
};
struct qq {
    int l,r,id;
}Q[maxn];
struct node1 {
    int l,r,sum;
}tr[maxn << 2];
int pre[10 * maxn] = {0};
vector <node> v[maxn];
int res[maxn] = {0};
bool cmp (const qq a,const qq b) {
    return a.r < b.r;
}
void build (int l,int r,int root) {
    tr[root].l = l;
    tr[root].r = r;
    tr[root].sum = 0;
    if (l == r) return ;
    int mid = (l + r) >> 1;
    build(l, mid, root << 1);
    build(mid + 1, r, root << 1 | 1);
}

void update (int pos,int add,int root) {
    if (tr[root].l == tr[root].r) {
        tr[root].sum += add;
        return ;
    }
    int mid = (tr[root].l + tr[root].r) >> 1;
    if (mid >= pos) update(pos, add, root << 1);
    else update(pos, add, root << 1 | 1);
    tr[root].sum = tr[root << 1].sum + tr[root << 1 | 1].sum;
}

int query (int l,int r,int root) {
    if (l <= tr[root].l && r >= tr[root].r)
        return tr[root].sum;
    int mid = (tr[root].l + tr[root].r) >> 1;
    int temp1 = 0,temp2 = 0;
    if (mid >= l) temp1 = query(l, r, root << 1);
    if (mid < r) temp2 = query(l, r, root << 1 | 1);
    return temp1 + temp2;
}

int main () {
    ios_base :: sync_with_stdio(false);
    while (cin >> n >> q) {
        v[0].clear();
        for (int i = 1;i <= n; ++ i) {
            cin >> a[i];
            v[i].clear();
        }
        memset (res,0,sizeof (res));
        build(0, n + 1, 1);
        for (int i = 1;i <= n ; ++ i) {
            int x = a[i];
            int y = i;
            int k = v[i - 1].size();
            for (int j = 0;j < k; ++ j) {
                int u = v[i - 1][j].gcd;
                if (__gcd(u,x) != x) {
                    v[i].push_back({y,x});
                    y = v[i - 1][j].pos;
                    x = __gcd (u,x);
                }
            }
            v[i].push_back({y,x});
        }
        for (int i = 1;i <= q; ++ i) {
            int l,r;
            cin >> l >> r;
            Q[i].id = i;
            Q[i].l = l;
            Q[i].r = r;
        }
        sort (Q + 1,Q + q + 1,cmp);
        int now = 1;
        memset (pre,0,sizeof (pre));
        for (int i = 1;i <= n; ++ i) {
            int k = v[i].size();
            for (int j = 0;j < k; ++ j) {
                int pos = v[i][j].pos;
                int gcd = v[i][j].gcd;
                update(pre[gcd], -1, 1);
                pre[gcd] = pos;
                update(pos, 1, 1);
            }
            while (Q[now].r == i) {
                int id = Q[now].id;
                res[id] = query(Q[now].l, i, 1);
                now ++;
            }
            if (now > q) break;
        }
        for (int i = 1;i <= q; ++ i)
            cout << res[i] << endl;
    }
    return 0;
}