日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

【bzoj3343】教主的魔法 (分块 + 二分)

發布時間:2025/3/15 编程问答 37 豆豆
生活随笔 收集整理的這篇文章主要介紹了 【bzoj3343】教主的魔法 (分块 + 二分) 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

傳送門(權限題)

題目分析

題意為:給定一個數列,修改和查詢兩種操作,修改每次給定一個區間,區間的所有元素都加上一個給定值,查詢詢問一段區間的數權值大于等于給定值的數有多少個。 首先對原序列分塊,然后將分塊后的數組每個塊內的數字進行排序,這樣查詢時就可以暴力枚舉散塊,并二分枚舉每個整塊。對于修改,對于整塊部分只需修改tag,然后暴力修改散塊的原序列值,然后對散塊元素所在塊進行重排序即可。

code

3048?ms

#include<iostream> #include<cstdio> #include<cstdlib> #include<cstring> #include<cmath> #include<algorithm> #include<vector> using namespace std;const int N = 1e6 + 5; int n, Q, h[N], H[N], tag[N], S, blo; #define bug(x) cout<<#x<<":"<<x<<endl inline int read(){int i = 0, f = 1; char ch = getchar();for(; (ch < '0' || ch > '9') && ch != '-'; ch = getchar());if(ch == '-') f = -1, ch = getchar();for(; ch >= '0' && ch <= '9'; ch = getchar())i = (i << 3) + (i << 1) + (ch - '0');return i * f; }inline void wr(int x){if(x < 0) putchar('-'),x = -x;if(x > 9) wr(x / 10);putchar(x % 10 + '0'); }inline void add(int x, int y, int v){if(y - x + 1 <= 2 * S){for(int i = x; i <= y; i++)h[i] += v;}int Bx = x / S + (x % S ? 1 : 0), By = y / S + (y % S ? 1 : 0);int L = Bx + 1, R = By - 1;if(x == (Bx - 1) * S + 1) L--;if(y == min(n, By * S)) R++;int ans = 0;for(int i = x; i < (L - 1) * S; i++)h[i] += v;int l = (Bx - 1) * S + 1, r = min(n, Bx * S);for(int i = l; i <= r; i++)H[i] = h[i];sort(H + l, H + r + 1);for(int i = min(n, R * S) + 1; i <= y; i++)h[i] += v;l = (By - 1) * S + 1, r = min(n, By * S);for(int i = l; i <= r; i++)H[i] = h[i];sort(H + l, H + r + 1);for(int i = L; i <= R; i++)tag[i] += v; }inline int query(int x, int y, int v){int ans = 0;if(y - x + 1 <= 2 * S){for(int i = x; i <= y; i++)if(h[i] + tag[i / S + (i % S ? 1 : 0)] >= v) ans++;return ans;}int Bx = x / S + (x % S ? 1 : 0), By = y / S + (y % S ? 1 : 0); // bug(Bx), bug(By);int L = Bx + 1, R = By - 1;if(x == (Bx - 1) * S + 1) L--;if(y == min(n, By * S)) R++; // bug(L), bug(R);for(int i = x; i <= (L - 1) * S; i++)if(h[i] + tag[Bx] >= v) ans++;for(int i = min(n, R * S) + 1; i <= y; i++)if(h[i] + tag[By] >= v) ans++;for(int i = L; i <= R; i++){int l = (i - 1) * S + 1, r = min(n, i * S), len = r - l + 1;int tmp = lower_bound(H + l, H + r + 1, v - tag[i]) - (H + l - 1);ans += len - (tmp - 1);}return ans; }int main(){n = read(), Q = read(), S = sqrt(n), blo = n / S + (n % S ? 1 : 0); // bug(S);for(int i = 1; i <= n; i++) h[i] = H[i] = read();for(int i = 1; i <= blo; i++){int l = (i - 1) * S + 1, r = min(i * S, n);sort(H + l, H + r + 1);}for(int i = 1; i <= Q; i++){char opt[5]; scanf("%s", opt + 1);if(opt[1] == 'M'){int l = read(), r = read(), w = read();add(l, r, w);}else if(opt[1] == 'A'){int l = read(), r = read(), c = read();wr(query(l, r, c)), putchar('\n');}}return 0; }

?

轉載于:https://www.cnblogs.com/CzYoL/p/7396352.html

總結

以上是生活随笔為你收集整理的【bzoj3343】教主的魔法 (分块 + 二分)的全部內容,希望文章能夠幫你解決所遇到的問題。

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