[POI2011]MET-Meteors

单调性是显然的，暴力也是显然的：即对于每个国家都进行一次二分答案。复杂度：nm logn

考虑对于一整段连续的国家整体二分。

在solve的过程中，对于1-mid次陨星雨增加的值，用线段树或树状数组来维护。进行区间累加以后，判断每一个国家是否能在前mid次陨星雨以后收集到相应的陨星。对于能与不能，分别放在左右区间，进行下一步递归。

这里注意两个小优化：

1.区间修改单点查询的线段树，可以利用差分思想后用树状数组维护。

2.对于每次solve，不用一遍又一遍地把1-mid次陨星雨都重新累加一遍。对于之前mid次无法收集到期望陨星的国家，在之前就不断累减去mid次陨星雨收集的陨星，这样的话在单次solve中，只需要把l-mid做一次即可。

#include <bits/stdc++.h>
#define int long long
#define lowbit(x) x&-x
using namespace std;
const int N=3e5+5;
int m,n,x,k;
vector<int>vec[N];
int sum[N],ans[N];
struct node{int k,id;}q[N],q1[N],q2[N];
struct number{int x,y,w;}num[N];

inline void change(int x,int v)
{
	while (x<=n+1)
	{
		sum[x]+=v;
		x+=lowbit(x);
	}
}
inline void modify(int l,int r,int w)
{
	if (l<=r)
	{
		change(l,w);
		change(r+1,-w);
	}
	else
	{
		change(l,w);
		change(n+1,-w);
		change(1,w);
		change(r+1,-w);
	}
}
inline int query(int x)
{
	int res=0;
	while (x)
	{
		res+=sum[x];
		x-=lowbit(x);
	}
	return res;
}

void solve(int ql,int qr,int l,int r)
{
	if (ql>qr) return;
	if (l==r)
	{
		for (register int i=ql; i<=qr; ++i) ans[q[i].id]=l;
		return;	
	}
	int mid=l+r>>1;
	for (register int i=l; i<=mid; ++i) modify(num[i].x,num[i].y,num[i].w);
	int p1=0,p2=0;
	for (register int i=ql; i<=qr; ++i)
	{
		int tot=0;
		bool jay=false;
		for (register int j=0; j<vec[q[i].id].size(); ++j)
		{
			tot+=query(vec[q[i].id][j]);
			if (tot>=q[i].k) {jay=true; break;}
		}
		if (jay) q1[++p1]=q[i];	
		else q[i].k-=tot,q2[++p2]=q[i];
	}
	for (register int i=l; i<=mid; ++i) modify(num[i].x,num[i].y,-num[i].w);
	for (register int i=1; i<=p1; ++i) q[ql+i-1]=q1[i];
	for (register int i=1; i<=p2; ++i) q[ql+p1+i-1]=q2[i];
	solve(ql,ql+p1-1,l,mid);
	solve(ql+p1,qr,mid+1,r);
}

inline int read()
{
	int ret=0,ff=1; char ch=getchar();
	while (!isdigit(ch)) {if (ch=='-') ff=-ff; ch=getchar();}
	while (isdigit(ch)) {ret=(ret<<3)+(ret<<1)+ch-'0'; ch=getchar();}
	return ret*ff;
}
inline void write(int x)
{
	if (x<0) {x=-x; putchar('-');}
	if (x>9) write(x/10);
	putchar(x%10+48); 
}
inline void writeln(int x)
{
	write(x);
	putchar('\n');
}

signed main(){
	m=read(); n=read();
	for (register int i=1; i<=n; ++i)
	{
		x=read();
		vec[x].push_back(i);
	}
	for (register int i=1; i<=m; ++i) q[i].k=read(),q[i].id=i;
	k=read();
	for (register int i=1; i<=k; ++i) num[i].x=read(),num[i].y=read(),num[i].w=read();
	k++; num[k].x=1; num[k].y=n; num[k].w=1e9;
	solve(1,m,1,k);
	for (register int i=1; i<=m; ++i) 
	{
		if (ans[i]==k) 
		{
			putchar('N'); putchar('I'); putchar('E'); putchar('\n');
		}
		else writeln(ans[i]);
	}
return 0;	
}