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[3],add,mul;
    void init(int _l,int _r){
        l=_l;r=_r;mul=1;add=0;
    }
    void Add(int a){
        sum[2]=(sum[2]+a*a%M*a%M*len()%M+3*a*sum[1]%M+3*a*a%M*sum[0]%M)%M;
        sum[1]=(sum[1]+a*a%M*len()%M+2*a*sum[0]%M)%M;
        sum[0]=(sum[0]+a*len()%M)%M;
        add=(add+a)%M;
    }
    void Mul(int a){
        sum[2]=sum[2]*a%M*a%M*a%M;
        sum[1]=sum[1]*a%M*a%M;
        sum[0]=sum[0]*a%M;
        add=add*a%M;
        mul=mul*a%M;
    }
    int len(){
        return r-l+1;
    }
    int mid(){
        return (l+r)>>1;
    }
}T[MAX*4];
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=0;
    T[i].mul=1;
}
void updateson(int i)
{
    for(int j=0;j<3;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,0,sizeof(T[i].sum));
    if(l==r) return ;
    int m=T[i].mid();
    built(lson,l,m);
    built(rson,m+1,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+1,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-1];
    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+1,r,p))%M;
}
int main(){
    int m,n,op,x,y,p;
    while(scanf("%d%d",&n,&m)&&(m+n)){
        built(1,1,n);
        while(m--){
            scanf("%d%d%d%d",&op,&x,&y,&p);
            if(op==1) update(1,x,y,1,p);
            else if(op==2) update(1,x,y,p,0);
            else if(op==3) update(1,x,y,0,p);
            else printf("%d
",query(1,x,y,p)%M);
        }
    }
    return 0;
}