【CF878E】Numbers on the blackboard(贪心)(并查集)

传送门


怎么说呢。。。这个贪心,难以形容。。。

题解:

首先考虑全为正数的时候,显然就是贪心从尾部开始向前合并,全为负数的时候,显然就是贪心每次取头部两个合并。

其实大胆猜想+证明可以得到这道题正确的贪心方式。

注意到贡献的表示形如 ∑ a i 2 k i \sum a_i2^{k_i} ai2ki,我们希望尽可能最大化正数的贡献。

对于一个序列,从最后一个位置开始考虑。

如果这个位置为正,与前面的合并,否则我们跳到前面的位置继续考虑。

显然这样最后合并出来除了第一个块可能是正数,剩下块全部都是负数,可以直接按照全部都是负数的情况合并。

我们发现这种情况好像可以拓展,考虑在已经得到答案的序列尾部加一个数可以直接按照上面的规则类似一个单调栈向前合并。

询问离线,右端暴力加点,并查集维护一下块的权值,以及块的哈希值的前缀和。

然后考虑左端点,在块内的决策显然是从左到右依次+1,那么类似字符串哈希维护一下哈希值就行了。


代码:

#include<bits/stdc++.h>
#include<tr1/unordered_map>
#define ll long long
#define re register
#define cs const

namespace IO{
	inline char gc(){
		static cs int Rlen=1<<22|1;
		static char buf[Rlen],*p1,*p2;
		return (p1==p2)&&(p2=(p1=buf)+fread(buf,1,Rlen,stdin),p1==p2)?EOF:*p1++;
	}
	template<typename T>T get(){
		char c;T num;bool f=0;
		while(!isdigit(c=gc()))f=c=='-';num=c^48;
		while(isdigit(c=gc()))num=(num+(num<<2)<<1)+(c^48);
		return f?-num:num;
	}inline int gi(){return get<int>();}
}
using namespace IO;

using std::cerr;
using std::cout;

cs int mod=1e9+7;
inline int add(int a,int b){return a+b>=mod?a+b-mod:a+b;}
inline int dec(int a,int b){return a-b<0?a-b+mod:a-b;}
inline int mul(int a,int b){ll r=(ll)a*b;return r>=mod?r%mod:r;}
inline int po(int a,int b){
	int r=1;for(;b;b>>=1,a=mul(a,a))
	if(b&1)r=mul(r,a);return r;
}
inline void Inc(int &a,int b){a+=b-mod;a+=a>>31&mod;}
inline void Dec(int &a,int b){a-=b;a+=a>>31&mod;}
inline void Mul(int &a,int b){a=mul(a,b);}

cs int N=5e5+7;
int n,Q,a[N],h[N],pw[N],ans[N];
struct Query{int l,r,i;}q[N];
int pre[N],len[N],bel[N],s[N];ll sum[N];
inline int get_h(int l,int r){return dec(h[l],mul(h[r+1],pw[r-l+1]));}
int gf(int x){while(x!=bel[x])x=bel[x]=bel[bel[x]];return x;}
void merge(int x,int y){
	bel[x]=y,pre[y]=pre[x];int len=x-pre[x];
	if(len>30&&sum[y]>0||sum[x]+sum[y]*(1ll<<len)>mod)sum[y]=mod;
	else sum[y]=sum[x]+sum[y]*(1ll<<len);
}

void work(){
	n=gi(),Q=gi();
	for(int re i=1;i<=n;++i)a[i]=gi();
	for(int re i=1;i<=Q;++i){int l=gi(),r=gi();q[i]=(Query){l,r,i};}
	std::sort(q+1,q+Q+1,[](cs Query &a,cs Query &b){return a.r<b.r;});
	for(int re i=n;i;--i)h[i]=add(add(h[i+1],h[i+1]),dec(a[i],0));
	for(int re i=1;i<=n;++i)bel[i]=i,pre[i]=i-1,sum[i]=a[i];
	pw[0]=1;for(int re i=1;i<=n;++i)pw[i]=add(pw[i-1],pw[i-1]);
	for(int re i=1,j=0;i<=n&&j<Q;++i){
		while(pre[i]&&sum[i]>=0)merge(pre[i],i);
		s[i]=add(s[pre[i]],add(get_h(pre[i]+1,i),get_h(pre[i]+1,i)));
		while(q[j+1].r==i){
			int b=gf(q[++j].l);
			ans[q[j].i]=add(dec(s[i],s[b]),get_h(q[j].l,b));
		}
	}
	for(int re i=1;i<=Q;++i)cout<<ans[i]<<"\n";
}

#undef zxyoi
void file(){
#ifdef zxyoi
	freopen("border.in","r",stdin);
#else
#ifndef ONLINE_JUDGE
	freopen("border.in","r",stdin);freopen("border.out","w",stdout);
#endif
#endif
}
signed main(){file();work();return 0;}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值