【ATcoder s8pc_3 F】 壽司

http://s8pc-3.contest.atcoder.jp/tasks/s8pc_3_f (題目連接)

題意

　　有一個長度爲$N$的數列$A$，初始爲$0$。$Q$次操做，每次兩個參數$x,y$。

1. 在$A[1]...A[x]$中找出最小的數，若是有多個找編號最小的，假設爲$u$。
2. $A[u]++$。
3. 重複這個過程$y$次。

　　輸出最後的$A$序列。

Solution

　　能夠發現，數列$A$確定是單調不降的，那就很是好作了。

　　用一個線段樹維護數列，區間賦值，支持區間求和。在詢問區間中二分找到一個位置$pos$，$pos$到$x$的數改爲與$pos-1$大小相同所須要花費的次數不超過$y$。剩餘的次數所有加上就好。

細節

　　注意特判$x=1$和$A[x-1]-A[x]>=y$的狀況。

代碼

// ATcoder 
#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<cmath>
#define LL long long
#define inf (1ll<<30)
#define free(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout)
using namespace std;

const int maxn=100010;
int n,Q;
struct node {int l,r;LL s,tag;}tr[maxn<<2];

void pushdown(int k) {
	int l=k<<1,r=k<<1|1;LL w=tr[k].tag;tr[k].tag=0;
	tr[l].s=(tr[l].r-tr[l].l+1)*w;tr[l].tag=w;
	tr[r].s=(tr[r].r-tr[r].l+1)*w;tr[r].tag=w;
}
void modify(int k,int s,int t,LL val) {
	int l=tr[k].l,r=tr[k].r,mid=(l+r)>>1;
	if (l==s && r==t) {tr[k].s=val*(tr[k].r-tr[k].l+1);tr[k].tag=val;return;}
	if (tr[k].tag) pushdown(k);
	if (t<=mid) modify(k<<1,s,t,val);
	else if (s>mid) modify(k<<1|1,s,t,val);
	else modify(k<<1,s,mid,val),modify(k<<1|1,mid+1,t,val);
	tr[k].s=tr[k<<1].s+tr[k<<1|1].s;
}
void build(int k,int s,int t) {
	tr[k].l=s;tr[k].r=t;
	if (s==t) return;
	int mid=(s+t)>>1;
	build(k<<1,s,mid);
	build(k<<1|1,mid+1,t);
}
LL query(int k,int s,int t) {
	int l=tr[k].l,r=tr[k].r,mid=(l+r)>>1;
	if (l==s && r==t) return tr[k].s;
	if (tr[k].tag) pushdown(k);
	if (t<=mid) return query(k<<1,s,t);
	else if (s>mid) return query(k<<1|1,s,t);
	else return query(k<<1,s,mid)+query(k<<1|1,mid+1,t);
}

int main() {
	scanf("%d%d",&n,&Q);
	build(1,1,n);
	for (int x,i=1;i<=Q;i++) {
		LL y;
		scanf("%d%lld",&x,&y);
		if (x==1) {modify(1,1,1,query(1,1,1)+y);continue;}
		if (query(1,x-1,x-1)-query(1,x,x)>=y) {modify(1,x,x,query(1,x,x)+y);continue;}
		int l=2,r=x,pos;
		while (l<=r) {
			int mid=(l+r)>>1;
			if (query(1,mid-1,mid-1)*(x-mid+1)-query(1,mid,x)<=y) r=mid-1,pos=mid;
			else l=mid+1;
		}
		LL val=query(1,pos-1,pos-1),num=x-pos+1;
		y-=val*num-query(1,pos,x);
		modify(1,pos,x,val);++num;
		if (y/num) modify(1,pos-1,x,val+y/num);
		if (y%num) modify(1,pos-1,pos-2+y%num,val+y/num+1);
	}
	for (int i=1;i<=n;i++) printf("%lld\n",query(1,i,i));
	return 0;
}