【BZOJ4293 】【PA2015】Siano(思维+线段树区间覆盖与线段树上二分)

有个弟弟性质,一开始没注意到,就是所有草的高度排名是永远不变的

有了这个性质后,很容易想到将原序列通过二分找到那个临界点。然后全局加,答案是右半边的和减去原本的和。再对右边进行区间覆盖

注意,代码细节很多,比如线段树上二分,只需要维护一个max就能省很多事,还有pushdown时,f与cover标记的处理不是常规的先f在cover,原因留给你们思考一下。

#include<bits/stdc++.h>
#define re register
#define ll long long
const ll N=500005;
using namespace std;
template<class T>
inline void read(T &x)
{
	x=0; int f=1;
	static char ch=getchar();
	while(!isdigit(ch))	{if(ch=='-') f=-1;ch=getchar();} 
	while(isdigit(ch))	x=x*10+ch-'0',ch=getchar();
	x*=f;
}
inline void write(ll x)
{
	if(x<0) { putchar('-'); x=-x;}
	if(x>9)	write(x/10);
	putchar(x%10+'0');
}
ll n,m;
ll a[N],sum[N];
struct SegmentTree
{
	ll l,r;
	ll maxv,cover,f,sum;
}tree[N*4];
inline void pushup(ll now)
{
	tree[now].sum=tree[now<<1].sum+tree[now<<1|1].sum;
	tree[now].maxv=tree[now<<1|1].maxv;
}
inline void build(ll now,ll l,ll r)
{
	tree[now].l=l; tree[now].r=r; tree[now].cover=-1;
	if(l==r) return; 
	ll m=(l+r)>>1;
	build(now*2,l,m);
	build(now*2+1,m+1,r);
}
inline void modify(ll now,ll t)
{
	tree[now].f+=t;
	tree[now].maxv+=a[tree[now].r]*t;
	tree[now].sum+=(sum[tree[now].r]-sum[tree[now].l-1])*t;
}
inline void pushdown(ll now)
{
	ll f=tree[now].f; ll cover=tree[now].cover;
	if(tree[now].cover!=-1)
	{
		tree[now<<1].maxv=tree[now<<1|1].maxv=cover;
		
		tree[now<<1].sum=cover*(tree[now<<1].r-tree[now<<1].l+1); 
		tree[now<<1|1].sum=cover*(tree[now<<1|1].r-tree[now<<1|1].l+1);
		
		tree[now<<1].cover=tree[now<<1|1].cover=cover;
		tree[now<<1].f=tree[now<<1|1].f=0;
		tree[now].cover=-1;
	}
	if(f)
	{
		modify(now<<1,f); modify(now<<1|1,f);
		tree[now].f=0;
	}
}
inline ll find(ll now,ll x)
{
	if(tree[now].l==tree[now].r) return tree[now].l;
	pushdown(now);
	if(tree[now<<1].maxv>=x)	return find(now<<1,x);
	else return find(now<<1|1,x);
}
inline ll Query(ll now,ll l,ll r)
{
	if(l<=tree[now].l&&tree[now].r<=r)	return tree[now].sum;
	pushdown(now); ll m=(tree[now].l+tree[now].r)>>1; ll ans=0;
	if(l<=m) ans+=Query(now<<1,l,r);
	if(r>m) ans+=Query(now<<1|1,l,r);
	return ans;
} 
inline void Cover(ll now,ll l,ll r,ll x)
{
	if(l<=tree[now].l&&tree[now].r<=r)
	{
		tree[now].maxv=x;
		tree[now].cover=x;
		tree[now].sum=(tree[now].r-tree[now].l+1)*x;
		tree[now].f=0;
		return;
	}
	pushdown(now);
	ll m=(tree[now].l+tree[now].r)>>1;
	if(l<=m) Cover(now<<1,l,r,x);
	if(r>m) Cover(now<<1|1,l,r,x);
	pushup(now);
}
signed main()
{
	read(n); read(m);
	for(re int i=1;i<=n;i++) read(a[i]);
	sort(a+1,a+n+1);
	for(re int i=1;i<=n;i++) sum[i]=sum[i-1]+a[i];
	build(1,1,n);
	ll last=0;
	for(int i=1;i<=m;i++)
	{
		ll d,b;
		read(d); read(b);
		modify(1,d-last); last=d; 
		ll pos=find(1,b);
		if(tree[1].maxv<b)
		{
			puts("0");
			continue;
		}
		ll ans=Query(1,pos,n);
		write(ans-(b*(n-pos+1))); putchar('\n');
		Cover(1,pos,n,b);
	}
	return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值