线段树毒瘤题——序列操作（清华集训）

蒟蒻的垂死挣扎

这道题毒瘤之处在于，你要暴力拆一波式子。

反正最终你会得到一个优秀的式子，然后合并贡献直接平方暴力合并即可。

 

新博客里会讲的很清楚的（吧）。

 

（没时间解释了，放学了，后面会去维护新博客咯！）

 

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
#include<map>
#include<vector>
#include<queue>
#include<ctime>
using namespace std;
#define ll long long
#define N 50100
#define MOD 19940417
#define ls (now<<1)
#define rs ((now<<1)|1)
#define RG register
inline ll read(){
    RG ll x=0,t=1;RG char ch=getchar();
    while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
    if(ch=='-')t=-1,ch=getchar();
    while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
    return x*t;
}

int n,m,val[N],C[N][30],t[30];
struct SeqTree{
	int sum[N<<2][21],lz1[N<<2],lz2[N<<2],ans[21],tmp[21];
	inline void Pushup(int now){
		for(RG int i=0;i<=20;++i) sum[now][i]=0;
		for(RG int i=0,k=20;i<=20;++i,--k)
			for(RG int j=0;j<=k;++j)
				sum[now][i+j]+=1LL*sum[ls][i]*sum[rs][j]%MOD;
		for(RG int i=1;i<=20;++i) sum[now][i]%=MOD;
	}
	inline void Pushlazy(int now,int plu,int mul,int siz){
		int to=20;
		if(mul!=1){
			for(RG int i=1;i<=to;i+=2) sum[now][i]*=-1;
			lz1[now]*=-1,lz2[now]*=-1;
		}
		if(plu){   
			for(RG int i=0;i<=to;++i) tmp[i]=sum[now][i];
			for(RG int i=1;i<=to;++i) t[i]=1LL*t[i-1]*plu%MOD;
			for(RG int i=0;i<=to;++i)
				for(RG int j=0;j<i;++j)
					sum[now][i]=(1LL*t[i-j]*tmp[j]%MOD*C[siz-j][i-j]+sum[now][i])%MOD;
			for(RG int i=0;i<=20;++i) tmp[i]=0;	lz1[now]+=plu;
		}   return ;
	}
	inline void Pushdown(int now,int siz){
		if(!lz1[now]&&lz2[now]==1) return ; int Siz=siz>>1;
		Pushlazy(ls,lz1[now],lz2[now],siz-Siz),Pushlazy(rs,lz1[now],lz2[now],Siz);
		lz1[now]=0,lz2[now]=1;
	}
	inline void Build(int l,int r,int now){
		sum[now][0]=lz2[now]=1;
		if(l==r) { sum[now][1]=val[l]; return ; }
		int mid=l+r>>1;
		Build(l,mid,ls),Build(mid+1,r,rs);
		Pushup(now);
	}
	inline void Modify(int l,int r,int now,int L,int R,int plu,int mul){
		if(l>=L&&r<=R) { Pushlazy(now,plu,mul,r-l+1); return ; }
		int mid=l+r>>1; Pushdown(now,r-l+1);
		if(L<=mid) Modify(l,mid,ls,L,R,plu,mul);
		if(R>mid) Modify(mid+1,r,rs,L,R,plu,mul);
		Pushup(now);
	}
	inline void Query(int l,int r,int now,int L,int R){
		if(l>=L&&r<=R){
			for(RG int i=0,k=20;i<=20;++i,--k)
				for(RG int j=0;j<=k;++j)
					tmp[i+j]+=1LL*sum[now][i]*ans[j]%MOD;
			for(RG int i=0;i<=20;++i) ans[i]=tmp[i]%MOD,tmp[i]=0;
			return ;
		}   int mid=l+r>>1; Pushdown(now,r-l+1);
		if(L<=mid) Query(l,mid,ls,L,R);
		if(R>mid) Query(mid+1,r,rs,L,R);
		Pushup(now);
	}
} T;

int main(){
	n=read(),m=read();
	for(RG int i=1;i<=n;++i) val[i]=read()%MOD;
	T.Build(1,n,1),t[0]=1;
	for(RG int i=0;i<=n;++i) C[i][0]=1;
	for(RG int i=1;i<=n;++i)
		for(RG int j=1;j<=min(i,20);++j)
			C[i][j]=(C[i-1][j]+C[i-1][j-1])%MOD;
	for(RG int i=1;i<=m;++i){
		char op[5]; scanf("%s",op);
		if(op[0]=='I'){
			int l=read(),r=read(),w=read();
			T.Modify(1,n,1,l,r,w,1);
		}
		if(op[0]=='Q'){
			int l=read(),r=read(),w=read();
			for(RG int i=0;i<=20;++i) T.ans[i]=0;
			T.ans[0]=1,T.Query(1,n,1,l,r);
			printf("%d\n",(T.ans[w]+MOD)%MOD);
		}
		if(op[0]=='R'){
			int l=read(),r=read();
			T.Modify(1,n,1,l,r,0,-1);
		}
	}
}