hdu 4578 Transformation

         比较裸的线段树，但是写起来比较烦，昨天无限wrong+TLE。后来发现写的实在太麻烦，只好推倒重写……

         这题主要问题就是处理乘和加的操作，其实这两个操作可以看成一个，乘看成a*c+0，加看成a*1+c，这样的话就可以把两个操作当成一个操作了，这么写应该是没问题的，但是我写挂了Orz……后来又换了下方法，每次都先处理乘的操作，再处理加的操作，但是先乘再加和先加再乘是不同的，如果已经有加的操作，比如是（a*k+b），这种情况再乘就变成（a*k+b）*c=a*k*c+b*c了，这样的话乘的时候再处理一下加的操作就ok了，另外，set操作可以看成*0+c。

 

代码：

 

#include <iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<map>
#include<queue>
#include<stack>
#include<cmath>
#include<vector>
#define inf 2139062143
#define Inf 0x3FFFFFFFFFFFFFFFLL
#define eps 1e-9
#define pi acos(-1.0)
using namespace std;
typedef long long ll;
const int maxn=100000+10;
const int mod=10007;
ll sum[maxn<<2],sum2[maxn<<2],sum3[maxn<<2],mulv[maxn<<2],addv[maxn<<2];
ll tmp1,tmp2;
inline void PushUp(int rt)
{
    sum[rt]=(sum[rt<<1]+sum[rt<<1|1])%mod;
    sum2[rt]=(sum2[rt<<1]+sum2[rt<<1|1])%mod;
    sum3[rt]=(sum3[rt<<1]+sum3[rt<<1|1])%mod;
}
inline void sumcal(int p,ll ad,int len)
{
    tmp1=sum[p],tmp2=sum2[p];
    ll v=(ad*ad*ad)%mod;
    sum[p]=(tmp1+len*ad)%mod;
    sum2[p]=(tmp2+len*ad*ad+2*ad*tmp1)%mod;
    sum3[p]=(sum3[p]+len*v+3*tmp2*ad+3*ad*ad*tmp1)%mod;
}
inline void mulcal(int p,ll v)
{
    ll c=v;
    sum[p]=(sum[p]*c)%mod;
    c*=v;c%=mod;
    sum2[p]=(sum2[p]*c)%mod;
    c*=v;c%=mod;
    sum3[p]=(sum3[p]*c)%mod;
}
void PushDown(int l,int r,int rt)
{
    int ls=rt<<1,rs=rt<<1|1;
    int m=(l+r)>>1;
    if(mulv[rt]!=1)
    {
        mulv[ls]=(mulv[ls]*mulv[rt])%mod;
        mulv[rs]=(mulv[rs]*mulv[rt])%mod;
        addv[ls]=(addv[ls]*mulv[rt])%mod;
        addv[rs]=(addv[rs]*mulv[rt])%mod;
        mulcal(ls,mulv[rt]);
        mulcal(rs,mulv[rt]);
        mulv[rt]=1;
    }
    if(addv[rt]!=0)
    {
        addv[ls]=(addv[rt]+addv[ls])%mod;
        addv[rs]=(addv[rt]+addv[rs])%mod;
        sumcal(ls,addv[rt],m-l+1);
        sumcal(rs,addv[rt],r-m);
        addv[rt]=0;
    }
}
void build(int l,int r,int rt)
{
    sum[rt]=sum2[rt]=sum3[rt]=0;
    addv[rt]=0;
    mulv[rt]=1;
    if(l==r) return ;
    int m=(l+r)>>1;
    build(l,m,rt<<1);
    build(m+1,r,rt<<1|1);
}
ll Query(int L,int R,int l,int r,int rt,int p)
{
    if(l>=L&&r<=R)
    {
        if(p==1) return sum[rt];
        else if(p==2) return sum2[rt];
        else return sum3[rt];
    }
    PushDown(l,r,rt);
    int m=(l+r)>>1;
    ll sum=0;
    if(m>=L)
      sum+=Query(L,R,l,m,rt<<1,p);
    if(m<R)
      sum+=Query(L,R,m+1,r,rt<<1|1,p);
    return sum%mod;
}
void Update(int L,int R,int l,int r,int rt,ll c,int op)
{
    if(l>=L&&r<=R)
    {
        if(op==3)
        {
            addv[rt]=c;mulv[rt]=0;
            sum[rt]=((r-l+1)*c)%mod;
            sum2[rt]=(sum[rt]*c)%mod;
            sum3[rt]=(sum2[rt]*c)%mod;
            return ;
        }
        else if(op==1)
        {
            addv[rt]+=c;addv[rt]%=mod;
            sumcal(rt,c,r-l+1);
            return ;
        }
        else
        {
            mulv[rt]*=c;mulv[rt]%=mod;
            addv[rt]*=c;addv[rt]%=mod;
            mulcal(rt,c);
            return ;
        }
    }
    PushDown(l,r,rt);
    int m=(l+r)>>1;
    if(m>=L)
      Update(L,R,l,m,rt<<1,c,op);
    if(m<R)
      Update(L,R,m+1,r,rt<<1|1,c,op);
    PushUp(rt);
}
int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    int n,m;
    while(~scanf("%d%d",&n,&m))
    {
        if(n==0&&m==0) break;
        build(1,n,1);
        int op,x,y;
        ll c;
        for(int i=0;i<m;++i)
        {
            scanf("%d%d%d%I64d",&op,&x,&y,&c);
            if(op==4)
            {
                ll ans=Query(x,y,1,n,1,c);
                printf("%I64d\n",ans);
            }
            else
                Update(x,y,1,n,1,c,op);
        }
    }
    return 0;
}