HDU 5152 线段树+欧拉函数

题意：
本题是求对一个长为$n$序列进行$m$此操作。$1$、区间加，$2$、节点$l—r$，$x$值变为$2^x$,3、求$l--r$的和。

解题思路：
对于操作二，利用公式 当$x >= Phi(C)$, $A^x$ =$A$^ $(x%Phi(C)$ + $Phi(C)) (mod C)$ 对于$2333333$这个模数来说，求$18$次欧拉函数后就变成了$1$，所以只需要保存$19$层第三次操作的加数即可，然后就直接是线段树区间更新和询问操作了。 Ps:很多人把题目想简单了，操作二不能直接简单的更新，不满足$2$^$2^x$ mod P= $2$^$( 2^x mod P) mod P$。
(确实不知道这个定理，重新认识了欧拉函数T_T，保留了官方给出的题解以后留着看吧。)

#include <cstdio>
#include <cstring>
#include <vector>
#include <algorithm>
#include <ctime>
#include <iostream>
using namespace std;

typedef __int64 ll;
const int MOD = 2333333;
const int N = 50000+5;

int mo[20] = {2333333,2196720,580608,165888,55296,18432,
6144,2048,1024,512,256,128,64,32,16,8,4,2,1};
int pow2[33];

vector<ll> vt[N];
struct Segtree {
#define lson rt<<1, l, mid
#define rson rt<<1|1, mid+1, r
    ll mark[N<<2];
    int length[N<<2], sum[N<<2];
    void build(int rt, int l, int r) {
        mark[rt] = 0;
        length[rt] = r-l+1;
        if(l == r) {
            vt[l].clear();
            int x;
            scanf("%d", &x);
            vt[l].push_back(x);
            sum[rt] = x%MOD;
            return ;
        }
        int mid = l+r>>1;
        build(lson); build(rson);
        up(rt);
    }

    inline void Add(int &x, int y) {
        x += y;
        if(x >= MOD)    x -= MOD;
    }

    void up(int rt) {
        sum[rt] = sum[rt<<1] + sum[rt<<1|1];
        Add(sum[rt], 0);
    }

    void down(int rt) {
        if(mark[rt]) {
            mark[rt<<1] += mark[rt];
            mark[rt<<1|1] += mark[rt];
            Add(sum[rt<<1], 1LL*length[rt<<1]*mark[rt]%MOD);
            Add(sum[rt<<1|1], 1LL*length[rt<<1|1]*mark[rt]%MOD);
            mark[rt] = 0;
        }
    }

    void update(int rt, int l, int r, int L, int R, int v) {
        if(L <= l && R >= r) {
            mark[rt] += v;
            Add(sum[rt], 1LL*length[rt]*v%MOD);
            return ;
        }
        down(rt);
        int mid = l+r>>1;
        if(L <= mid)    update(lson, L, R, v);
        if(R > mid) update(rson, L, R, v);
        up(rt);
    }

    int pow_mod(int x, int n, int mod) {
        int ret = 1%mod;
        while(n) {
            if(n&1) ret = 1LL*ret*x%mod;
            x = 1LL*x*x%mod;
            n >>= 1;
        }
        return ret;
    }

    int cal(vector<ll> &v) {
        if(v.size() < 19) {
            int pos = v.size()-1;
            ll ret = v[0];
            bool flag = false;
            if(v[0] >= mo[pos]) {
                flag = true;
                ret = ret%mo[pos] + mo[pos];
            }
            pos--;
            for(int i = 1;i < v.size(); i++) {
                if(flag) {
                    ret = (pow_mod(2, ret, mo[pos]) +v[i])%mo[pos] + mo[pos];
                }
                else {
                    if(ret >= 30) {
                        flag = true;
                        ret = (pow_mod(2, ret, mo[pos]) + v[i])%mo[pos]+mo[pos];
                    }
                    else if(pow2[ret] >= mo[pos]) {
                        flag = true;
                        ret = (pow_mod(2, ret, mo[pos]) + v[i])%mo[pos] + mo[pos];
                    }
                    else {
                        ret = pow2[ret] + v[i];
                        if(ret >= mo[pos]) {
                            flag = true;
                            ret = ret%mo[pos] + mo[pos];
                        }
                    }
                }
                pos--;
            }
            return ret%MOD;
        }
        else {
            ll ret = 1;
            int pos = 17;
            bool flag = true;
            for(int i = v.size()-18;i < v.size(); i++) {
                ret = (pow_mod(2, ret, mo[pos]) + v[i])%mo[pos] + mo[pos];
                pos--;
            }
            return ret%MOD;
        }
    }

    void modify(int rt, int l, int r, int x) {
        if(l == r) {
            if(mark[rt]) {
                vt[l][vt[l].size()-1] += mark[rt];
                mark[rt] = 0;
            }
            vt[l].push_back(0);
            sum[rt] = cal(vt[l]);
            return ;
        }
        down(rt);
        int mid = l+r>>1;
        if(x <= mid)    modify(lson, x);
        else    modify(rson, x);
        up(rt);
    }

    int query(int rt, int l, int r, int L, int R) {
        if(L <= l && R >= r)    return sum[rt];
        down(rt);
        int mid = l+r>>1, ret = 0;
        if(L <= mid)    Add(ret, query(lson, L, R));
        if(R > mid) Add(ret, query(rson, L, R));
        up(rt);
        return ret;
    }
}tree;

int n, m;

void init() {
    for(int i = 0;i <= 30; i++) pow2[i] = 1<<i;
}

int main() {
    init();
    int op, l, r, x;
    while(scanf("%d%d", &n, &m) == 2) {
        tree.build(1, 1, n);
        while(m--) {
            scanf("%d", &op);
            if(op == 1) {
                scanf("%d%d", &l, &r);
                printf("%d\n", tree.query(1, 1, n, l, r));
            }
            else if(op == 2) {
                scanf("%d", &x);
                tree.modify(1, 1, n, x);
            }
            else {
                scanf("%d%d%d", &l, &r, &x);
                tree.update(1, 1, n, l, r, x);
            }
        }
    }
    return 0;
}