hdu 5869 区间内不同的GCD数（离线+树状数组）

题目：http://acm.split.hdu.edu.cn/showproblem.php?pid=5869
题意：

给定n个数的序列，有q个询问，每次询问一个区间中所有子区间所形成不同的gcd的数量。

分析：

这个题的本质和HDU 3333是一样的，HDU 3333要求计算区间内不同的数字有几个。
这题稍微变了一下，相当于原来扫描到i的之后是更新a[i]的情况，现在是更新log级别个数的数字（因为以i为结尾的区间，最多只有log级别种不同的gcd）。

代码：

#include <cstdio>#include <iostream>#include <vector>#include <algorithm>#include <cstring>#include <string>#include <map>#include <cmath>#include <queue>#include <set>using namespace std;typedef long long ll;const int N = 1000000 + 9;int c[N], a[N], pre[N], last[N], ans[N], n, q;int gcd (int a, int b) {    return b == 0 ? a : gcd (b, a % b);}int lowbit (int x) {    return x & (-x);}void update (int x, int val) {    for (int i = x; i <= n; i += lowbit (i) )        c[i] += val;}int getsum (int x) {    int ans = 0;    for (int i = x; i >= 1; i -= lowbit (i) )        ans += c[i];    return ans;}int query (int l, int r) {    return getsum (r) - getsum (l - 1);}struct node {    int l, id;};vector<node>inv[N];int main() {    while (~scanf ("%d%d", &n, &q) ) {        int l, r;        for (int i = 1; i <= n; i++) scanf ("%d", &a[i]), inv[i].clear();        for (int i = 1; i <= q; i++) {            scanf ("%d%d", &l, &r);            inv[r].push_back ( (node) {l, i});        }        memset (last, 0, sizeof (last) );        memset (c, 0, sizeof (c) );        for (int i = 1; i <= n; i++) {            for (int j = i, x = a[j]; j > 0; j --, x = gcd (x, a[j]) ) {                if (j > last[x]) {                    if (last[x] != 0) update (last[x], -1);                    update (j, 1);                    last[x] = j;                }                if (x == 1) break;            }            for (int j = 0; j < inv[i].size(); j++)                ans[inv[i][j].id] = query (inv[i][j].l, i);        }        for (int i = 1; i <= q; i++)            printf ("%d\n", ans[i]);    }    return 0;}/*Sample Input5 31 3 4 6 93 52 51 5Sample Output666*/