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+;

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

 struct Node{
     int op,v;
 };

 int Pow[][][],ans[][][][];

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

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

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

 struct SegT{
     int trans[maxn<<][],mod;

     inline void update(int p){
         int a=p<<,b=p<<|;
         for(int i=;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=;i<mod;i++) trans[p][i]=cal(i,y,mod);
         }else{
             int m=(l+r)>>;
             if(x<=m) change(p<<,l,m,x,y);
             else change(p<<|,m+,r,x,y);
             update(p);
         }
     }
 }tr[];

 int N,M;

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

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

 int main(){
     int i,j,k;
     tr[].mod=,tr[].mod=,tr[].mod=,tr[].mod=;
     init_pow(),init_pow(),init_pow(),init_pow();
     for(i=;i<;i++){
         ans[i%][i%][i%][i%]=i;
     }
     for(int T=rd(),t=;t<=T;t++){
         for(i=;i<;i++){
             for(j=;j<=*N;j++){
                 for(k=;k<tr[i].mod;k++)
                     tr[i].trans[j][k]=k;
             }
         }
         printf("Case #%d:\n",t);
         N=rd(),M=rd();
         for(i=;i<=N;i++){
             Node x=getnode();
             for(j=;j<;j++) tr[j].change(,,N,i,x);
         }
         for(i=;i<=M;i++){
             int a=rd();
             if(a==){
                 printf("%d\n",solve(rd()));
             }else{
                 int x=rd();Node y=getnode();
                 for(j=;j<;j++) tr[j].change(,,N,x,y);
             }
         }
     }
     return ;
 }