HDU 4578 Transformation（2013 杭州赛区全国邀请赛 1003）

http://acm.hdu.edu.cn/showproblem.php?pid=4578

Transformation

Time Limit: 15000/8000 MS (Java/Others)    Memory Limit: 65535/65536 K (Java/Others)
Total Submission(s): 1483    Accepted Submission(s): 341

Problem Description

Yuanfang is puzzled with the question below: 
There are n integers, a ₁, a ₂, …, a _n. The initial values of them are 0. There are four kinds of operations.
Operation 1: Add c to each number between a _x and a _y inclusive. In other words, do transformation a _k<---a _k+c, k = x,x+1,…,y.
Operation 2: Multiply c to each number between a _x and a _y inclusive. In other words, do transformation a _k<---a _k×c, k = x,x+1,…,y.
Operation 3: Change the numbers between a _x and a _y to c, inclusive. In other words, do transformation a _k<---c, k = x,x+1,…,y.
Operation 4: Get the sum of p power among the numbers between a _x and a _y inclusive. In other words, get the result of a _x ^p+a _x+1 ^p+…+a _y  ^p.
Yuanfang has no idea of how to do it. So he wants to ask you to help him. 

 

Input

There are no more than 10 test cases.
For each case, the first line contains two numbers n and m, meaning that there are n integers and m operations. 1 <= n, m <= 100,000.
Each the following m lines contains an operation. Operation 1 to 3 is in this format: "1 x y c" or "2 x y c" or "3 x y c". Operation 4 is in this format: "4 x y p". (1 <= x <= y <= n, 1 <= c <= 10,000, 1 <= p <= 3)
The input ends with 0 0.

 

Output

For each operation 4, output a single integer in one line representing the result. The answer may be quite large. You just need to calculate the remainder of the answer when divided by 10007.

思路：看题意可以很明显地看出是一道线段树题，也是一道很典型的线段树，但是在现场赛上我还真不能确定自己能不能把它搞定，题意相当简单，但是处理起来还是挺麻烦的，写了半天，WA了两发才过的。。。

题意相当明显，就不多说了，有三种操作，这三种操作顺序不一样的话得到的结果是不一样的，对于这种多操作的问题，我们要做的事尽量将多种操作合并成一种操作。

我们可将区间中的数看成是 ax+b的形式，对于加c操作，则变成 ax+b+c（b->b+c）,对于乘c操作，则变成 acx+bc,（a->ac，b->bc）对于赋值c操作，则变成c，即（a->1，x->c,b->0）,则我们可以在线段数中加以下标记，a,b,x分别是以上提到的变量，sum[3]，表示答案，sum[0],sum[1],sum[2]分别表示1次方和，平方和，立方和。对于更新和查询操作我们用pushdown函数向下更新。对于维护平方和还有立方和的问题，我们只要将平方，立方展开再利用更新之前的值即可维护，具体方法这里不多说了。

代码写的有点乱，仅供参考。。。

#include <iostream>
#include <string.h>
#include <stdio.h>
#include <algorithm>
#define maxn 100010
#define mod 10007
#define mid ((t[p].l+t[p].r)>>1)
#define ls (p<<1)
#define rs (ls|1)
using namespace std;
struct tree
{
    int l,r;
    int a,b,c;
    int sum[3];
}t[maxn<<2];
void pushup(int p)
{
    int i;
    for(i=0;i<3;i++)
    t[p].sum[i]=(t[ls].sum[i]+t[rs].sum[i])%mod;
}
void pushdown(int p)
{
    if(t[p].l==t[p].r)
    return;
     int a=t[p].a,b=t[p].b,c=t[p].c;
     if(c)
     {
         t[ls].a=t[rs].a=a;
         t[ls].b=t[rs].b=b;
         t[ls].c=t[rs].c=c;
         t[ls].sum[2]=(t[ls].r-t[ls].l+1)*(a*c%mod+b)%mod*(a*c%mod+b)%mod*(a*c%mod+b)%mod;
         t[ls].sum[1]=(t[ls].r-t[ls].l+1)*(a*c%mod+b)%mod*(a*c%mod+b)%mod;
         t[ls].sum[0]=(t[ls].r-t[ls].l+1)*(a*c%mod+b)%mod;
         t[rs].sum[2]=(t[rs].r-t[rs].l+1)*(a*c%mod+b)%mod*(a*c%mod+b)%mod*(a*c%mod+b)%mod;
         t[rs].sum[1]=(t[rs].r-t[rs].l+1)*(a*c%mod+b)%mod*(a*c%mod+b)%mod;
         t[rs].sum[0]=(t[rs].r-t[rs].l+1)*(a*c%mod+b)%mod;
     }
     else
     {
         t[ls].a=(t[ls].a*a)%mod;
         t[ls].b=(t[ls].b*a+b)%mod;
         t[rs].a=(t[rs].a*a)%mod;
         t[rs].b=(t[rs].b*a+b)%mod;
         t[ls].sum[2]=(a*a%mod*a%mod*t[ls].sum[2]%mod+3*a%mod*a%mod*b%mod*t[ls].sum[1]%mod+3*a%mod*b%mod*b%mod*t[ls].sum[0]%mod+b*b%mod*b%mod*(t[ls].r-t[ls].l+1)%mod)%mod;
         t[ls].sum[1]=(a*a%mod*t[ls].sum[1]%mod+b*b%mod*(t[ls].r-t[ls].l+1)%mod+2*a*b%mod*t[ls].sum[0]%mod)%mod;
         t[ls].sum[0]=(a*t[ls].sum[0]%mod+b*(t[ls].r-t[ls].l+1)%mod)%mod;
         t[rs].sum[2]=(a*a%mod*a%mod*t[rs].sum[2]%mod+3*a%mod*a%mod*b%mod*t[rs].sum[1]%mod+3*a%mod*b%mod*b%mod*t[rs].sum[0]%mod+b*b%mod*b%mod*(t[rs].r-t[rs].l+1)%mod)%mod;
         t[rs].sum[1]=(a*a%mod*t[rs].sum[1]%mod+b*b%mod*(t[rs].r-t[rs].l+1)%mod+2*a*b%mod*t[rs].sum[0]%mod)%mod;
         t[rs].sum[0]=(a*t[rs].sum[0]%mod+b*(t[rs].r-t[rs].l+1)%mod)%mod;
     }
     t[p].b=t[p].c=0;
     t[p].a=1;
}
void build(int p,int l,int r)
{
    t[p].l=l,t[p].r=r;
    t[p].a=1,t[p].b=t[p].c=0;
    t[p].sum[0]=t[p].sum[1]=t[p].sum[2]=0;
    if(l<r)
    {
        build(ls,l,mid);
        build(rs,mid+1,r);
    }
}
void change(int p,int l,int r,int val,int flag)
{
    if(t[p].l==l&&t[p].r==r)
    {
        if(flag==0)
        {
            t[p].b=(t[p].b+val)%mod;
            t[p].sum[2]=(t[p].sum[2]+3*val%mod*t[p].sum[1]%mod+3*val%mod*val%mod*t[p].sum[0]%mod+val*val%mod*val%mod*(t[p].r-t[p].l+1)%mod)%mod;
            t[p].sum[1]=(t[p].sum[1]+val*val%mod*(t[p].r-t[p].l+1)%mod+2*val*t[p].sum[0]%mod)%mod;
            t[p].sum[0]=(t[p].sum[0]+val*(t[p].r-t[p].l+1))%mod;
        }
        else if(flag==1)
        {
            t[p].a=(t[p].a*val)%mod;
            t[p].b=(t[p].b*val)%mod;
            t[p].sum[2]=val*val%mod*val%mod*t[p].sum[2]%mod;
            t[p].sum[1]=val*val%mod*t[p].sum[1]%mod;
            t[p].sum[0]=val*t[p].sum[0]%mod;
        }
        else if(flag==2)
        {
            t[p].a=1;
            t[p].b=0;
            t[p].c=val;
            t[p].sum[2]=(t[p].r-t[p].l+1)%mod*val%mod*val%mod*val%mod;
            t[p].sum[1]=(t[p].r-t[p].l+1)%mod*val%mod*val%mod;
            t[p].sum[0]=(t[p].r-t[p].l+1)*val%mod;
        }
        return;
    }
    pushdown(p);
    if(l>mid)
    change(rs,l,r,val,flag);
    else if(r<=mid)
    change(ls,l,r,val,flag);
    else
    {
        change(ls,l,mid,val,flag);
        change(rs,mid+1,r,val,flag);
    }
    pushup(p);
}
int query(int p,int l,int r,int flag)
{
    if(t[p].l==l&&t[p].r==r)
    return t[p].sum[flag];
    pushdown(p);
    if(l>mid)
    return query(rs,l,r,flag);
    else if(r<=mid)
    return query(ls,l,r,flag);
    else
    {
        return (query(ls,l,mid,flag)+query(rs,mid+1,r,flag))%mod;
    }
}
int main()
{
    freopen("dd.txt","r",stdin);
    int n,m;
    while(scanf("%d%d",&n,&m)&&(n+m))
    {
        build(1,1,n);
        int q,x,y,c;
        while(m--)
        {
            scanf("%d%d%d%d",&q,&x,&y,&c);
            if(q==4)
            {
                printf("%d\n",query(1,x,y,c-1));
            }
            else
            {
                change(1,x,y,c,q-1);
            }
        }
    }
    return 0;
}