HDU - 4578 Transformation（线段树）

异常变态的正常线段树。给你三种操作，三种输出方式：

1.让l~r区间的所有值加c；

2.让l~r区间的所有值乘c；

3.让l~r区间的所有值变成c；

4.让l~r区间的所有值都先求p（1~3）次方，然后求和输出。

别人都是一看题懵逼好简单操作好复杂，我就比较优秀，一看题就懵逼感觉很难不知道怎么做，我对不起这几天努力学的线段树，对不起。其实我考虑到分三种lazy数组，三种ans数组（正常大家都能想到。。。），但是我不知道数组到开2次方，3次方怎么做。再此跟初中老师道歉，大家看下面。

sum1 = a1 + a2 + a3;

sum2 = a1^2 + a2^2 + a3^2;

sum3 = a1^3 + a2^3 + a3^3;

如果每个元素都加c呢？

sum1 = a1 + a2 + a3 + 3*c;

sum2 = (a1 + c)^2 + (a2 + c) ^2 + (a3 + c)^2;

sum3 = (a1 + c)^3 + (a2 + c) ^3 + (a3 + c)^3;

真的傻，我竟然不知道怎么处理sum2，式子分解一下不就好了

sum2 = a1^2 + a2^2 + a3^2 + 2 * c * (a + b + c) + 3 * c^2;

sum2 = 上一状态的sum2 + 2 * c * 上一状态的sum1 + 3 * c^2;

sum3 = .....

每个元素乘c，这个就很简单了，大家都懂，每个元素变成c也很好操作。

这只是其中一种通俗易懂代码写起来不要命的操作，让我写代码不可能的，代码转自：https://www.cnblogs.com/H-Vking/p/4297973.html

#include <iostream>
#include <cstdio>
#include <vector>
#include <cmath>
#include <string>
#include <string.h>
#include <algorithm>
using namespace std;
#define LL __int64
#define eps 1e-8
#define INF INT_MAX
#define lson l , m , rt << 1
#define rson m + 1 , r , rt << 1 | 1
const int MOD = 10007; 
const int maxn = 100000 + 5;
const int N = 12;
LL add[maxn << 2] , set[maxn << 2] , mul[maxn << 2];
LL sum1[maxn << 2] , sum2[maxn << 2] , sum3[maxn << 2];
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] = set[rt] = 0;
    mul[rt] = 1;
    if(l == r) {
        sum1[rt] = sum2[rt] = sum3[rt] = 0;
        return;
    }
    int m = (l + r) >> 1;
    build(lson);
    build(rson);
    PushUp(rt);
}
void PushDown(int rt , int len)
{
    if(set[rt]) {
        set[rt << 1] = set[rt << 1 |