HDU 4578 Transformation 解题报告（线段树，2013杭州邀请赛）

    题目大意：对一组数据有4种操作：对区间[l,r]的加法，乘法，覆盖，以及查询1次方，2次方，3次方和。

    解题报告：当时杭州邀请赛中的一道线段树裸题。记得是当时我们都卡在别的题目上，都挺低落的。最后做这题时时间也不够了，导致最后打铁了。当然，还是实力问题。

    线段树裸题，要注意的是两点：一是加法更新时先更新3次方和，再更新2次方和，最后更新1次方和。自己列个公式很容易推出来。

    二是标记传递时的处理顺序，先是覆盖标记，再是乘法标记，最后是加法标记。

    然后注意细节，随时取模，很容易就A了。

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;

#define LL long long
#define defm int m=(l+r)>>1;
#define lson l,m,pos<<1
#define rson m+1,r,pos<<1|1

const int MOD=10007;
const int maxn=111111;
int add[maxn<<2],mul[maxn<<2],cov[maxn<<2];
int sum[maxn<<2][3];

void updateFather(int pos)
{
	for(int i=0;i<3;i++)
		sum[pos][i]=(sum[pos<<1][i]+sum[pos<<1|1][i])%MOD;
}

void Add(int v,int amount,int pos)
{
	sum[pos][2]=sum[pos][2]+3*((v*sum[pos][1])%MOD)+3*(((v*v)%MOD*sum[pos][0])%MOD)+(((v*v)%MOD)*((v*amount)%MOD))%MOD;
	sum[pos][2]%=MOD;

	sum[pos][1]=sum[pos][1]+2*(v*sum[pos][0]%MOD)+(((v*v)%MOD)*amount)%MOD;
	sum[pos][1]%=MOD;

	sum[pos][0]=sum[pos][0]+(v*amount)%MOD;
	sum[pos][0]%=MOD;
	add[pos]=(add[pos]+v)%MOD;
}

void Mul(int v,int amount,int pos)
{
	sum[pos][2]=(((v*v)%MOD)*((v*sum[pos][2])%MOD))%MOD;
	sum[pos][1]=(((v*v)%MOD)*sum[pos][1])%MOD;;
	sum[pos][0]=(v*sum[pos][0])%MOD;

	mul[pos]=(mul[pos]*v)%MOD;
	add[pos]=(add[pos]*v)%MOD;
}

void Cov(int v,int amount,int pos)
{
	sum[pos][0]=(v*amount)%MOD;
	sum[pos][1]=(((v*v)%MOD)*amount)%MOD;
	sum[pos][2]=(((v*v)%MOD)*((v*amount)%MOD))%MOD;
	cov[pos]=v;
	mul[pos]=1;
	add[pos]=0;
}

void updateSon(int amount,int pos)
{
	if(~cov[pos])
	{
		Cov(cov[pos],amount-amount/2,pos<<1);
		Cov(cov[pos],amount/2,pos<<1|1);
		cov[pos]=-1;
	}
	if(mul[pos]!=1)
	{
		Mul(mul[pos],amount-amount/2,pos<<1);
		Mul(mul[pos],amount/2,pos<<1|1);
		mul[pos]=1;
	}
	if(add[pos])
	{
		Add(add[pos],amount-amount/2,pos<<1);
		Add(add[pos],amount/2,pos<<1|1);
		add[pos]=0;
	}
}

void update(int L,int R,int op,int v,int l,int r,int pos)
{
	if(L<=l && r<=R)
	{
		switch(op)
		{
		case 1: Add(v,r-l+1,pos);break;
		case 2: Mul(v,r-l+1,pos);break;
		case 3: Cov(v,r-l+1,pos);break;
		}
		return ;
	}
	updateSon(r-l+1,pos);
	defm;
	if(L<=m)
		update(L,R,op,v,lson);
	if(m<R)
		update(L,R,op,v,rson);
	updateFather(pos);
}

int query(int L,int R,int p,int l,int r,int pos)
{
	if(L<=l && r<=R)
		return sum[pos][p];
	updateSon(r-l+1,pos);
	defm;
	return ((L<=m?query(L,R,p,lson):0)+(m<R?query(L,R,p,rson):0))%MOD;
}

int main()
{
	int n,m;
	while(~scanf("%d%d",&n,&m) && (n|m))
	{
		update(1,n,3,0,1,n,1);

		while(m--)
		{
			int a,b,c,d;
			scanf("%d%d%d%d",&a,&b,&c,&d);
			if(a==4)
				printf("%d\n",query(b,c,d-1,1,n,1));
			else
				update(b,c,a,d,1,n,1);
		}
	}
}