Description
给出n个数字,每次询问一个区间[l,r],对这个区间内部的点进行操作。
每次操作可以合并相邻两个数x,y,将它们变成x+2y
对于每次询问输出当最后只剩下一个数字时,这个数字的最大值。
询问互相独立,答案对1e9+7取模。
n<=1e5
Solution
首先我们考虑怎么O(n)解决询问
容易发现最后的答案一定是形如
∑ni=1ai2ki
的式子
其中k[1]=0,k[i+1]<=k[i]+1
那么贪心的发现,我们的k会形成一段段连续的块,每块内k[i+1]=k[i]+1并且块头的k=1(除了第一块)
接下来我们考虑如何维护这个块,每次我们在末尾加入一个数,如果它是负数,那么自成一块,否则和前面那个块合并。
注意我们合并的时候可以引发一连串的合并。
接下来我们考虑多组询问,我们把询问离线按右端点排序,每次询问用并查集/二分找到左端点所属的块然后统计答案即可。
注意虽然答案要取模但我们有必要保留负数,注意判断overflow的情况。
Code
#include <cmath>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
typedef long long ll;
typedef double db;
int read() {
char ch;int sig=1;
for(ch=getchar();ch<'0'||ch>'9';ch=getchar()) if (ch=='-') sig=-1;
int x=ch-'0';
for(ch=getchar();ch>='0'&&ch<='9';ch=getchar()) x=x*10+ch-'0';
return x*sig;
}
void write(ll x) {
if (!x) {puts("0");return;}
if (x<0) {putchar('-');x=-x;}
char ch[20];int tot=0;
for(;x;x/=10) ch[++tot]=x%10;
fd(i,tot,1) putchar(ch[i]+'0');
puts("");
}
const int N=1e5+5,Mo=1e9+7;
struct que{int l,r,id;}ask[N];
bool cmp(que x,que y) {return x.r<y.r;}
int pwr(int x,int y) {
int z=1;
for(;y;y>>=1,x=(ll)x*x%Mo)
if (y&1) z=(ll)z*x%Mo;
return z;
}
int n,q,pre[N],fa[N],a[N];
ll sum[N],S[N],an[N],suf[N],two[N];
int get(int x) {return (fa[x]!=x)?fa[x]=get(fa[x]):x;}
void merge(int x,int y) {
fa[x]=fa[y];pre[y]=pre[x];
int len=x-pre[x];
if (len>30&&sum[y]>0||sum[x]+(sum[y]<<len)>Mo) sum[y]=Mo;
else sum[y]=sum[x]+(sum[y]<<len);
}
ll query(int l,int r) {
return (suf[l]-suf[r+1]*two[r-l+1]%Mo+Mo)%Mo;
}
int main() {
n=read();q=read();
fo(i,1,n) a[i]=read();
fo(i,1,q) ask[i].l=read(),ask[i].r=read(),ask[i].id=i;
sort(ask+1,ask+q+1,cmp);
fo(i,1,n) {
fa[i]=i;pre[i]=i-1;
sum[i]=a[i];
}
fd(i,n,1) suf[i]=(suf[i+1]*2+a[i]+Mo)%Mo;
two[0]=1;fo(i,1,n) two[i]=two[i-1]*2%Mo;
int j=0;
fo(i,1,n) {
while (pre[i]&&sum[i]>=0) merge(pre[i],i);
S[i]=(S[pre[i]]+query(pre[i]+1,i)*2)%Mo;
while (ask[j+1].r==i) {
int x=get(ask[++j].l);
an[ask[j].id]=(S[i]-S[x]+query(ask[j].l,x)+Mo)%Mo;
}
}
fo(i,1,q) write(an[i]);
return 0;
}