傳送門

?

以為對于這類問題線段樹都能解決，分塊比線段樹菜，結果培訓完才知道線段樹是一種特殊的分塊方法，有的分塊的題線段樹不能做，看來分塊還是有必要學的。

對于這個題，先分塊，然后另開一個數組對于每個塊內排序。

區間加的話，加一個標記，每一個整塊區間加，里面的數的相對大小不變，而左右兩邊零散的塊直接暴力重構。

查詢可以對于每個塊二分查找。

時間復雜度應該是 nlogn + Q√nlog√n，剛好卡過。。

?

注意：第10個點會被卡，手寫二分比stl的lower_bound快一點，可以避免被卡。

　　　也可以在 S = sqrt(n) 后面 + x，x 不要太大也不要太小，不知道為什么，速度也會快點，玄學啊。

%%%有些dalao不知道怎么寫的，1000ms

?

——代碼（最終只能優化到1600ms）

1 #include <cmath> 2 #include <cstdio> 3 #include <iostream> 4 #include <algorithm> 5 #define LL long long 6 7 const int MAXN = 1000001; 8 int n, q, S, C; 9 int belong[MAXN], st[MAXN], ed[MAXN]; 10 LL a[MAXN], b[MAXN], add[MAXN]; 11 12 inline LL read() 13 { 14 LL x = 0, f = 1; 15 char ch = getchar(); 16 for(; !isdigit(ch); ch = getchar()) if(ch == '-') f = -1; 17 for(; isdigit(ch); ch = getchar()) x = x * 10 + ch - '0'; 18 return x * f; 19 } 20 21 inline int find(int x, int y, int z) 22 { 23 int mid; 24 while(x < y) 25 { 26 mid = (x + y) >> 1; 27 if(b[mid] >= z) y = mid; 28 else x = mid + 1; 29 } 30 return x; 31 } 32 33 inline void retreat(int x, int y) 34 { 35 int i; 36 for(i = x; i <= y; i++) b[i] = a[i]; 37 std::sort(b + x, b + y + 1); 38 } 39 40 inline void init() 41 { 42 int i, j; 43 S = int(sqrt(n)) + 10; 44 for(i = 1; i <= n; i++) scanf("%d", &a[i]); 45 for(i = 1; i <= n; i += S) 46 { 47 st[++C] = i; 48 ed[C] = std::min(i + S - 1, n); 49 } 50 for(i = 1; i <= C; i++) 51 for(j = st[i]; j <= ed[i]; j++) 52 belong[j] = i; 53 for(i = 1; i <= C; i++) retreat(st[i], ed[i]); 54 } 55 56 inline void update(int x, int y, LL z) 57 { 58 int i, l = belong[x], r = belong[y]; 59 if(l == r) 60 { 61 for(i = x; i <= y; i++) a[i] += z; 62 retreat(st[l], ed[r]); 63 } 64 else 65 { 66 for(i = x; i <= ed[l]; i++) a[i] += z; 67 for(i = l + 1; i <= r - 1; i++) add[i] += z; 68 for(i = st[r]; i <= y; i++) a[i] += z; 69 retreat(st[l], ed[l]); 70 retreat(st[r], ed[r]); 71 } 72 } 73 74 inline int query(int x, int y, LL z) 75 { 76 int i, l = belong[x], r = belong[y], ans = 0; 77 if(l == r) return y + 1 - find(x, y + 1, z - add[l]); 78 ans += ed[l] - find(x, ed[l] + 1, z - add[l]) + 1; 79 for(i = l + 1; i <= r - 1; i++) ans += ed[i] - find(st[i], ed[i] + 1, z - add[i]) + 1; 80 ans += y - find(st[r], y + 1, z - add[r]) + 1; 81 return ans; 82 } 83 84 int main() 85 { 86 int i, j, x, y; 87 LL z; 88 char ch; 89 n = read(); 90 q = read(); 91 init(); 92 for(i = 1; i <= q; i++) 93 { 94 while ((ch=getchar()) < 'A' || ch > 'Z'); 95 x = read(); 96 y = read(); 97 z = read(); 98 if(ch == 'M') update(x, y, z); 99 else printf("%d\n", query(x, y, z)); 100 } 101 return 0; 102 } View Code

?

轉載于:https://www.cnblogs.com/zhenghaotian/p/6823527.html

總結

以上是生活随笔為你收集整理的[luoguP2801] 教主的魔法（二分 + 分块）的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。