[CF878E]Numbers on the blackboard

Description

给出n个数字，每次询问一个区间[l,r]，对这个区间内部的点进行操作。
每次操作可以合并相邻两个数x,y，将它们变成x+2y
对于每次询问输出当最后只剩下一个数字时，这个数字的最大值。
询问互相独立，答案对1e9+7取模。
n<=1e5

Solution

首先我们考虑怎么O(n)解决询问
容易发现最后的答案一定是形如$\sum_{i=1}^{n}ai2^{ki}$的式子
其中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;
}