hdu 5239 Doom（线段树+数论）

题意：

给出n个数和一个初始值为0的答案。每次操作给出一个区间[l,r]，把区间所有的数加到答案中，之后把区间的每个数都平方。每次操作都需要输出答案 mod 9223372034707292160(2 ^ 63 - 2 ^ 31)

思路

这题的做法和 hdu 4027 类似
首先，应该想到这个题可以用线段树来维护区间的和，然后对于(2 ^ 63 - 2 ^ 31)这个数，即任意数的平方MOD此数，重复操作至多29次就会进入一个不变的数。
因此，可以在线段树维护时，用一个标记数组cover来维护。如果发现当前的是叶区间且当前数字的平方去模等于当前数字cover[o] = true。
向上维护cover[o] = cover[ls] & cover[rs]。表示向上维护的区间是否还需要更新。

AC代码：

#include <cstdio>
#include <cstring>
#include <algorithm>
#define ls (o*2)
#define rs (o*2+1)
using namespace std;
typedef unsigned __int64 ll;
const int N = 100005;
const ll MOD = 9223372034707292160ll;
ll sumv[N<<2];
bool cover[N<<2];
int n, q;

ll modmul(ll a, ll k) {
    ll c = 0;
    while(k) {
        if(k & 1) c = (c+a) % MOD;
        a = (a+a) % MOD;
        k >>= 1;
    }
    return c;
}

void build(int o, int L, int R) {
    sumv[o] = cover[o] = 0;
    if(L == R) {
        scanf("%I64u", &sumv[o]);
        return ;
    }
    int M = (L+R)/2;
    build(ls, L, M);
    build(rs, M+1, R);
    sumv[o] = (sumv[ls] + sumv[rs]) % MOD;
}

int ql, qr;
ll query(int o, int L, int R) {
    if(ql <= L && R <= qr) {
        return sumv[o];
    }
    int M = (L+R)/2;
    ll ret = 0;
    if(ql <= M) ret = (ret + query(ls, L, M)) % MOD; 
    if(qr > M) ret = (ret + query(rs, M+1, R)) % MOD;
    return ret;
}

void modify(int o, int L, int R) {
    if(cover[o] && ql <= L && R <= qr)
        return ;
    if(L == R) {
        ll tmp = modmul(sumv[o], sumv[o]);
        if(tmp == sumv[o])
            cover[o] = true;
        sumv[o] = tmp;
        return ;
    }
    int M = (L+R)/2;
    if(ql <= M) modify(ls, L, M);
    if(qr > M) modify(rs, M+1, R);
    cover[o] = cover[ls] & cover[rs];
    sumv[o] = (sumv[ls] + sumv[rs]) % MOD;
}

int main() {
    int T, cas = 1;
    scanf("%d", &T);
    while(T--) {
        scanf("%d%d", &n, &q);
        build(1, 1, n);
        printf("Case #%d:\n", cas++);
        ll ans = 0;
        while(q--) {
            scanf("%d%d", &ql, &qr);
            ans = (ans + query(1, 1, n)) % MOD;
            printf("%I64u\n", ans);
            modify(1, 1, n);
        }
    }
    return 0;
}