HDU 4578 线段树复杂题

题目大意：

题意：有一个序列，有四种操作：

1：区间[l，r]内的数全部加c。

2：区间[l，r]内的数全部乘c。

3：区间[l，r]内的数全部初始为c。

4：询问区间[l，r]内所有数的P次方之和。

这里p可以等于1,2,3三种情况，所以我们需要建立3个数组，当然这里其实只用一个sum[4*N][3]的2维数组其实更好

to[],add[],mul[]为三个懒惰标记to[]先于另外两个，每次做完to[]，另外两个标记就要将其初始化

对于mul[]来说，每次执行，add[cur]*=mul[cur]也要随之增加

因为这道题目数过大，需要求mod 10007的值

因为这个我把代码从上午看到了下午简直要疯了

int最大值2147483647

你必须每执行一步操作都需要mod一次，*2后，10007*2*10007 int还不至于爆掉，但也差不多

但是若*3之后那肯定要是再碰到一个超大数，那就GG了

所以每次乘法算完都mod一次吧，每次做mod带上括号是一个好习惯~~

 //这鬼题目里面最好每执行一次乘法就mod一次，不然可能int要爆掉导致报错，尤其是我在fun_add()中改错改了一天才发现是这个原因
 #include <cstdio>
 #include <cstring>
 using namespace std;
 #define N 100005
 #define mod 10007
 int sum[*N][],add[*N],to[*N],mul[*N];
 void update(int cur)
 {
     sum[cur][]=(sum[cur<<][]+sum[cur<<|][])%mod;
     sum[cur][]=(sum[cur<<][]+sum[cur<<|][])%mod;
     sum[cur][]=(sum[cur<<][]+sum[cur<<|][])%mod;
 }
 void build(int cur,int x,int y)
 {
     int mid=(x+y)/,ls=cur<<,rs=cur<<|;
     add[cur]=,to[cur]=,mul[cur]=;
     if(x==y){
         sum[cur][]=sum[cur][]=sum[cur][]=;
         return;
     }
     build(ls,x,mid);
     build(rs,mid+,y);
     update(cur);
 }
 void fun_mul(int cur,int x,int y,int val)
 {
     mul[cur]*=val;
     mul[cur]%=mod;
     add[cur]*=val;
     add[cur]%=mod;

     int t1=val*val%mod;
     int t2=t1*val%mod;
     sum[cur][]=sum[cur][]*val;
     sum[cur][]%=mod;

     sum[cur][]=sum[cur][]*t1;
     sum[cur][]%=mod;

     sum[cur][]=sum[cur][]*t2;
     sum[cur][]%=mod;
 }
 void fun_add(int cur,int x,int y,int val)
 {
     add[cur]+=val;
     add[cur]%=mod;
     int t1=sum[cur][];
     int t2=sum[cur][];
     int tmp1=val*val%mod;
     int tmp2=val*tmp1%mod;
     sum[cur][]+=(y-x+)*val;
     sum[cur][]%=mod;

     sum[cur][]+=(y-x+)*tmp1%mod+((*val)%mod)*t1%mod;
     sum[cur][]%=mod;

     sum[cur][]+=((*val)%mod)*t2%mod+*tmp1*t1%mod+(y-x+)*tmp2%mod;
     sum[cur][]%=mod;
 }
 void fun_to(int cur,int x,int y,int val)
 {
     mul[cur]=,add[cur]=,to[cur]=val%mod;
     int t1=val*val%mod;
     int t2=t1*val%mod;
     sum[cur][]=(y-x+)*val%mod;
     sum[cur][]=(y-x+)*t1%mod;
     sum[cur][]=(y-x+)*t2%mod;
 }
 void pushdown(int cur,int x,int y)
 {
     int mid=(x+y)/,ls=cur<<,rs=cur<<|;
     if(to[cur]){
         //to[ls]=to[rs]=to[cur];
         fun_to(ls,x,mid,to[cur]);
         fun_to(rs,mid+,y,to[cur]);
         to[cur]=;
     }
     if(mul[cur]>){
         //mul[ls]*=mul[cur],mul[rs]*=mul[cur];
         fun_mul(ls,x,mid,mul[cur]);
         fun_mul(rs,mid+,y,mul[cur]);
         mul[cur]=;
     }
     if(add[cur]){
         //add[ls]+=add[cur],add[rs]+=add[cur];
         fun_add(ls,x,mid,add[cur]);
         fun_add(rs,mid+,y,add[cur]);
         add[cur]=;
     }
 }
 void change(int cur,int x,int y,int s,int t,int op,int v)
 {
     int mid=(x+y)/,ls=cur<<,rs=cur<<|;
     if(x>=s&&y<=t){
         if(op==) fun_add(cur,x,y,v);
         if(op==) fun_mul(cur,x,y,v);
         if(op==) fun_to(cur,x,y,v);
         return;
     }
     pushdown(cur,x,y);
     if(mid>=s) change(ls,x,mid,s,t,op,v);
     if(mid+<=t) change(rs,mid+,y,s,t,op,v);
     update(cur);
 }
 void query(int cur,int x,int y,int s ,int t,int p,int &ans)
 {
     int mid=(x+y)/,ls=cur<<,rs=cur<<|;
     if(x>=s&&y<=t){
         ans+=sum[cur][p-];
         ans%=mod;
         return;
     }
     pushdown(cur,x,y);
     if(mid>=s) query(ls,x,mid,s,t,p,ans);
     if(mid+<=t) query(rs,mid+,y,s,t,p,ans);
 }
 int main()
 {
     int n,m,op,x,y,p;
     while(scanf("%d%d",&n,&m),n+m){
         build(,,n);
         for(int i=;i<m;i++){
             scanf("%d%d%d%d",&op,&x,&y,&p);
             if(op==){
                 int ans=;
                 //printf("%d\n",0);
                 query(,,n,x,y,p,ans);
                 printf("%d\n",ans%mod);
             }
             else change(,,n,x,y,op,p);
         }
     }
 }