HDU 5239 Doom（线段树）

题意：给出一个长度为n的序列，和m次查询，每次查询一个区间的和，但查询之后把区间内每个点的值平方一次，所有的结果模上9223372034707292160。

其实这题一点也不难。。。

一个数又平方又模的次数多了以后，会维持在一个固定的值不动。打表可以发现规律，但是仔细想想也是可以想明白的。如果x的素因子的集合包含模的素因子集合，那么经过足够多的次数之后，取模会直接为0；反之，据费马小定理，则会维持一个固定的值不变。

另外模没有超过ULL的一半，不需要CRT。。。（其实是写不对ULL的CRT。。似乎是扩展欧几里得那里有点问题

（用了输入输出挂，，另外把模改成减法了）

#include <algorithm>
#include <iostream>
#include <string.h>
#include <stdlib.h>
#include <assert.h>
#include <stdio.h>
#include <math.h>
using namespace std;
#define ULL unsigned long long

const ULL MOD = 9223372034707292160LL;
const int MAXN = 100010;
const int MAXL = 30;

template <class T>
inline int RD(T &x) {
x = 0;
char ch = getchar();
while(!isdigit(ch)) { ch = getchar(); if(ch == EOF) return 0; }
while(isdigit(ch)) { x *= 10; x += ch - '0'; ch = getchar(); }
return 1;
}

template <class T>
void PT(T x) {
if(x / 10) PT(x / 10);
putchar(x % 10 + '0');
}

ULL val[MAXN][MAXL];
ULL sum[MAXN << 2];
int cnt[MAXN << 2];

void Build(int l, int r, int rt) {
cnt[rt] = 0;
if(l >= r) {
sum[rt] = val[l][0];
return ;
}
int m = (l + r) >> 1;
Build(l, m, rt << 1);
Build(m + 1, r, rt << 1 | 1);
sum[rt] = sum[rt << 1] + sum[rt << 1 | 1];
if(sum[rt] > MOD) sum[rt] -= MOD;
}

void Update(int l, int r, int rt) {
cnt[rt]++;
if(l >= r) {
sum[rt] = val[l][cnt[rt]];
return ;
}
int m = (l + r) >> 1;
if(cnt[rt << 1] < MAXL - 1) Update(l, m, rt << 1);
if(cnt[rt << 1 | 1] < MAXL - 1) Update(m + 1, r, rt << 1 | 1);
sum[rt] = sum[rt << 1] + sum[rt << 1 | 1];
if(sum[rt] > MOD) sum[rt] -= MOD;
}

ULL Query(int L, int R, int l, int r, int rt) {
if(L <= l && r <= R) {
ULL ret = sum[rt];
if(cnt[rt] < MAXL - 1)
Update(l, r, rt);
return ret;
}
int m = (l + r) >> 1;
ULL lch = 0, rch = 0, ret;
if(L <= m) lch = Query(L, R, l, m, rt << 1);
if(R  > m) rch = Query(L, R, m + 1, r, rt << 1 | 1);
sum[rt] = sum[rt << 1] + sum[rt << 1 | 1];
if(sum[rt] > MOD) sum[rt] -= MOD;
ret = lch + rch;
if(ret > MOD) ret -= MOD;
return ret;
}

ULL Add_Mod(ULL x, ULL y) {
ULL ret = 0, p = x, k = y;
while(k) {
if(k % 2 == 1) {
ret += p;
if(ret > MOD) ret -= MOD;
}
p += p;
if(p > MOD) p -= MOD;
k /= 2;
}
return ret;
}

int main() {
//    freopen("D.in", "r", stdin);

int t, cas = 0;
ULL cd;
RD(t);
while(t--) {
putchar('C');
putchar('a');
putchar('s');
putchar('e');
putchar(' ');
putchar('#');
PT(++cas);
putchar(':');
puts("");
int n, m, a, b;
RD(n), RD(m);
for(int i = 1; i <= n; i++) {
RD(val[i][0]);
for(int j = 1; j < MAXL; j++)
val[i][j] = Add_Mod(val[i][j - 1], val[i][j - 1]);
}
Build(1, n, 1);
cd = 0;
while(m--) {
RD(a), RD(b);
cd = cd + Query(a, b, 1, n, 1);
if(cd > MOD) cd -= MOD;
PT(cd); puts("");
}
}
return 0;
}