整体二分初探 两类区间第K大问题 poj2104 & hdu5412
看到好多讲解都把整体二分和$CDQ$分治放到一起讲 不过自己目前还没学会$CDQ$分治 就单独谈谈整体二分好了
先推荐一下$XHR$的 <浅谈数据结构题的几个非经典解法> 整体二分在当中有较为详细的讲解
先来说一下静态第$K$小的整体二分解法 $(POJ2104)$
题目链接:http://poj.org/problem?id=2104
所谓整体二分 就是利用所有的询问相互独立而把它们$($此题没有修改操作$)$通过二分把它们分治进行处理
不妨先来考虑下一个简单易懂的$O(NlogS)$的排序算法$(S$为数值范围$)$
这个方法是自己在思考整体二分的时候$yy$的 虽然在实际应用上没什么意义 但是有助于理解整体二分的分治过程
我们假设当前处理的区间里最小值不小于$L$ 最大值不大于$R$ 令$MID = (L + R) / 2$
然后把当前区间扫描一遍 如果一个数不大于$MID$就放到左子区间 否则放到右子区间
如此下去 直到区间内只剩一个数或者$L$ 与 $R$相等 排序就完成了
现在回到静态区间第$K$小问题 和刚才那个排序算法类似 我们先二分一个答案$MID$
如果区间内小于等于$MID$的数的个数$($记为$CNT)$不超过$K$ 那么最终答案显然也是不超过$MID$的
这类询问我们把它们划分到左子区间
而对于$CNT$大于$K$的 我们则把它们划分到右子区间 并且把$K$减去$CNT$
换句话说就是把小于等于$MID$的数的贡献全部算上后之后就不用考虑了
可以发现这样划分的层数是$logS$ 而每一层的询问个数是$Q$个 再加上算贡献时用到的$BIT$ 所以复杂度是$O(QlogNlogS)$
以下是$poj2104$参考代码
//#include <bits/stdc++.h>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
const int N = 1e5 + , Q = 5e3 + , lim = 1e9;
struct point
{
int x, num;
}a[N];
struct query
{
int x, y, k, cnt, num;
}q[Q], b[Q];
int sum[N], ans[Q];
int n, m;
bool cmp(const point &aa, const point &bb)
{
return aa.x < bb.x;
}
void calc(int ll, int rr, int rawL, int mid)
{
int L = , R = n + , MID;
while(L < R)
{
MID = (L + R) >> ;
if(a[MID].x >= rawL)
R = MID;
else
L = MID + ;
}
for(int i = R; i <= n && a[i].x <= mid; ++i)
for(int j = a[i].num; j <= n; j += (j & -j))
++sum[j];
for(int i = ll; i <= rr; ++i)
{
q[i].cnt = ;
for(int j = q[i].y; j; j -= (j & -j))
q[i].cnt += sum[j];
for(int j = q[i].x - ; j; j -= (j & -j))
q[i].cnt -= sum[j];
}
for(int i = R; i <= n && a[i].x <= mid; ++i)
for(int j = a[i].num; j <= n; j += (j & -j))
--sum[j];
}
void divide(int ll, int rr, int rawL, int rawR)
{
if(rawL == rawR)
{
for(int i = ll; i <= rr; ++i)
ans[q[i].num] = rawR;
return;
}
int mid = rawL + ((rawR - rawL) >> );
calc(ll, rr, rawL, mid);
int now1 = ll, now2 = rr;
for(int i = ll; i <= rr; ++i)
{
if(q[i].cnt >= q[i].k)
b[now1++] = q[i];
else
{
q[i].k -= q[i].cnt;
b[now2--] = q[i];
}
}
for(int i = ll; i <= rr; ++i)
q[i] = b[i];
if(now1 != ll)
divide(ll, now1 - , rawL, mid);
if(now2 != rr)
divide(now2 + , rr, mid + , rawR);
}
int main()
{
scanf("%d%d", &n, &m);
for(int i = ; i <= n; ++i)
{
scanf("%d", &a[i].x);
a[i].num = i;
}
a[n + ].x = 2e9;
sort(a + , a + + n, cmp);
for(int i = ; i <= m; ++i)
{
scanf("%d%d%d", &q[i].x, &q[i].y, &q[i].k);
q[i].num = i;
}
divide(, m, -lim, lim);
for(int i = ; i <= m; ++i)
printf("%d\n", ans[i]);
return ;
}
然后就是整体二分真正突出存在意义的问题 动态区间第$K$大了$(hdu5412)$
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5412
动态区间第$K$大就是额外多了修改操作 修改不仅对询问有影响 对修改也有影响 因此直观上看起来要麻烦很多
我们知道 区间询问之所以能排序后分治 主要是所有的询问相互独立
对于修改操作 如果我们希望它们相互独立该怎么做呢 那就算贡献好了
实际上 每次把一个位置的数修改为另一个数相当于使这个位置删除一个原来的数 并插入一个新的数 而初始的数则相当于仅有插入操作
这样只要保证同一区间内的查询和删除或插入操作的相对顺序不变 即相对时间顺序不变 最后算贡献结果也是一样的
并且由于对大于$MID$的数的删除或插入操作 即便时间顺序是在划分到左子区间的询问之前 也不会造成影响 因此它们可以划分到右子区间
而对小于$MID$的数的删除或插入操作 对划分到右子区间的询问的影响可以直接在划分前计入贡献之后不再考虑 因此它们可以划分到左子区间
这样按照类似静态区间第$K$大问题计算复杂度会发现是$O((Q+N)logNlogS)$
以下是$hdu5412$的参考代码
#include <bits/stdc++.h>
using namespace std;
const int N = 1e5 + , Q = 1e5 + , A = N + Q * , lim = 1e9;
struct operate
{
int x, y, k, cnt, num;
operate(){}
operate(int _x, int _y, int _k, int _cnt, int _num)
{
x = _x;
y = _y;
k = _k;
cnt = _cnt;
num = _num;
}
}a[A], q1[A], q2[A];
int raw[N], ans[Q], sum[N];
int n, m, len, l1, l2;
void update(int x, int y)
{
for(int i = x; i <= n; i += (i & -i))
sum[i] += y;
}
int query(int x)
{
int re = ;
for(int i = x; i; i -= (i & -i))
re += sum[i];
return re;
}
void calc(int ll, int rr, int rawl, int mid)
{
for(int i = ll; i <= rr; ++i)
{
if(a[i].k)
a[i].cnt = query(a[i].y) - query(a[i].x - );
else if(a[i].y <= mid)
update(a[i].x, a[i].cnt);
}
for(int i = ll; i <= rr; ++i)
if(!a[i].k && a[i].y <= mid)
update(a[i].x, -a[i].cnt);
l1 = l2 = ;
for(int i = ll; i <= rr; ++i)
if(a[i].k)
{
if(a[i].k <= a[i].cnt)
q1[++l1] = a[i];
else
{
a[i].k -= a[i].cnt;
q2[++l2] = a[i];
}
}
else
{
if(a[i].y <= mid)
q1[++l1] = a[i];
else
q2[++l2] = a[i];
}
int now = ll;
for(int i = ; i <= l1; ++i)
a[now++] = q1[i];
for(int i = ; i <= l2; ++i)
a[now++] = q2[i];
}
void divide(int ll, int rr, int rawl, int rawr)
{
if(rawl == rawr)
{
for(int i = ll; i <= rr; ++i)
if(a[i].k)
ans[a[i].num] = rawl;
return;
}
int mid = (rawl + rawr) >> ;
calc(ll, rr, rawl, mid);
int tmp = l1;
if(tmp)
divide(ll, ll + tmp - , rawl, mid);
if(ll + tmp <= rr)
divide(ll + tmp, rr, mid + , rawr);
}
int main()
{
while(scanf("%d", &n) != EOF)
{
len = ;
for(int i = ; i <= n; ++i)
{
scanf("%d", &raw[i]);
a[++len] = operate(i, raw[i], , , );
}
scanf("%d", &m);
int op, x, y, z;
for(int i = ; i <= m; ++i)
{
scanf("%d", &op);
if(op & )
{
scanf("%d%d", &x, &y);
a[++len] = operate(x, raw[x], , -, );
a[++len] = operate(x, y, , , );
raw[x] = y;
ans[i] = ;
}
else
{
scanf("%d%d%d", &x, &y, &z);
a[++len] = operate(x, y, z, , i);
}
}
divide(, len, , lim);
for(int i = ; i <= m; ++i)
if(ans[i])
printf("%d\n", ans[i]);
}
return ;
}
整体二分初探 两类区间第K大问题 poj2104 & hdu5412的更多相关文章
- 整体二分(模板) 求区间第k小
整体二分,将询问与初值一起放入一个结构体里,然后每次二分判断询问在哪边,树状数组维护,时间复杂度O((n+Q)lognlogMAX_a[i] 代码 #include<iostream> # ...
- 【BZOJ3110】【整体二分+树状数组区间修改/线段树】K大数查询
Description 有N个位置,M个操作.操作有两种,每次操作如果是1 a b c的形式表示在第a个位置到第b个位置,每个位置加入一个数c 如果是2 a b c形式,表示询问从第a个位置到第b个位 ...
- 【bzoj3110】[Zjoi2013]K大数查询 整体二分+树状数组区间修改
题目描述 有N个位置,M个操作.操作有两种,每次操作如果是1 a b c的形式表示在第a个位置到第b个位置,每个位置加入一个数c.如果是2 a b c形式,表示询问从第a个位置到第b个位置,第C大的数 ...
- [NBUT 1458 Teemo]区间第k大问题,划分树
裸的区间第k大问题,划分树搞起. #pragma comment(linker, "/STACK:10240000") #include <map> #include ...
- 区间第k大问题 权值线段树 hdu 5249
先说下权值线段树的概念吧 权值平均树 就是指区间维护值为这个区间内点出现次数和的线段树 用这个加权线段树 解决第k大问题就很方便了 int query(int l,int r,int rt,int k ...
- 算法笔记--CDQ分治 && 整体二分
参考:https://www.luogu.org/blog/Owencodeisking/post-xue-xi-bi-ji-cdq-fen-zhi-hu-zheng-ti-er-fen 前置技能:树 ...
- Dynamic Rankings(整体二分)
Dynamic Rankings(整体二分) 带修区间第k小.\(n,q\le 10^4\). 这次我们旧瓶装新酒,不用带修主席树.我们用整体二分!整体二分是啥东西呢? 二分答案可以解决一次询问的问题 ...
- BZOJ_3110_[Zjoi2013]K大数查询_整体二分+树状数组
BZOJ_3110_[Zjoi2013]K大数查询_整体二分+树状数组 Description 有N个位置,M个操作.操作有两种,每次操作如果是1 a b c的形式表示在第a个位置到第b个位置,每个位 ...
- 一篇自己都看不懂的CDQ分治&整体二分学习笔记
作为一个永不咕咕咕的博主,我来更笔记辣qaq CDQ分治 CDQ分治的思想还是比较简单的.它的基本流程是: \(1.\)将所有修改操作和查询操作按照时间顺序并在一起,形成一段序列.显然,会影响查询操作 ...
随机推荐
- 【JZOJ 3918】蛋糕
题面: 正文: 根据数据\(4\leq R,C\leq 75\)我们大概可以先枚举切横的刀再二分答案. 更具体的: 假设我们已经枚举到这样横切: 再假设我们已经二分到最小的巧克力是\(7\) 康康第一 ...
- P2523 [HAOI2011]Problem c
传送门 先考虑如何判断无解,设 $sum[i]$ 表示确定的人中,编号大于 $i$ 的人的人数 如果 $sum[i]>n-i+1$ 则无解,进一步考虑设 $f[i][j]$ 表示当前确定完编号大 ...
- 在django中使用循环与条件语言
{% if not Article_type_id %} <li class="active"><a href="/">全部</a ...
- 什么情况下会出现undefined
1.函数定义形参不传值,2.预解释,只声明不定义时输出变量3.对象取属性值,属性值不存在
- 神经网络中的反向传播法--bp【转载】
from: 作者:Charlotte77 出处:http://www.cnblogs.com/charlotte77/ 一文弄懂神经网络中的反向传播法——BackPropagation 最近在看深度学 ...
- 负载均衡的3种模型(httpd+lvs几十万并发的负载均衡搭建)
一.几种常见的负载均衡模型 二.搭建httpd+lvs LVS .node01 启动一块eth0:2的网卡子接口(ifconfig eth0: down可以把网卡down掉) ifconfig eth ...
- 用element-ui进行文件的上传
在进项的项目中,我使用vue配合element搭建的项目,现在需要用到上传文件的写法!现在列举下我的项目中所使用的上传文件的方法! <el-upload style="display: ...
- diff 比较两个文件的不同
1.命令功能 diff 逐行比较文件内容,并输出文件差异. 2.语法格式 diff option file1 file2 diff 选项 文件1 文件2 参数说明 参数 参数说明 ...
- 前端每日实战:14# 视频演示如何用纯 CSS 创作一种侧立图书的特效
效果预览 按下右侧的"点击预览"按钮在当前页面预览,点击链接全屏预览. https://codepen.io/zhang-ou/pen/deVgRM 可交互视频教程 此视频是可以交 ...
- java:序列化Serializable 接口
java:序列化Serializable 接口 public class SerializePerson implements Serializable { private String name; ...