hdu 4578 Transformation 线段树
没什么说的裸线段树，注意细节就好了！！！
代码如下：
1 #include<iostream> 2 #include<stdio.h> 3 #include<algorithm> 4 #include<iomanip> 5 #include<cmath> 6 #include<cstring> 7 #include<vector> 8 #define ll __int64 9 #define pi acos(-1.0) 10 #define MAX 100003 11 #define M 10007 12 #define lson i<<1 13 #define rson i<<1|1 14 using namespace std; 15 struct tree 16 { 17 int l,r,sum[3],add,mul; 18 void init(int _l,int _r){ 19 l=_l;r=_r;mul=1;add=0; 20 } 21 void Add(int a){ 22 sum[2]=(sum[2]+a*a%M*a%M*len()%M+3*a*sum[1]%M+3*a*a%M*sum[0]%M)%M; 23 sum[1]=(sum[1]+a*a%M*len()%M+2*a*sum[0]%M)%M; 24 sum[0]=(sum[0]+a*len()%M)%M; 25 add=(add+a)%M; 26 } 27 void Mul(int a){ 28 sum[2]=sum[2]*a%M*a%M*a%M; 29 sum[1]=sum[1]*a%M*a%M; 30 sum[0]=sum[0]*a%M; 31 add=add*a%M; 32 mul=mul*a%M; 33 } 34 int len(){ 35 return r-l+1; 36 } 37 int mid(){ 38 return (l+r)>>1; 39 } 40 }T[MAX*4]; 41 void down(int i) 42 { 43 T[lson].Mul(T[i].mul); 44 T[rson].Mul(T[i].mul); 45 T[lson].Add(T[i].add); 46 T[rson].Add(T[i].add); 47 T[i].add=0; 48 T[i].mul=1; 49 } 50 void updateson(int i) 51 { 52 for(int j=0;j<3;j++) 53 T[i].sum[j]=(T[lson].sum[j]+T[rson].sum[j])%M; 54 } 55 void built(int i,int l,int r) 56 { 57 T[i].init(l,r); 58 memset(T[i].sum,0,sizeof(T[i].sum)); 59 if(l==r) return ; 60 int m=T[i].mid(); 61 built(lson,l,m); 62 built(rson,m+1,r); 63 } 64 void update(int i,int l,int r,int mul,int add) 65 { 66 if(T[i].l==l&&T[i].r==r){ 67 T[i].Mul(mul); 68 T[i].Add(add); 69 return ; 70 } 71 down(i); 72 int m=T[i].mid(); 73 if(r<=m) update(lson,l,r,mul,add); 74 else if(l>m) update(rson,l,r,mul,add); 75 else{ 76 update(lson,l,m,mul,add); 77 update(rson,m+1,r,mul,add); 78 } 79 updateson(i); 80 } 81 int query(int i,int l,int r,int p) 82 { 83 if(T[i].l==l&&T[i].r==r) return T[i].sum[p-1]; 84 down(i); 85 int m=T[i].mid(); 86 if(r<=m) return query(lson,l,r,p); 87 else if(l>m) return query(rson,l,r,p); 88 else return (query(lson,l,m,p)+query(rson,m+1,r,p))%M; 89 } 90 int main(){ 91 int m,n,op,x,y,p; 92 while(scanf("%d%d",&n,&m)&&(m+n)){ 93 built(1,1,n); 94 while(m--){ 95 scanf("%d%d%d%d",&op,&x,&y,&p); 96 if(op==1) update(1,x,y,1,p); 97 else if(op==2) update(1,x,y,p,0); 98 else if(op==3) update(1,x,y,0,p); 99 else printf("%d\n",query(1,x,y,p)%M); 100 } 101 } 102 return 0; 103 }
