题意:
给出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;
}