hdu5238 calculator (线段树+crt)

（并不能）发现29393不是质数，而是等于7*13*17*19

于是可以用四个线段树分别维护模意义下，对x进行一个区间的操作后的值

最后再把这四个的答案用crt拼起来

也可以不crt，而是预处理0~29392的每个情况

为了降低复杂度，预处理模7/13/17/19的幂

注意询问时，要把询问对7/13/17/19先取模

复杂度$O((7+13+17+19)nlogn)$

#include<cstdio>
#include<cstring>
#include<algorithm>
#define pa pair<int,int>
#define CLR(a,x) memset(a,x,sizeof(a))
using namespace std;
const int maxn=5e4+10;

inline int rd(){
    int x=0;char c=getchar();int neg=1;
    while(c<'0'||c>'9'){if(c=='-') neg=-1;c=getchar();}
    while(c>='0'&&c<='9') x=x*10+c-'0',c=getchar();
    return x*neg;
}

struct Node{
    int op,v;
};

int Pow[20][20][20],ans[10][15][20][20];

inline void init_pow(int x){
    for(int i=0;i<x;i++){
        Pow[x][i][0]=1;
        for(int j=1;j<x-1;j++) Pow[x][i][j]=Pow[x][i][j-1]*i%x;
    }
}

inline int fpow(int x,int y,int mod){
    return Pow[mod][x%mod][y%(mod-1)];
}

inline int cal(int x,Node y,int mod){
    if(y.op==1) return (x+y.v)%mod;
    else if(y.op==2) return x*y.v%mod;
    else return fpow(x,y.v,mod);
}

struct SegT{
    int trans[maxn<<2][20],mod;
    
    inline void update(int p){
        int a=p<<1,b=p<<1|1;
        for(int i=0;i<mod;i++){
            trans[p][i]=trans[b][trans[a][i]];
        }
    }
    inline void change(int p,int l,int r,int x,Node y){
        if(l==r){
            for(int i=0;i<mod;i++) trans[p][i]=cal(i,y,mod);
        }else{
            int m=(l+r)>>1;
            if(x<=m) change(p<<1,l,m,x,y);
            else change(p<<1|1,m+1,r,x,y);
            update(p);
        }
    }
}tr[4];

int N,M;

inline Node getnode(){
    char s[10];Node x;
    scanf("%s",s);
    if(s[0]=='+') x.op=1;
    else if(s[0]=='*') x.op=2;
    else x.op=3;
    x.v=0;
    for(int k=1,j=strlen(s);k<j;k++) x.v=x.v*10+s[k]-'0';
    return x;
}

inline int solve(int x){
    int a[5];
    for(int i=0;i<4;i++) a[i]=tr[i].trans[1][x%tr[i].mod];
    return ans[a[0]][a[1]][a[2]][a[3]];
}

int main(){
    int i,j,k;
    tr[0].mod=7,tr[1].mod=13,tr[2].mod=17,tr[3].mod=19;
    init_pow(7),init_pow(13),init_pow(17),init_pow(19);
    for(i=0;i<29393;i++){
        ans[i%7][i%13][i%17][i%19]=i;
    }
    for(int T=rd(),t=1;t<=T;t++){
        for(i=0;i<4;i++){
            for(j=1;j<=4*N;j++){
                for(k=0;k<tr[i].mod;k++)
                    tr[i].trans[j][k]=k;
            }
        }
        printf("Case #%d:
",t);
        N=rd(),M=rd();
        for(i=1;i<=N;i++){
            Node x=getnode();
            for(j=0;j<4;j++) tr[j].change(1,1,N,i,x);
        }
        for(i=1;i<=M;i++){
            int a=rd();
            if(a==1){
                printf("%d
",solve(rd()));
            }else{
                int x=rd();Node y=getnode();
                for(j=0;j<4;j++) tr[j].change(1,1,N,x,y);
            }
        }
    }
    return 0;
}