[HDOJ4027]Can you answer these queries?（线段树，特殊成段更新，成段查询）

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4027

RT，该题要求每次更新是更新所有节点，分别求平方根，查询是求和。昨晚思前想后找有没有一个数学上的开平方的和等于和的开平方之类的规律。但是想了想发现这不就是小的时候，如果你这么想那老师就会骂死你的那个- -!

所以显然这个题是无法按套路成段更新了，懒惰标记也是没有用了，我们应该去按照区间更新每一个节点。结果TLE了一发，这说明这题不是这么搞，一定还有规律的。注意到题目给数据规模是2^63以及题目红字所述：开根号全都取下整，除了考虑用longlong以外，我们其实还可以想一下对一个longlong的数据开平方，最终都会变成1。那么在更新这个节点至1的那次update里，更新结束后一定会从叶子rt向上更新父亲，如果rt的兄弟节点也是1，那么说明rt的父亲已经不再需要更新了（因为1开平方还是1），这时候rt的父亲存的结果是1+1=2。也就是(rt+1-rt+1)。推广到更高的节点，换成区间来表示就是(r-l+1)。我们遇到这种节点就不需要更新了。

特别需要注意的是题目中x和y的大小…会出现x>y的情况……以及，每个case最后都要有一个额外的\n…………（PE到死WA到死）

 1 #include <algorithm>
 2 #include <iostream>
 3 #include <iomanip>
 4 #include <cstring>
 5 #include <climits>
 6 #include <complex>
 7 #include <fstream>
 8 #include <cassert>
 9 #include <cstdio>
10 #include <bitset>
11 #include <vector>
12 #include <deque>
13 #include <queue>
14 #include <stack>
15 #include <ctime>
16 #include <set>
17 #include <map>
18 #include <cmath>
19
20 using namespace std;
21
22 #define fr first
23 #define sc second
24 #define pb(a) push_back(a)
25 #define Rint(a) scanf("%d", &a)
26 #define Rll(a) scanf("%I64d", &a)
27 #define Rs(a) scanf("%s", a)
28 #define FRead() freopen("in", "r", stdin)
29 #define FWrite() freopen("out", "w", stdout)
30 #define Rep(i, len) for(LL i = 0; i < (len); i++)
31 #define For(i, a, len) for(LL i = (a); i < (len); i++)
32 #define Cls(a) memset((a), 0, sizeof(a))
33 #define Full(a) memset((a), 0x7f7f, sizeof(a))
34
35 typedef long long LL;
36 #define lrt rt << 1
37 #define rrt rt << 1 | 1
38 const LL maxn = 100100;
39 LL sum[maxn<<2];
40 LL n, q;
41
42 void pushUP(LL rt) {
43     sum[rt] = sum[lrt] + sum[rrt];
44 }
45
46 void build(LL l, LL r, LL rt) {
47     if(l == r) {
48         Rll(sum[rt]);
49         return;
50     }
51     LL m = (l + r) >> 1;
52     build(l, m, lrt);
53     build(m+1, r, rrt);
54     pushUP(rt);
55 }
56
57 void update(LL L, LL R, LL l, LL r, LL rt) {
58     if(sum[rt] == r - l + 1) return;
59     if(l == r) {
60         sum[rt] = LL(double(sqrt(sum[rt])));
61         return;
62     }
63     LL m = (l + r) >> 1;
64     if(m >= L) update(L, R, l, m, lrt);
65     if(m < R) update(L,R, m+1, r, rrt);
66     pushUP(rt);
67 }
68
69 LL query(LL L, LL R, LL l, LL r, LL rt) {
70     if(l >= L && R >= r) return sum[rt];
71     LL m = (l + r) >> 1;
72     LL ret = 0;
73     if(m >= L) ret += query(L, R, l, m, lrt);
74     if(m < R) ret += query(L, R, m+1, r, rrt);
75     return ret;
76 }
77
78 int main() {
79     // FRead();
80     LL orz = 1;
81     LL a, b, c;
82     while(~Rll(n)) {
83         Cls(sum);
84         printf("Case #%I64d:\n", orz++);
85         build(1, n, 1);
86         Rll(q);
87         while(q--) {
88             Rll(a); Rll(b); Rll(c);
89             if(b > c) swap(b, c);
90             if(a == 0) update(b, c, 1, n, 1);
91             if(a == 1) cout << query(b, c, 1, n, 1) << endl;
92         }
93         printf("\n");
94     }
95     return 0;
96 }