HDU 5238 线段树+数论

原题：http://acm.hdu.edu.cn/showproblem.php?pid=5238、

题解:给你长度为n的操作序列，和m组操作求每组操作的模29393的值。这道题直接显然是没有前途的，观察29393，有29393=7*13*17*19，显然可以分别处理这些小质数，在用CRT合并。处理方法可以用线段树。如设：

为第rt个质数x区间输入为i的答案，那么就有

这样就可以用线段树树进行单点修改。

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 50005
#define M 29393
#define L 5
using namespace std;
struct pa{int k,v;}a[N];
int cas,n,k;
int m[5]={0,7,13,17,19},b[5],ans[5];
char str[10];
int f[5][N<<2][20],T;
inline int ls(int x){return x<<1;} 
inline int rs(int x){return x<<1|1;}
inline pa make(int k,int v) {pa o; o.k=k;o.v=v;return o;}
int exgcd(int a,int b,int &x,int &y){
	if(a==0){
		x=0;y=1;return b;
	}
	int tx,ty,d;
	d=exgcd(b%a,a,tx,ty);
	x=ty-(b/a)*tx;
	y=tx; 
	return d;
}
inline int getnum(char *s){
	int x=0;
	for(int i=0;i<strlen(s);i++)
		if(s[i]>='0' && s[i]<='9') x=(x<<1)+(x<<3)+s[i]-'0'; 
	return x;
}
inline int check(char s){
	if(s=='+') return 1;
	if(s=='*') return 2;
	if(s=='^') return 3;
	return 0;
}
inline int calc(int key,int x,int y,int rt){
	if(key==0) return 0;
	if(key==1) return (x+y)%m[rt]; 
	if(key==2) return (x*y)%m[rt];
	if(key==3){
		int  a=x;int b=y;int ans=1;
		while(b){
			if(b&1) ans=ans*a%m[rt];
			b>>=1; a=a*a%m[rt];
		}
		return ans%m[rt];
	}
}
void build(int x,int l,int r,int rt){
	if(l==r) {
		for(int i=0;i<m[rt];i++) f[rt][x][i]=calc(a[l].k,i,a[l].v,rt);
		return ;
	}
	int mid=(l+r)>>1;
	build(ls(x),l,mid,rt);
	build(rs(x),mid+1,r,rt);
	for(int i=0;i<m[rt];i++) f[rt][x][i]=f[rt][rs(x)][f[rt][ls(x)][i]];
}
int change(int x,int l,int r,int key,int val,int pos,int rt){
	if(l==pos && r==pos)
		if(key==0) {
			return val;
		}
		else{
			for(int i=0;i<m[rt];i++) f[rt][x][i]=calc(key,i,val,rt);return -1;
		}
		

	int mid=(l+r)>>1;int lans=0;
	if(pos<=mid) lans=change(ls(x) ,l,mid,key,val,pos,rt );
	if(mid<pos)  change(rs(x) , mid+1,r,key,val,pos,rt);
	
	if(key!=0){
		for(int i=0;i<m[rt];i++) f[rt][x][i]=f[rt][rs(x)][f[rt][ls(x)][i]];return -1;
	}
 	return (f[rt][rs(x)][lans]);
}
inline int crt(){
	int s=0;
	for(int i=1;i<=4;i++)
		s=(s+b[i]*ans[i])%M; 
	return s;
}
int main(){
	//freopen("g.in","r",stdin);
	scanf("%d",&cas);
	for(int i=1,A,B,x,y;i<=4;i++){
		exgcd(M/m[i],m[i],x,y);
		b[i]=(M/m[i])*((x+m[i])%m[i]);
	}
	while(cas--){
		scanf("%d%d",&n,&k);
		a[1]=make(0,0);
		for(int i=1,key,val;i<=n;i++) {
			scanf("%s",&str);
			key=check(str[0]);val=getnum(str);
			a[i+1]=make(key,val);
		}
		
		n++;for(int i=1;i<L;i++) build(1,1,n,i);
		
		printf("Case #%d:\n",++T);
		for(int i=1,op,x,key,val;i<=k;i++){
			scanf("%d%d",&op,&x);
			if(op==1){
				for(int j=1;j<=4;j++){
					ans[j]=change(1,1,n,0,x%m[j],1,j);
				}
				printf("%d\n",crt());
			}else {
				scanf("%s",&str) ;
				key=check(str[0]);val=getnum(str);
				for(int j=1;j<=4;j++)
					change(1,1,n,key,val,x+1,j);
			}		
		}	
	}
	return 0;
}