CodeForces 272F - Ant colony线段树

题意：给出长度为n的串a，然后 给出q个询问 。对于每个询问 l,r，有一个计算其中对应的数的能量值的方法，比如power[l] 的计算方法就是 在区间[l,r]中能被 a[l]整除的数的个数。问这个区间中除去能量值等于r-l的个数还剩多少个。

因为个数 等于r-l  ，就相当于最小的数是 这个区间的最大公约数。就是求如果存在最小的数等于这个区间的最大公约数，最小的数的个数。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <climits>
#include <string>
#include <iostream>
#include <map>
#include <cstdlib>
#include <list>
#include <set>
#include <queue>
#include <stack>
#include <math.h>
using namespace std;
typedef long long LL;
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
const int maxn = 111111;
const int INF = 1000000000 + 2;
int sum[maxn << 2], Min[maxn << 2], ret[maxn << 2];
int a[maxn];
int gcd(int a, int b)
{
    if (a%b == 0) return b;
    return gcd(b, a%b);
}
void up(int rt)
{
    Min[rt] = min(Min[rt << 1], Min[rt << 1 | 1]);
    sum[rt] = gcd(sum[rt << 1], sum[rt << 1 | 1]);
    if (Min[rt << 1] == Min[rt << 1 | 1]) ret[rt] = ret[rt << 1] + ret[rt << 1 | 1];
    if (Min[rt << 1] > Min[rt << 1 | 1]) ret[rt] = ret[rt << 1 | 1];
    if (Min[rt << 1] < Min[rt << 1 | 1]) ret[rt] = ret[rt << 1];
}

void build(int l, int r, int rt)
{
    if (l == r){
        Min[rt] = sum[rt] = a[l]; ret[rt] = 1;
        return;
    }
    int mid = (l + r) >> 1;
    build(lson);
    build(rson);
    up(rt);
}

int askmin(int L, int R, int l, int r, int rt)
{
    int ans = INF;
    if (L <= l&&r <= R) return Min[rt];
    int mid = (l + r) >> 1;
    if (L <= mid) ans = askmin(L, R, lson);
    if (R > mid) ans = min(ans, askmin(L, R, rson));
    return ans;
}

int askgcd(int L, int R, int l, int r, int rt)
{
    if (L <= l&&r <= R) return sum[rt];
    int ans1 = -1; int ans2 = -1;
    int mid = (l + r) >> 1;
    if (L <= mid) ans1 = askgcd(L, R, lson);
    if (R>mid) ans2 = askgcd(L, R, rson);
    if (ans1 == -1) return ans2;
    if (ans2 == -1) return ans1;
    return gcd(ans1, ans2);
}

int askret(int L, int R, int m, int l, int r, int rt)
{
    int ans = 0;
    if (L <= l&&r <= R){
        if (m == Min[rt]) return ret[rt];
        return 0;
    }
    int mid = (l + r) >> 1;
    if (L <= mid) ans += askret(L, R, m, lson);
    if (R>mid) ans += askret(L, R, m, rson);
    return ans;
}

int main()
{
    int n, l, r, q;
    cin >> n;
    for (int i = 1; i <= n; i++)
        scanf("%d", &a[i]);
    build(1, n, 1);
    cin >> q;
    for (int i = 0; i < q; i++){
        scanf("%d%d", &l, &r);
        int m = askmin(l, r, 1, n, 1);
        int Gcd = askgcd(l, r, 1, n, 1);
        if (Gcd == m)printf("%d\n", r - l + 1 - askret(l, r, m, 1, n, 1));
        else printf("%d\n", r - l + 1);
    }
    return 0;
}