传送门
怎么说呢。。。这个贪心,难以形容。。。
题解:
首先考虑全为正数的时候,显然就是贪心从尾部开始向前合并,全为负数的时候,显然就是贪心每次取头部两个合并。
其实大胆猜想+证明可以得到这道题正确的贪心方式。
注意到贡献的表示形如 ∑ 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;}