BZOJ5177: [Jsoi2013]贪心的导游

Description

南京有一条著名的购物街。购物街嘛，就是一排整齐的商店啦~

导游小Z每次都会把游客团带到购物街里走一段，然后选择一个商店进去购物。小Z接待的游客都是购物狂，他们恨

不得将店内的商品洗劫一空，也就是说，只要他们能买，就一定会继续买（钱够不够你不用考虑，他们都有信用卡

可以透支）。但是有一点，他们都非常讲究平等、很谦虚，每个人都不能忍受比别人多买什么东西或者少买什么东

西，于是他们每个人最后买的商品数量都是一样的。这虽然导致他们没办法每次都把商店搬空，但是每次已经给店

家带来一大笔生意了，店家已经非常感谢了！为了表示感谢，店家决定把游客们买完之后剩下来那几件没卖掉的商

品就送给导游小Z了。贪心的小Z自然希望自己能获赠的商品数量越大越好啦~现在告诉你这一排共n个商店（标号为

0到n-1）每个商店里的商品总数，每次小Z会带一批共p个游客的旅游团，到其中u号商店和v号商店之间逛一逛，请

你帮小Z在所逛的商店区间内选择一个，告诉小Z他最多能获赠多少件商品。

Input

第一行，包含两个整数n、m，分别表示商店个数、小Z带来的旅游团个数。

接下来一行，包含n个整数ai（i=0,1,……,n-1），表示第i个商店的商品总数。

接下来m行，每行三个整数u、v、p（0≤u,v≤n-1，2≤p≤1000）

表示这个旅游团逛u号商店和v号商店之间的商店（包含u、v），且这个旅游团的人数为p。

n≤1000000， m≤50000， 0≤ai≤1000， 2≤p≤1000

Output

共输出m行，每行一个整数，第i行输出第i个旅游团购物后，小Z最多能获赠的商品数量。

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
第一个旅游团， 2 个人， 0 号商店到 1 号商店的区间。若去 1 号商店，共 2 件商品，每人买 1 件，剩 0 件
。若去 2 号商店，共 4 件商品，每人买 2 件，剩 0 件。所以，小 Z 最多获赠 0 件。
第二个旅游团， 3 个人，小 Z 选择带他们去 4 号商店，共 8 件商品，每人买 2 件商品（因为每人 3 件不够）
，剩下 2 件，小 Z 最多获赠就是 2 件。 可以验证去其它商店小 Z 最多获赠的商品不会达到 2 件。

Solution

考虑到a_i和p均<=1000.所以可以暴力枚举p的倍数，对每个倍数求那段区间内的前驱。
那么问题就变成了询问区间[l,r]的前驱。这个东西可以用主席树+二分维护。（二分查找前驱，主席树维护有没有这个数）。
复杂度是(O(500mlogn))。
不过第一次打主席树+二分，基本是照着别人的代码和自己的板子打的...还是得多练几道题。

#include <bits/stdc++.h>
#define ll long long
#define inf 0x3f3f3f3f
#define il inline

namespace io {

#define in(a) a = read()
#define out(a) write(a)
#define outn(a) out(a), putchar('
')

#define I_int ll
inline I_int read() {
    I_int x = 0, f = 1;
    char c = getchar();
    while (c < '0' || c > '9') {
        if (c == '-') f = -1;
        c = getchar();
    }
    while (c >= '0' && c <= '9') {
        x = x * 10 + c - '0';
        c = getchar();
    }
    return x * f;
}
char F[200];
inline void write(I_int x) {
    if (x == 0) return (void) (putchar('0'));
    I_int tmp = x > 0 ? x : -x;
    if (x < 0) putchar('-');
    int cnt = 0;
    while (tmp > 0) {
        F[cnt++] = tmp % 10 + '0';
        tmp /= 10;
    }
    while (cnt > 0) putchar(F[--cnt]);
}
#undef I_int

}
using namespace io;

using namespace std;

#define N 12000000
const int M = 1000;

int n, m, a[1000010], b[1000010];
int lc[N], rc[N], root[1000010], cnt, sum[N];

int upd(int l, int r, int pos, int rt) {
	int now = ++cnt;
	lc[now] = lc[rt]; rc[now] = rc[rt]; sum[now] = sum[rt];
	if(l == r) {sum[now]++; return now;}
	int mid = (l + r) >> 1;
	if(pos <= mid) lc[now] = upd(l, mid, pos, lc[rt]);
	else rc[now] = upd(mid + 1, r, pos, rc[rt]);
	sum[now] = sum[lc[now]] + sum[rc[now]];
	return now;
}

int ans = 0; bool flag = 0;

void find(int L, int R, int l, int r, int pos) {
	while(l < r) {
		int mid = (l + r) >> 1, x = sum[rc[R]] - sum[rc[L]];
		if(x > 0) l = mid + 1, L = rc[L], R = rc[R];
		else r = mid, L = lc[L], R = lc[R];
	}
	ans = max(ans, l - pos);
}

void solve(int L, int R, int l, int r, int vl, int vr, int pos) {
	if(sum[R] - sum[L] == 0) return;
	if(l == vl && r == vr) { flag = 1; find(L, R, l, r, pos); return; }
	int mid = (l + r) >> 1;
	/*
	if(vr >= mid) solve(rc[L], rc[R], mid + 1, r, vl, vr, pos);
	if(flag) return;
	if(vl < mid) solve(lc[L], lc[R], l, mid, vl, vr, pos);
	*/
	if(vl > mid) solve(rc[L], rc[R], mid + 1, r, vl, vr, pos);
	else if(vr <= mid) solve(lc[L], lc[R], l, mid, vl, vr, pos);
	else {
		solve(rc[L], rc[R], mid + 1, r, mid + 1, vr, pos);
		if(flag) return;
		solve(lc[L], lc[R], l, mid, vl, mid, pos);
	}
}

int main() {
	in(n), in(m);
	for(int i = 1; i <= n; ++i) in(a[i]), root[i] = upd(0, M, a[i], root[i-1]);
	while(m--) {
		int l = read()+1, r = read()+1, p = read();
		if(l > r) swap(l, r);
		int k = 0; ans = 0;
		while(k <= M) {
			flag = 0;
			solve(root[l-1], root[r], 0, M, k, min(M, k+p-1), k);
			k += p;
			if(ans == p - 1) break;
		}
		outn(ans);
	}
//	outn(cnt);
//	for(int i = 1; i <= cnt; ++i)printf("%d ", sum[i]);
}