【算法学习】整体二分

我们开门见山，讲讲一道sb题：

给你一个数组，查这个数组的第x大元素。

排序？可以

二分？怎么做啊？

二分出一个mid，判断这个数组中有多少个数小于等于mid，如果个数大于等于x，就递归到[l,mid]区间，否则是[mid+1,r]区间，这样递归下去就能得到结果。

怎么计算小于等于mid的个数？

for一遍原数组？其实只要存储数组中值在[l,r]中的数，然后计算就好，而对于[1,l-1]，可以预先存储。

怎么计算？树状数组，范围太大就离散化。

这样很麻烦吧？排序更好吧？

但是当我给你多个询问时，你就不会这么觉得了。

回顾刚才的过程，我们做上图的类比。

把询问看成球，答案的区间就是这个矩形，每一次我们将上一层的球判断它会掉到下层的哪一块中，掉到底层就是获得了答案。

刚刚的过程就是球只有1个的情况。

球有多个时，就是整体二分！

就是对于当前的mid，对每一个询问都计算答案是大了还是小了，然后决定分配到下一层的哪边。

要注意保证时间复杂度哦！

一道例题：HDU 5412。

查询区间K小，要支持单点修改。

看代码：

 1 #include<cstdio>
 2 int n,q,a[100001],sum[300001],ans[300001],que1[300001],que2[300001],bit[300001];
 3 int I[300001],type[300001],ql[300001],qr[300001],k[300001],cnt;
 4 inline void ins(int t,int l,int r,int x){type[++cnt]=t,ql[cnt]=l,qr[cnt]=r,k[cnt]=x,I[cnt]=cnt;}
 5 inline void Ins(int i,int x){for(;i<=n;bit[i]+=x,i+=i&-i);}
 6 inline int Qur(int i){int Sum=0;for(;i;Sum+=bit[i],i-=i&-i);return Sum;}
 7 void divide(int l,int r,int low,int upp){
 8     if(l>r) return;
 9     if(low==upp) {for(int i=l;i<=r;++i) if(type[I[i]]==3) ans[I[i]]=low; return;}
10     int mid=low+upp>>1,cl=0,cr=0;
11     for(int i_=l,i;i_<=r;++i_){
12         i=I[i_];
13         if(type[i]==1) {if(k[i]<=mid) Ins(ql[i],1), que1[++cl]=i; else que2[++cr]=i;}
14         if(type[i]==2) {if(k[i]<=mid) Ins(ql[i],-1), que1[++cl]=i; else que2[++cr]=i;}
15         if(type[i]==3){
16             int tmp=Qur(qr[i])-Qur(ql[i]-1);
17             if(k[i]<=sum[i]+tmp) que1[++cl]=i;
18             else que2[++cr]=i, sum[i]+=tmp;
19         }
20     }
21     for(int i_=l,i;i_<=r;++i_){
22         i=I[i_];
23         if(type[i]==1) if(k[i]<=mid) Ins(ql[i],-1);
24         if(type[i]==2) if(k[i]<=mid) Ins(ql[i],1);
25     }
26     for(int i=1;i<=cl;++i) I[l+i-1]=que1[i];
27     for(int i=1;i<=cr;++i) I[l+cl+i-1]=que2[i];
28     divide(l,l+cl-1,low,mid); divide(l+cl,r,mid+1,upp);
29 }
30 int main(){
31     scanf("%d",&n);
32     for(int i=1;i<=n;++i) scanf("%d",a+i), ins(1,i,i,a[i]);
33     scanf("%d",&q);
34     for(int i=1,t,x,y,z;i<=q;++i){
35         scanf("%d",&t);
36         if(t==1) scanf("%d%d",&x,&y), ins(2,x,x,a[x]), ins(1,x,x,y), a[x]=y;
37         else scanf("%d%d%d",&x,&y,&z), ins(3,x,y,z);
38     }
39     divide(1,cnt,1,1000000000);
40     for(int i=1;i<=cnt;++i) if(ans[i]) printf("%d\n",ans[i]);
41     return 0;
42 }

View Code

一道例题：洛谷 P3332。

查询区间K大，要支持区间加入。

看代码：

 1 #include <cstdio>
 2
 3 typedef long long LL;
 4 const int MN = 50005;
 5
 6 int N, Q;
 7
 8 int opt[MN], L[MN], R[MN]; LL V[MN];
 9 int P[MN];
10 int s1[MN], s2[MN], t1, t2;
11
12 LL Sum[MN];
13 int Ans[MN];
14
15 LL b1[MN], b2[MN];
16 inline void Add(LL *b, int i, LL x) { for(; i <= N; i += i & -i) b[i] += x; }
17 inline LL Qur(LL *b, int i) { LL A = 0; for(; i; i -= i & -i) A += b[i]; return A; }
18 inline void Add(int l, int r, LL x) {
19     ++r;
20     Add(b1, l,  x), Add(b2, l,  (l - 1) * x);
21     Add(b1, r, -x), Add(b2, r, -(r - 1) * x);
22 }
23 inline LL Qur(int l, int r) {
24     --l;
25     return r * Qur(b1, r) - Qur(b2, r) - l * Qur(b1, l) + Qur(b2, l);
26 }
27
28 void Solve(int l, int r, int lb, int rb) {
29     if (lb == rb) {
30         for (int i = l; i <= r; ++i)
31             if (opt[P[i]] == 2)
32                 Ans[P[i]] = lb;
33         return ;
34     }
35     int mid = lb + rb >> 1;
36     t1 = t2 = 0;
37     for (int i = l; i <= r; ++i) {
38         if (opt[P[i]] == 1) {
39             if (V[P[i]] <= mid)
40                 s1[++t1] = P[i];
41             else {
42                 Add(L[P[i]], R[P[i]], 1);
43                 s2[++t2] = P[i];
44             }
45         }
46         else {
47             LL num = Sum[P[i]] + Qur(L[P[i]], R[P[i]]);
48             if (num >= V[P[i]])
49                 s2[++t2] = P[i];
50             else {
51                 Sum[P[i]] = num;
52                 s1[++t1] = P[i];
53             }
54         }
55     }
56     for (int i = l; i <= r; ++i) {
57         if (opt[P[i]] == 2) continue;
58         if (V[P[i]] > mid)
59             Add(L[P[i]], R[P[i]], -1);
60     }
61     int pos = l;
62     for (int i = 1; i <= t1; ++i)
63         P[pos++] = s1[i];
64     for (int i = 1; i <= t2; ++i)
65         P[pos++] = s2[i];
66     int M = l + t1 - 1;
67     Solve(l, M, lb, mid);
68     Solve(M + 1, r, mid + 1, rb);
69 }
70
71 int main() {
72     scanf("%d%d", &N, &Q);
73     for (int i = 1; i <= Q; ++i) {
74         scanf("%d%d%d%lld", opt + i, L + i, R + i, V + i);
75         P[i] = i;
76     }
77     Solve(1, Q, -N, N);
78     for (int i = 1; i <= Q; ++i) {
79         if (opt[i] == 1) continue;
80         printf("%d\n", Ans[i]);
81     }
82     return 0;
83 }

View Code

转载于:https://www.cnblogs.com/PinkRabbit/p/8039083.html

【算法学习】整体二分相关推荐

经典题:poj2104-区间第k小整体二分学习
写在前面区间第k小可以说是一个很经典的数据结构题了,这道题有很多种解法比如莫队离线.主席树.整体二分等等. 之前用莫队和主席树写过这道题,今天来学习一个以前不会的算法--整体二分. 因为最近遇到一 ...
数据结构与算法学习⑤(BFS和DFS 贪心算法二分查找)
数据结构与算法学习⑤ 数据结构与算法学习⑤ 1.BFS和DFS 1.1.深度优先搜索算法 1.2.广度优先搜索算法面试实战 102. 二叉树的层序遍历 104. 二叉树的最大深度 515. 在每个树 ...
【学习笔记】整体二分
文章目录引整体二分几道模板题 Dynamic Rankings [ZJOI2013]K大数查询 [国家集训队]矩阵乘法 [THUPC2017] 天天爱射击 [CTSC2018]混合果汁引例1 ...
【cdq分治】cdq分治与整体二分学习笔记Part2.cdq分治
上午的学习学会了整体二分,下午学了cdq分治发现了二者的区别: 整体二分的主体是在不断地二分答案(把所有询问二分),而cdq分治则是在不断地二分操作. 当然同样的,cdq分治的复杂度也是与区间长度正 ...
【cdq分治】cdq分治与整体二分学习笔记Part1.整体二分
之所以把cdq分治和整体二分放在一起学习,是因为他们两个实在太像了-不管是做法还是代码- 感觉整体二分可能会比cdq分治稍微简单那么一点点?所以先学整体二分. 整体二分是对答案进行二分,其具体操作如下 ...
数据结构和算法学习指南
点击上方蓝字设为星标下面开始今天的学习- 这篇文章会涵盖之前的所有内容,并且会举很多代码的实例,谈谈如何使用框架思维,并且给对于算法无从下手的朋友给一点具体可执行的刷题建议. 首先,这里讲的都是普通 ...
[总结]CDQ分治整体二分
从昨天到现在除了90%的颓废时间一直在研究一些分治的姿势,主要就是CDQ分治和整体二分. 首先推荐一些学习资料: 陈丹琦 <从 < Cash > 谈一类分治算法的应用> 许昊然 ...
算法学习之道，应有三重境界
https://www.toutiao.com/a6712297555167805966/ 王国维先生在<人间词话>中写道:古今之成大事业.大学问者,必经过三种境界:"昨夜西风凋 ...
C++初级算法-学习笔记
目录算法的基本概念什么是算法什么是数据结构模拟与高精度模拟算法高精度运算算法评价与算法复杂度算法"评价员" 如何评价算法暴力枚举枚举子集枚举排列枚举排序 ...

【算法学习】整体二分

【算法学习】整体二分相关推荐

最新文章

热门文章