「BJOI2018」链上二次求和

「BJOI2018」链上二次求和

https://loj.ac/problem/2512

我说今天上午写博客吧。怕自己写一上午，就决定先写道题。

然后我就调了一上午线段树。

花了2h找到lazy标记没有清空。我tm清空了有没有标记没清空标记本身。

又花25min找到某个乘法爆int了。int真的淡疼，要不是longlong自带巨无霸常数，这辈子都不想用int。

一个上午就没有了。

 //Achen
 #include<bits/stdc++.h>
 #define For(i,a,b) for(int i=(a);i<=(b);i++)
 #define Rep(i,a,b) for(int i=(a);i>=(b);i--)
 #define Formylove return 0
 const int N=,p=1e9+,inv2=5e8+;
 typedef long long LL;
 typedef double db;
 using namespace std;
 int n,m,a[N],sum[N];

 template<typename T> void read(T &x) {
     char ch=getchar(); T f=; x=;
     while(ch!='-'&&(ch<''||ch>'')) ch=getchar();
     if(ch=='-') f=-,ch=getchar();
     for(;ch>=''&&ch<='';ch=getchar()) x=x*+ch-''; x*=f;
 }

 int mo(int x) { return x<?x+p:(x>=p?x-p:x); }

 int sg[N<<],si0[N<<],si1[N<<],si2[N<<];
 bool tg[N<<];
 struct tag {
     int a,b,c;
 }lz[N<<],pls;
 tag operator +(const tag&A,const tag&B) { return (tag){mo(A.a+B.a),mo(A.b+B.b),mo(A.c+B.c)}; }

 #define lc (x<<1)
 #define rc ((x<<1)|1)
 #define mid ((l+r)>>1)
 void addtag(int x,tag y) {
     //printf("%d %d %d %d\n",x,y.a,y.b,y.c);
     sg[x]=mo(mo(mo(sg[x]+(LL)y.a*si2[x]%p)+(LL)y.b*si1[x]%p)+(LL)y.c*si0[x]%p);
     lz[x]=lz[x]+y; tg[x]=;
     //printf("%d %d %d\n",lz[x].a,lz[x].b,lz[x].c);
 }

 void down(int x) {
     if(!tg[x]) return ;
     addtag(lc,lz[x]);
     addtag(rc,lz[x]);
     tg[x]=; lz[x]=(tag){,,};/////////
 }

 void build(int x,int l,int r) {
     if(l==r) {
         sg[x]=sum[l]; si0[x]=; si1[x]=l;
         si2[x]=(LL)l*l%p; return;
     }
     build(lc,l,mid); build(rc,mid+,r);
     sg[x]=mo(sg[lc]+sg[rc]);
     si0[x]=mo(si0[lc]+si0[rc]);
     si1[x]=mo(si1[lc]+si1[rc]);
     si2[x]=mo(si2[lc]+si2[rc]);
 }

 void upd(int x,int l,int r,int ql,int qr) {
     if(l>=ql&&r<=qr) { addtag(x,pls); return; }
     down(x);
     if(ql<=mid) upd(lc,l,mid,ql,qr);
     if(qr>mid) upd(rc,mid+,r,ql,qr);
     sg[x]=mo(sg[lc]+sg[rc]);
 }

 int qry(int x,int l,int r,int ql,int qr) {
     if(ql>qr) return ;
     if(l>=ql&&r<=qr) return sg[x];
     down(x);
     if(qr<=mid) return qry(lc,l,mid,ql,qr);
     if(ql>mid) return qry(rc,mid+,r,ql,qr);
     return mo(qry(lc,l,mid,ql,qr)+qry(rc,mid+,r,ql,qr));
 }

 int main() {
     //freopen("1.in","r",stdin);
     //freopen("1.out","w",stdout);
     read(n); read(m);
     For(i,,n) { read(a[i]); sum[i]=mo(sum[i-]+a[i]); }
     For(i,,n) sum[i]=mo(sum[i-]+sum[i]);
     build(,,n);
     For(i,,m) {
         int o,l,r,v;
         read(o); read(l); read(r);
         if(o==) {
             read(v); if(l>r) swap(l,r);
             int len=r-l+,a=(LL)v*inv2%p;
             pls=(tag){a,(LL)a*(-*l+p)%p,((LL)l*l-*l+)%p*a%p}; upd(,,n,l,r);
             if(r<n) {
                 pls=(tag){,(LL)len*v%p,((LL)len*(len+)/%p*v%p-(LL)len*r%p*v%p+p)%p};
                 upd(,,n,r+,n);
             }
         }
         else if(o==) {
             LL ans=(((LL)r-l+)*qry(,,n,n,n)%p-qry(,,n,max(,l-),r-)+p-qry(,,n,max(,n-r),n-l)+p)%p;
             printf("%lld\n",ans);
         }
     }
     Formylove;
 }