Description

　　有一個長度為n的序列，有三個操作1.I a b c表示將[a,b]這一段區間的元素集體增加c，2.R a b表示將[a,b]區間內所有元素變成相反數，3.Q a b c表示詢問[a,b]這一段區間中選擇c個數相乘的所有方案的和mod 19940417的值。

Input

　　第一行兩個數n，q表示序列長度和操作個數。
　　第二行n個非負整數，表示序列。
　　接下來q行每行輸入一個操作I a b c或者 R a b或者Q a b c意義如題目描述。

Output

　　對於每個詢問，輸出選出c個數相乘的所有方案的和mod19940417的值。

Sample Input

Sample Output

HINT

　　100%的數據n<=50000,q<=50000,初始序列的元素的絕對值<=109，I a b c中保證[a,b]是一個合法區間，|c|<=109，R a b保證[a,b]是個合法的區間。Q a b c中保證[a,b]是個合法的區間1<=c<=min(b-a+1,20)。

Source

中國國家隊清華集訓 2012-2013 第三天

#include<bits/stdc++.h>
#define ll long long
#define maxn 50005
#define ha 19940417
using namespace std;
int C[maxn][25],n,m,tp,num;
int a[maxn],le,ri,w,ci[25];

inline int add(int x,int y){
    x+=y;
    return x>=ha?x-ha:x;
}

inline void init(){
    C[0][0]=1;
    for(int i=1;i<=n;i++){
        tp=min(i,20);
        C[i][0]=1;
        for(int j=1;j<=tp;j++) C[i][j]=add(C[i-1][j-1],C[i-1][j]);
    }
}

struct node{
    int tag,ad;
    int c[25];
        
    inline void clear(){
        memset(c,0,sizeof(c));
        tag=ad=0;
    }
    
    node operator +(const node &U)const{
        node r;
        r.clear();
        for(int i=0;i<=20;i++) if(c[i])
            for(int j=i;j<=20;j++) if(U.c[j-i])
                r.c[j]=add(r.c[j],c[i]*(ll)U.c[j-i]%ha);
        
        return r;
    }
}b[maxn<<2|1],AN;
char ch;

void build(int o,int l,int r){
    if(l==r){
        b[o].c[0]=1;
        b[o].c[1]=a[l];
        return;
    }
    
    int mid=l+r>>1,lc=o<<1,rc=(o<<1)|1;
    build(lc,l,mid);
    build(rc,mid+1,r);
    
    b[o]=b[lc]+b[rc];
}

inline void work(int o,int d,int len){
    tp=min(len,20);
    ci[0]=1;
    for(int i=1;i<=tp;i++) ci[i]=ci[i-1]*(ll)d%ha;
    b[o].ad=add(b[o].ad,d);
    
    for(int i=20;i;i--){
        int S=i,T=len;
        for(int j=0;j<i;j++,S--,T--){
            b[o].c[i]=add(b[o].c[i],b[o].c[j]*(ll)ci[i-j]%ha*(ll)C[T][S]%ha);
        }
    }
}

inline void rotate(int o,int len){
    b[o].tag^=1,b[o].ad=ha-b[o].ad;
    for(int i=1;i<=20;i+=2) b[o].c[i]=ha-b[o].c[i];
}

inline void pushdown(int o,int l,int r){
    int mid=l+r>>1,lc=o<<1,rc=(o<<1)|1;
    if(b[o].tag){
        rotate(lc,mid-l+1);
        rotate(rc,r-mid);
        b[o].tag=0;
    }
    if(b[o].ad){
        work(lc,b[o].ad,mid-l+1);
        work(rc,b[o].ad,r-mid);
        b[o].ad=0;
    }
}

void updateX(int o,int l,int r){
    if(l>=le&&r<=ri){
        work(o,w,r-l+1);		
        return;
    }
    
    int mid=l+r>>1,lc=o<<1,rc=(o<<1)|1;
    pushdown(o,l,r);
    if(le<=mid) updateX(lc,l,mid);
    if(ri>mid) updateX(rc,mid+1,r);
        
    b[o]=b[lc]+b[rc];
}

void updateY(int o,int l,int r){
    if(l>=le&&r<=ri){
        rotate(o,r-l+1);
        return;
    }
    
    int mid=l+r>>1,lc=o<<1,rc=(o<<1)|1;
    pushdown(o,l,r);
    if(le<=mid) updateY(lc,l,mid);
    if(ri>mid) updateY(rc,mid+1,r);
    
    b[o]=b[lc]+b[rc];
}

void query(int o,int l,int r){
    if(l>=le&&r<=ri){
        AN=AN+b[o];
        return;
    }
    
    int mid=l+r>>1,lc=o<<1,rc=(o<<1)|1;
    pushdown(o,l,r);
    if(le<=mid) query(lc,l,mid);
    if(ri>mid) query(rc,mid+1,r);
}

int main(){
    scanf("%d%d",&n,&m);
    init();
    for(int i=1;i<=n;i++){
        scanf("%d",a+i);
        a[i]%=ha;
        if(a[i]<0) a[i]+=ha;
    }
    build(1,1,n);
    
    while(m--){
        ch=getchar();
        while(ch!=‘I‘&&ch!=‘R‘&&ch!=‘Q‘) ch=getchar();
        scanf("%d%d",&le,&ri);
        if(ch==‘R‘) updateY(1,1,n);
        else{
            scanf("%d",&w);
            if(ch==‘Q‘){
                AN.clear();
                AN.c[0]=1;
                query(1,1,n);
                printf("%d\n",AN.c[w]);
            }
            else{
                w%=ha;
                if(w<0) w+=ha;
                updateX(1,1,n);
            }
        }
    }
    
    return 0;
}

bzoj 2962: 序列操作