LonLife-ACM 1129 - 喵哈哈村的战斗魔法师丶坏坏い月

原题链接

1129 - 喵哈哈村的战斗魔法师丶坏坏い月

Time Limit：3s Memory Limit：256MByte

Submissions：381Solved：85

DESCRIPTION

坏坏い月是月大叔的ID，他是一个掌握者772002种魔法的物理系战士，最擅长的技能就是搞事。今天他又要开始搞事了。

给你$n$个数，你需要实现一下操作：

1. l r v ,在[l,r]区间内找到第一个大于等于v的数,输出这个数的下标,如果找不到的话，请输出-1噢

2. l r v，让[l,r]区间所有数增加v

INPUT

输入第一行包含一个正整数 $t(1\le t\le 100)$ ,表示有t组数据对于每组数据：第一行包含两个整数 $n(1\le n\le 100000)$, $q(1\le q\le 100000)$,表示数的个数，以及询问的个数。第二行包含 $n$个整数 $a_i(1\le a_i\le 1000000000)$接下来q行，每行四个整数 $opt(1\le opt\le 2),l,r(1\le l\le r\le n),v(1\le v\le 1000000000)$

OUTPUT

对于每个询问，输出一行表示答案.

SAMPLE INPUT

15 31 2 3 4 51 1 2 32 1 2 31 1 2 3

SAMPLE OUTPUT

-11

常见的数据结构中，我们选择使用分块来处理这道题。
我们将n个数，分成sqrt(n)块，每块里面的元素最多有sqrt(n)个元素。对于每一个块，我们维护Upd[i]表示这个块整个增加了多少，Max[i]表示这个块内的最大值是多少。
对于查询和更新，如果完全囊括了块的话，我们就只对Upd[i]和Max[i]进行更新，否则的话，我们就暴力更新这个块的元素即可。

#include <bits/stdc++.h>
#define maxn 100005
using namespace std;
typedef long long ll;

ll d[maxn], add[1005], maxs[1005];
int m;
int get(int k) {
	if(k <= m * m)
		return k % m ? k / m + 1 : k / m;
	return m;
}
void update(int l, int r, ll p, int h) {
	for(int i = l; i <= r; i++) {
		d[i] += p;
		maxs[h] = max(maxs[h], d[i]+add[h]);
	}
}
int find(int l, int r, ll v, int h) {
	for(int i = l; i <= r; i++) {
		if(d[i] + add[h] >= v)
		 return i;
	}
	return -1;
}
int main() {
//	freopen("in.txt", "r", stdin);
	int t;
	scanf("%d", &t);
	while(t--) {
		memset(add, 0, sizeof(add));
		memset(maxs, 0, sizeof(maxs));
		int q, n;
		scanf("%d%d", &n, &q); 
		for(int i = 1; i <= n; i++)
		 scanf("%lld", d+i);
		m = sqrt(n);
		for(int i = 1; i <= m; i++) {
			for(int j = (i-1)*m+1; j <= i*m; j++) {
				maxs[i] = max(maxs[i], d[j]);
			}
		}
		if(m * m < n) {
			for(int i = m * m + 1; i <= n; i++) {
				maxs[m] = max(maxs[m], d[i]);
			}
		}
		while(q--) {
			int t, l, r;
			ll v;
			scanf("%d%d%d%lld", &t, &l, &r, &v);
			int f1 = get(l);
			int f2 = get(r);
			if(t == 2) {
				if(f1 == f2) {
					update(l, r, v, f1);
				}
				else {
					update(l, f1*m, v, f1);
					update((f2-1)*m+1, r, v, f2);
					for(int i = f1+1; i < f2; i++) {
						add[i] += v;
						maxs[i] += v;
					}
				}
			}
			else {
				if(f1 == f2) {
					int k = find(l, r, v, f1);
					printf("%d\n", k);
				}
				else {
					int k = find(l, f1*m, v, f1);
					if(k != -1){
					 printf("%d\n", k);
					 continue;
					}
					int sign = 0;
					for(int i = f1+1; i < f2; i++) {
						if(maxs[i] >= v) {
							sign = 1;
							printf("%d\n", find((i-1)*m+1, i*m, v, i));
							break;
						}
					}
					if(sign == 0) {
						printf("%d\n", find((f2-1)*m+1, r, v, f2));
					}
				}
			}
		}
	}
	return 0;
}