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带上括号是一个好习惯~~

  1 //这鬼题目里面最好每执行一次乘法就mod一次，不然可能int要爆掉导致报错，尤其是我在fun_add()中改错改了一天才发现是这个原因
  2 #include <cstdio>
  3 #include <cstring>
  4 using namespace std;
  5 #define N 100005
  6 #define mod 10007
  7 int sum[4*N][3],add[4*N],to[4*N],mul[4*N];
  8 void update(int cur)
  9 {
 10     sum[cur][0]=(sum[cur<<1][0]+sum[cur<<1|1][0])%mod;
 11     sum[cur][1]=(sum[cur<<1][1]+sum[cur<<1|1][1])%mod;
 12     sum[cur][2]=(sum[cur<<1][2]+sum[cur<<1|1][2])%mod;
 13 }
 14 void build(int cur,int x,int y)
 15 {
 16     int mid=(x+y)/2,ls=cur<<1,rs=cur<<1|1;
 17     add[cur]=0,to[cur]=0,mul[cur]=1;
 18     if(x==y){
 19         sum[cur][0]=sum[cur][1]=sum[cur][2]=0;
 20         return;
 21     }
 22     build(ls,x,mid);
 23     build(rs,mid+1,y);
 24     update(cur);
 25 }
 26 void fun_mul(int cur,int x,int y,int val)
 27 {
 28     mul[cur]*=val;
 29     mul[cur]%=mod;
 30     add[cur]*=val;
 31     add[cur]%=mod;
 32 
 33     int t1=val*val%mod;
 34     int t2=t1*val%mod;
 35     sum[cur][0]=sum[cur][0]*val;
 36     sum[cur][0]%=mod;
 37 
 38     sum[cur][1]=sum[cur][1]*t1;
 39     sum[cur][1]%=mod;
 40 
 41     sum[cur][2]=sum[cur][2]*t2;
 42     sum[cur][2]%=mod;
 43 }
 44 void fun_add(int cur,int x,int y,int val)
 45 {
 46     add[cur]+=val;
 47     add[cur]%=mod;
 48     int t1=sum[cur][0];
 49     int t2=sum[cur][1];
 50     int tmp1=val*val%mod;
 51     int tmp2=val*tmp1%mod;
 52     sum[cur][0]+=(y-x+1)*val;
 53     sum[cur][0]%=mod;
 54 
 55     sum[cur][1]+=(y-x+1)*tmp1%mod+((2*val)%mod)*t1%mod;
 56     sum[cur][1]%=mod;
 57 
 58     sum[cur][2]+=((3*val)%mod)*t2%mod+3*tmp1*t1%mod+(y-x+1)*tmp2%mod;
 59     sum[cur][2]%=mod;
 60 }
 61 void fun_to(int cur,int x,int y,int val)
 62 {
 63     mul[cur]=1,add[cur]=0,to[cur]=val%mod;
 64     int t1=val*val%mod;
 65     int t2=t1*val%mod;
 66     sum[cur][0]=(y-x+1)*val%mod;
 67     sum[cur][1]=(y-x+1)*t1%mod;
 68     sum[cur][2]=(y-x+1)*t2%mod;
 69 }
 70 void pushdown(int cur,int x,int y)
 71 {
 72     int mid=(x+y)/2,ls=cur<<1,rs=cur<<1|1;
 73     if(to[cur]){
 74         //to[ls]=to[rs]=to[cur];
 75         fun_to(ls,x,mid,to[cur]);
 76         fun_to(rs,mid+1,y,to[cur]);
 77         to[cur]=0;
 78     }
 79     if(mul[cur]>1){
 80         //mul[ls]*=mul[cur],mul[rs]*=mul[cur];
 81         fun_mul(ls,x,mid,mul[cur]);
 82         fun_mul(rs,mid+1,y,mul[cur]);
 83         mul[cur]=1;
 84     }
 85     if(add[cur]){
 86         //add[ls]+=add[cur],add[rs]+=add[cur];
 87         fun_add(ls,x,mid,add[cur]);
 88         fun_add(rs,mid+1,y,add[cur]);
 89         add[cur]=0;
 90     }
 91 }
 92 void change(int cur,int x,int y,int s,int t,int op,int v)
 93 {
 94     int mid=(x+y)/2,ls=cur<<1,rs=cur<<1|1;
 95     if(x>=s&&y<=t){
 96         if(op==1) fun_add(cur,x,y,v);
 97         if(op==2) fun_mul(cur,x,y,v);
 98         if(op==3) fun_to(cur,x,y,v);
 99         return;
100     }
101     pushdown(cur,x,y);
102     if(mid>=s) change(ls,x,mid,s,t,op,v);
103     if(mid+1<=t) change(rs,mid+1,y,s,t,op,v);
104     update(cur);
105 }
106 void query(int cur,int x,int y,int s ,int t,int p,int &ans)
107 {
108     int mid=(x+y)/2,ls=cur<<1,rs=cur<<1|1;
109     if(x>=s&&y<=t){
110         ans+=sum[cur][p-1];
111         ans%=mod;
112         return;
113     }
114     pushdown(cur,x,y);
115     if(mid>=s) query(ls,x,mid,s,t,p,ans);
116     if(mid+1<=t) query(rs,mid+1,y,s,t,p,ans);
117 }
118 int main()
119 {
120     int n,m,op,x,y,p;
121     while(scanf("%d%d",&n,&m),n+m){
122         build(1,1,n);
123         for(int i=0;i<m;i++){
124             scanf("%d%d%d%d",&op,&x,&y,&p);
125             if(op==4){
126                 int ans=0;
127                 //printf("%d\n",0);
128                 query(1,1,n,x,y,p,ans);
129                 printf("%d\n",ans%mod);
130             }
131             else change(1,1,n,x,y,op,p);
132         }
133     }
134 }

 

 posted on 2014-08-07 14:21  Love风吟  阅读(199)  评论(0)  编辑  收藏  举报