[Jsoi2013]贪心的导游 主席树+暴力

Description
给n个数，m个询问，每次询问l~r区间内的数mod p的最大值。

---

Sample Input
5 5
2 4 6 8 10
0 1 2
1 4 3
2 4 2
1 1 9
0 4 7

---

Sample Output
0
2
0
4
6

---

这道题考虑主席树。
我先看了p的范围和ai范围，于是只能想到一个大暴力。。。
就是枚举p的倍数，在一段一段里求最大值。。。
想一下其实时间复杂度最大O(500*log n*m)
差不多5亿吧，一秒跑1亿他又开了40s，ok吧。
然后你再想一下就比如说那些2、3之类很小的常数最大值非常有限，你判断一下如果当前答案已经最优秀了，那就break。
跑得还挺快的，10s就跑出来了。。。
权值线段树找最大值应该会吧，就是你找一个最右的有权值的点，先优先往右找，然后往左找。

---

#include <cstdio>
#include <cstring>

using namespace std;
int _max(int x, int y) {return x > y ? x : y;}

struct node {
    int l, r, lc, rc, c;
} t[21000000]; int cnt;
int rt[1100000];

void Link(int &u, int l, int r, int p) {
    if(!u) u = ++cnt;
    t[u].c++;
    if(l == r) return;
    int mid = (l + r) / 2;
    if(p <= mid) Link(t[u].lc, l, mid, p);
    else Link(t[u].rc, mid + 1, r, p);
}

void Merge(int &u1, int u2) {
    if(!u1 || !u2) {u1 = u1 + u2; return;}
    t[u1].c += t[u2].c;
    Merge(t[u1].lc, t[u2].lc);
    Merge(t[u1].rc, t[u2].rc);
}

int calc(int u1, int u2, int l, int r, int ll, int rr) {
    if(l == r) return l;
    int ls = t[t[u1].lc].c - t[t[u2].lc].c, rs = t[t[u1].rc].c - t[t[u2].rc].c;
    int mid = (l + r) / 2;
    if(rr <= mid) {
        if(ls == 0) return 0;
        return calc(t[u1].lc, t[u2].lc, l, mid, ll, rr);
    }
    else if(ll > mid) {
        if(rs == 0) return 0;
        return calc(t[u1].rc, t[u2].rc, mid + 1, r, ll, rr);
    }
    else {
        int k1 = 0, k2 = 0;
        if(rs != 0) {
            k2 = calc(t[u1].rc, t[u2].rc, mid + 1, r, mid + 1, rr);
            if(k2) return k2;
        }
        if(ls != 0) {
            k1 = calc(t[u1].lc, t[u2].lc, l, mid, ll, mid);
            if(k1) return k1;
        }
    }
}

int main() {
    int n, m; scanf("%d%d", &n, &m);
    for(int i = 1; i <= n; i++) {
        int x; scanf("%d", &x);
        Link(rt[i], 0, 1000, x);
        Merge(rt[i], rt[i - 1]);
    }
    for(int i = 1; i <= m; i++) {
        int l, r, p; scanf("%d%d%d", &l, &r, &p); l++; r++;
        int maxx = 0;
        for(int j = 1; j <= 1000 / p; j++) {
            if(maxx == p - 1) break;
            maxx = _max(maxx, calc(rt[r], rt[l - 1], 0, 1000, (j - 1) * p, j * p - 1) - (j - 1) * p);
        }
        if(1000 % p != 0 && maxx != p - 1) maxx = _max(maxx, calc(rt[r], rt[l - 1], 0, 1000, 1000 / p * p, 1000) - 1000 / p * p);
        printf("%d\n", maxx);
    }
    return 0;
}