HDU - 4578 Transformation(区间的加，乘，更改三种操作，线段树如何维护)

题目链接：https://vjudge.net/contest/346897#problem/E

翻译：
给定n和m，分别表示1~n的区间(刚开始每个数都是0)，接下来m种操作
op x y z
op=1时，区间[x.y]内的每个数加z
op=2时，区间[x,y]内的每个数乘z
op=3时，区间[x,y]内的每个数更改为z
op=4时，区间[x,y]内的每个数的z次幂求和
分析：
sum1，sum2，sum3分别表示对应区间每个数的1，2，3次幂的和
对于op=1，2，3，这三种操作的优先级顺序
更改值>乘法>加法
对应线段树的维护需要三个标记
更改值标记更新时，对应其余两个标记也需要更新
乘法标记更新时，对应加法标记也需要更新
对于三次方，二次方代入数学公式即可.
代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
const int mod=10007;
const int maxx=100010;
int n,m,sum;
void pushup(int rt)
{
    sum1[rt]=(sum1[rt<<1]+sum1[rt<<1|1])%mod;
    sum2[rt]=(sum2[rt<<1]+sum2[rt<<1|1])%mod;
    sum3[rt]=(sum3[rt<<1]+sum3[rt<<1|1])%mod;
}
void build(int l,int r,int rt)
{
    add[rt]=same[rt]=0;
    mult[rt]=1;
    if(l==r)
    {
        sum1[rt]=sum2[rt]=sum3[rt]=0;
        return;
    }
    int mid=(l+r)>>1;
    build(l,mid,rt<<1);
    build(mid+1,r,rt<<1|1);
    pushup(rt);
}
void pushdown(int rt,int len)    //下传标记
{
    if(same[rt])
    {
        sum1[rt<<1]=(len-(len>>1))*same[rt]%mod;
        sum1[rt<<1|1]=(len>>1)*same[rt]%mod;
        sum2[rt<<1]=((len-(len>>1))*same[rt]%mod)*same[rt]%mod;
        sum2[rt<<1|1]=((len>>1)*same[rt]%mod)*same[rt]%mod;
        sum3[rt<<1]=((len-(len>>1))*same[rt]%mod)*(same[rt]*same[rt]%mod)%mod;
        sum3[rt<<1|1]=((len>>1)*same[rt]%mod)*(same[rt]*same[rt]%mod)%mod;
        same[rt<<1]=same[rt<<1|1]=same[rt];
        add[rt<<1]=add[rt<<1|1]=0;  /*加操作和乘操作的标记都需要改变*/
        mult[rt<<1]=mult[rt<<1|1]=1; 
        same[rt]=0;
    }
    if(mult[rt]!=1)
    {
        mult[rt<<1]=mult[rt<<1]*mult[rt]%mod;
        mult[rt<<1|1]=mult[rt<<1|1]*mult[rt]%mod;
        
        add[rt<<1]=add[rt<<1]*mult[rt]%mod; //改变
        add[rt<<1|1]=add[rt<<1|1]*mult[rt]%mod;  //改变
        
        sum1[rt<<1]=sum1[rt<<1]*mult[rt]%mod;
        sum1[rt<<1|1]=sum1[rt<<1|1]*mult[rt]%mod;
        sum2[rt<<1]=(sum2[rt<<1]*mult[rt]%mod)*mult[rt]%mod;
        sum2[rt<<1|1]=(sum2[rt<<1|1]*mult[rt]%mod)*mult[rt]%mod;
        sum3[rt<<1]=(sum3[rt<<1]*mult[rt]%mod)*(mult[rt]*mult[rt]%mod)%mod;
        sum3[rt<<1|1]=(sum3[rt<<1|1]*mult[rt]%mod)*(mult[rt]*mult[rt]%mod)%mod;
        mult[rt]=1;
    }
    if(add[rt])
    {
        add[rt<<1]=(add[rt<<1]+add[rt])%mod;
        add[rt<<1|1]=(add[rt<<1|1]+add[rt])%mod;
        
        sum3[rt<<1]=(sum3[rt<<1]+3*sum2[rt<<1]*add[rt]%mod+3*(sum1[rt<<1]*add[rt]%mod)*add[rt]%mod+(add[rt]*add[rt]%mod)*(add[rt]*(len-(len>>1))%mod))%mod;
        sum3[rt<<1|1]=(sum3[rt<<1|1]+3*sum2[rt<<1|1]*add[rt]%mod+3*(sum1[rt<<1|1]*add[rt]%mod)*add[rt]%mod+(add[rt]*add[rt]%mod)*(add[rt]*(len>>1)%mod))%mod;
        sum2[rt<<1]=(sum2[rt<<1]+2*sum1[rt<<1]*add[rt]%mod+(add[rt]*add[rt]%mod)*(len-(len>>1))%mod)%mod;
        sum2[rt<<1|1]=(sum2[rt<<1|1]+2*sum1[rt<<1|1]*add[rt]%mod+(add[rt]*add[rt]%mod)*(len>>1)%mod)%mod;
        sum1[rt<<1]=(sum1[rt<<1]+(len-(len>>1))*add[rt]%mod)%mod;
        sum1[rt<<1|1]=(sum1[rt<<1|1]+(len>>1)*add[rt]%mod)%mod;
        add[rt]=0;
    }
}
void updata(int L,int R, int val,int op,int l,int r,int rt)
{
    if(L<=l&&R>=r)
    {
        if(op==3)
        {
            same[rt]=val;
            mult[rt]=1;  //改变
            add[rt]=0;   //改变
            sum1[rt]=(r-l+1)*val%mod;
            sum2[rt]=((r-l+1)*val%mod)*val%mod;
            sum3[rt]=((r-l+1)*val%mod)*(val*val%mod)%mod;
        }
        else if(op==2)
        {
            mult[rt]=mult[rt]*val%mod;
            add[rt]=add[rt]*val%mod; //改变
            sum1[rt]=sum1[rt]*val%mod;
            sum2[rt]=(sum2[rt]*val%mod)*val%mod;
            sum3[rt]=(sum3[rt]*val%mod)*(val*val%mod)%mod;
        }
        else if(op==1)
        {
            add[rt]=(add[rt]+val)%mod;
            sum3[rt]=(sum3[rt]+3*sum2[rt]*val%mod+3*(sum1[rt]*val%mod)*val%mod+(((r-l+1)*val%mod)*(val*val%mod)%mod))%mod;
            sum2[rt]=(sum2[rt]+2*sum1[rt]*val%mod+((r-l+1)*val%mod)*val%mod)%mod;
            sum1[rt]=(sum1[rt]+(r-l+1)*val%mod)%mod;
        }
        return;
    }
    pushdown(rt,r-l+1);
    int mid=(l+r)>>1;
    if(L<=mid)
        updata(L,R,val,op,l,mid,rt<<1);
    if(R>mid)
        updata(L,R,val,op,mid+1,r,rt<<1|1);
    pushup(rt);
}
ll query(int L,int R,int val,int l,int r,int rt)
{
    if(L<=l&&R>=r)
    {
        if(val==1)
        {
            return sum1[rt]%mod;
        }
        else if(val==2)
        {
            return sum2[rt]%mod;
        }
        else
            return sum3[rt]%mod;
    }
    ll sum=0;
    pushdown(rt,r-l+1);
    int mid=(l+r)>>1;
    if(L<=mid)
        sum=(sum+query(L,R,val,l,mid,rt<<1))%mod;
    if(R>mid)
        sum=(sum+query(L,R,val,mid+1,r,rt<<1|1))%mod;
    return sum;
}
int main()
{
    while(~scanf("%d%d",&n,&m))
    {
        if(n==0&&m==0)
            break;
        build(1,n,1);
        while(m--)
        {
            int op,x,y,val;
            scanf("%d%d%d%d",&op,&x,&y,&val);
            if(op==4)
            {
                printf("%lld\n",query(x,y,val,1,n,1));
            }
            else
                updata(x,y,val,op,1,n,1);
        }
    }
}

业，精于勤，荒于嬉
行，成于思，毁于随