hdu 4578 Transformation 线段树

没什么说的裸线段树，注意细节就好了！！！

代码如下：

 #include<iostream>
 #include<stdio.h>
 #include<algorithm>
 #include<iomanip>
 #include<cmath>
 #include<cstring>
 #include<vector>
 #define ll __int64
 #define pi acos(-1.0)
 #define MAX 100003
 #define M 10007
 #define lson i<<1
 #define rson i<<1|1
 using namespace std;
 struct tree
 {
     int l,r,sum[],add,mul;
     void init(int _l,int _r){
         l=_l;r=_r;mul=;add=;
     }
     void Add(int a){
         sum[]=(sum[]+a*a%M*a%M*len()%M+*a*sum[]%M+*a*a%M*sum[]%M)%M;
         sum[]=(sum[]+a*a%M*len()%M+*a*sum[]%M)%M;
         sum[]=(sum[]+a*len()%M)%M;
         add=(add+a)%M;
     }
     void Mul(int a){
         sum[]=sum[]*a%M*a%M*a%M;
         sum[]=sum[]*a%M*a%M;
         sum[]=sum[]*a%M;
         add=add*a%M;
         mul=mul*a%M;
     }
     int len(){
         return r-l+;
     }
     int mid(){
         return (l+r)>>;
     }
 }T[MAX*];
 void down(int i)
 {
     T[lson].Mul(T[i].mul);
     T[rson].Mul(T[i].mul);
     T[lson].Add(T[i].add);
     T[rson].Add(T[i].add);
     T[i].add=;
     T[i].mul=;
 }
 void updateson(int i)
 {
     for(int j=;j<;j++)
         T[i].sum[j]=(T[lson].sum[j]+T[rson].sum[j])%M;
 }
 void built(int i,int l,int r)
 {
     T[i].init(l,r);
     memset(T[i].sum,,sizeof(T[i].sum));
     if(l==r) return ;
     int m=T[i].mid();
     built(lson,l,m);
     built(rson,m+,r);
 }
 void update(int i,int l,int r,int mul,int add)
 {
     if(T[i].l==l&&T[i].r==r){
         T[i].Mul(mul);
         T[i].Add(add);
         return ;
     }
     down(i);
     int m=T[i].mid();
     if(r<=m) update(lson,l,r,mul,add);
     else if(l>m) update(rson,l,r,mul,add);
     else{
         update(lson,l,m,mul,add);
         update(rson,m+,r,mul,add);
     }
     updateson(i);
 }
 int query(int i,int l,int r,int p)
 {
     if(T[i].l==l&&T[i].r==r) return T[i].sum[p-];
     down(i);
     int m=T[i].mid();
     if(r<=m) return query(lson,l,r,p);
     else if(l>m) return query(rson,l,r,p);
     else return (query(lson,l,m,p)+query(rson,m+,r,p))%M;
 }
 int main(){
     int m,n,op,x,y,p;
     while(scanf("%d%d",&n,&m)&&(m+n)){
         built(,,n);
         while(m--){
             scanf("%d%d%d%d",&op,&x,&y,&p);
             if(op==) update(,x,y,,p);
             else if(op==) update(,x,y,p,);
             else if(op==) update(,x,y,,p);
             else printf("%d\n",query(,x,y,p)%M);
         }
     }
     return ;
 }