[BZOJ]3110 K大数查询(ZJOI2013)
这大概是唯一一道小C重写了4次的题目。
姿势不对的树套树(Fail) → 分块(Fail) → 整体二分(Succeed) → 树套树(Succeed)。
让小C写点心得静静。
Description
有N个位置,M个操作。操作有两种,每次操作如果是1 a b c的形式表示在第a个位置到第b个位置,每个位置加入一个数c。如果是2 a b c形式,表示询问从第a个位置到第b个位置,第C大的数是多少。
Input
第一行N,M。
接下来M行,每行形如1 a b c或2 a b c。
Output
输出每个询问的结果
Sample Input
2 5
1 1 2 1
1 1 2 2
2 1 1 2
2 1 1 1
2 1 2 3
Sample Output
1
2
1
HINT
N,M<=50000,N,M<=50000
a<=b<=N
1操作中abs(c)<=N
2操作中c<=Maxlongint
Solution
这道题让小C见识了整体二分的妙用,它是一种离线做法。
整体二分的思想就是二分答案。以单次询问的区间K大为例,我们可以二分答案,通过统计比当前二分出的答案大的数的个数,与K比较,决定缩小或是增大答案。
而对于多次询问,我们整体二分所有询问的答案,对于每个询问,我们都通过比较,来决定它们应该进入小的那一半,还是大的那一半。
而对于每个修改,我们也要判断它会对哪个区间的答案产生影响,来决定它应该进入哪一半。
要注意,任何时候,每个询问和修改之间的相对顺序都不能变。
这样的时间复杂度为什么是科学的呢?我们脑补一下:
动态询问的本质就是边修改边询问的过程,对于每次修改,我们用数据结构维护要统计的信息。(这是一句废话)
因为有数据结构的存在,时间复杂度才是科学的。(还不理解的读者请将这一行无视,继续往下看)
What?还需要数据结构?那我要整体二分何用?我咋不大力树套树呢?
整体二分当然不是白写的。原来树套树要维护的是,值域区间[L,R]内,数组的信息,需要二维。
而整体二分时,我们只需要维护 当前 二分的值域区间[L,R]内,数组的信息,只需要一维。
因为做完 当前 这段值域区间内的事情后,这棵线段树/平衡树就彻底没用了,因为关于它的所有事情都做完了,清空就行。(理解了再往下看吧)
这就巧妙利用离线做法将数据结构降维,空间和时间(常数)上都有很大的提升。
最后证明一下时间复杂度:对于每个询问和修改,我们都在二分答案的 log 个区间内做了一遍,总时间复杂度。
当然具体实现还有一些细节,还不理解的可以参见小C的代码。
说说这道题的题解:
树套树:以值域为第一维,以数组下标为第二维,第二维维护区间和。对于每个询问,像二分答案一样在第一维上走到底。
整体二分:离线,把树套树做法降维。
整体二分:
#include <cstdio>
#include <cstring>
#include <algorithm>
#define l(a) son[a][0]
#define r(a) son[a][1]
#define MN 50005
#define ll long long
using namespace std;
struct meg{int g,x,y,z,pos,gs;}q[MN],q1[MN],q2[MN];
int son[MN<<][],tg[MN<<],din,rt;
int ans[MN],n,m,ain;
ll t[MN<<]; inline int read()
{
int n=,f=; char c=getchar();
while (c<'' || c>'') {if(c=='-')f=-; c=getchar();}
while (c>='' && c<='') {n=n*+c-''; c=getchar();}
return n*f;
} inline void mark(int &x,int L,int R,int z)
{
if (!x) {x=++din; l(x)=r(x)=t[x]=tg[x]=;}
t[x]+=1LL*z*(R-L+); tg[x]+=z;
}
inline void down(int x,int L,int R)
{
if (!tg[x]) return;
int mid=L+R>>;
mark(l(x),L,mid,tg[x]); mark(r(x),mid+,R,tg[x]);
tg[x]=;
} ll getsum(int x,int L,int R,int ql,int qr)
{
if (!x) return ;
if (ql==L&&qr==R) return t[x];
down(x,L,R);
int mid=L+R>>;
if (qr<=mid) return getsum(l(x),L,mid,ql,qr);
else if (ql>mid) return getsum(r(x),mid+,R,ql,qr);
else return getsum(l(x),L,mid,ql,mid)+getsum(r(x),mid+,R,mid+,qr);
}
void add(int &x,int L,int R,int ql,int qr)
{
if (!x) {x=++din; l(x)=r(x)=t[x]=tg[x]=;}
if (ql==L&&qr==R) {mark(x,ql,qr,); return;}
down(x,L,R);
int mid=L+R>>;
if (qr<=mid) add(l(x),L,mid,ql,qr);
else if (ql>mid) add(r(x),mid+,R,ql,qr);
else add(l(x),L,mid,ql,mid),add(r(x),mid+,R,mid+,qr);
t[x]=t[l(x)]+t[r(x)];
} void work(int L,int R,int hd,int tl)
{
if (hd>tl) return;
register int i;
if (L==R) {for (i=hd;i<=tl;++i) if (q[i].g==) ans[q[i].pos]=L; return;}
int mid=L+R>>,p1=,p2=;
ll lt;
rt=din=;
for (i=hd;i<=tl;++i)
{
if (q[i].g==)
{
if (q[i].z>mid) add(rt,,n,q[i].x,q[i].y),q2[++p2]=q[i];
else q1[++p1]=q[i];
}
else if (q[i].g==)
{
lt=getsum(rt,,n,q[i].x,q[i].y);
if (lt+q[i].gs>=q[i].z) q2[++p2]=q[i];
else q[i].gs+=lt,q1[++p1]=q[i];
}
}
for (i=;i<=p1;++i) q[hd+i-]=q1[i];
for (i=;i<=p2;++i) q[hd+p1+i-]=q2[i];
work(L,mid,hd,hd+p1-); work(mid+,R,hd+p1,tl);
} int main()
{
register int i;
n=read(); m=read();
for (i=;i<=m;++i) {q[i].g=read(); q[i].x=read(); q[i].y=read(); q[i].z=read(); if (q[i].g==) q[i].pos=++ain;}
work(-n,n,,m);
for (i=;i<=ain;++i) printf("%d\n",ans[i]);
}
树套树:
#include <cstdio>
#include <cstring>
#include <algorithm>
#define l(a) (a<<1)
#define r(a) (a<<1|1)
#define L(a) (son[a][0])
#define R(a) (son[a][1])
#define TMN 15000005
#define MN 50005
#define ll long long
using namespace std;
ll t[TMN];
int son[TMN][],tg[TMN],rt[MN<<],din;
int n,m; inline int read()
{
int n=,f=; char c=getchar();
while (c<'' || c>'') {if(c=='-')f=-; c=getchar();}
while (c>='' && c<='') {n=n*+c-''; c=getchar();}
return n*f;
} ll Tgetsum(int x,int L,int R,int ql,int qr)
{
if (!x) return ;
if (ql==L&&qr==R) return t[x];
int mid=L+R>>;
ll lt=1LL*tg[x]*(qr-ql+);
if (qr<=mid) return lt+Tgetsum(L(x),L,mid,ql,qr);
else if (ql>mid) return lt+Tgetsum(R(x),mid+,R,ql,qr);
else return lt+Tgetsum(L(x),L,mid,ql,mid)+Tgetsum(R(x),mid+,R,mid+,qr);
}
void Tadd(int &x,int L,int R,int ql,int qr)
{
if (!x) x=++din;
if (ql==L&&qr==R) {++tg[x]; t[x]+=(ll)R-L+; return;}
int mid=L+R>>;
if (qr<=mid) Tadd(L(x),L,mid,ql,qr);
else if (ql>mid) Tadd(R(x),mid+,R,ql,qr);
else Tadd(L(x),L,mid,ql,mid),Tadd(R(x),mid+,R,mid+,qr);
t[x]=t[L(x)]+t[R(x)]+1LL*tg[x]*(R-L+);
} int getk(int x,int L,int R,int z,int yl,int yr)
{
if (L==R) return L;
ll lt;
int mid=L+R>>;
if ((lt=Tgetsum(rt[r(x)],,n,yl,yr))<z) return getk(l(x),L,mid,z-lt,yl,yr);
else return getk(r(x),mid+,R,z,yl,yr);
}
void add(int x,int L,int R,int q,int yl,int yr)
{
Tadd(rt[x],,n,yl,yr);
if (L==R) return;
int mid=L+R>>;
if (q<=mid) add(l(x),L,mid,q,yl,yr);
else add(r(x),mid+,R,q,yl,yr);
} int main()
{
register int g,x,y,z;
n=read(); m=read();
while (m--)
{
g=read(); x=read(); y=read(); z=read();
if (g==) add(,-n,n,z,x,y);
else if (g==) printf("%d\n",getk(,-n,n,z,x,y));
}
}
Last Word
说说前两次打这题的惨痛经历吧:
第一次树套树把(相对于正解)第一维和第二维搞反了,导致空间爆炸。
第二次分块,脑补了一个时间和空间都是O(n*sqrt(n)*log(n))的做法(权值线段树+替罪羊树瞎写),空间还是爆炸。
两次做法的思路都是二分答案,而且都把数组下标看作是第一维。
不过所幸打的这几次都没怎么需要Debug。打得挺爽。
[BZOJ]3110 K大数查询(ZJOI2013)的更多相关文章
- BZOJ 3110 K大数查询 | 整体二分
BZOJ 3110 K大数查询 题面 有N个位置,M个操作.操作有两种,每次操作如果是1 a b c的形式表示在第a个位置到第b个位置,每个位置加入一个数c 如果是2 a b c形式,表示询问从第a个 ...
- BZOJ 3110 k大数查询 & 树套树
题意: 有n个位置,每个位置可以看做一个集合,现在要求你实现一个数据结构支持以下功能: 1:在a-b的集合中插入一个数 2:询问a-b集合中所有元素的第k大. SOL: 调得火大! 李建说数据结构题能 ...
- bzoj 3110 K大数查询
第一道整体二分,因为只需要知道每个询问区间中比mid大的数有多少个,就可以直接用线段树区间加,区间求和了. #include<iostream> #include<cstdio> ...
- BZOJ 3110: [Zjoi2013]K大数查询 [树套树]
3110: [Zjoi2013]K大数查询 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 6050 Solved: 2007[Submit][Sta ...
- 树套树专题——bzoj 3110: [Zjoi2013] K大数查询 & 3236 [Ahoi2013] 作业 题解
[原题1] 3110: [Zjoi2013]K大数查询 Time Limit: 20 Sec Memory Limit: 512 MB Submit: 978 Solved: 476 Descri ...
- bzoj 3110: [Zjoi2013]K大数查询 树状数组套线段树
3110: [Zjoi2013]K大数查询 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 1384 Solved: 629[Submit][Stat ...
- BZOJ 3110: [Zjoi2013]K大数查询( 树状数组套主席树 )
BIT+(可持久化)权值线段树, 用到了BIT的差分技巧. 时间复杂度O(Nlog^2(N)) ---------------------------------------------------- ...
- BZOJ 3110([Zjoi2013]K大数查询-区间第k大[段修改,在线]-树状数组套函数式线段树)
3110: [Zjoi2013]K大数查询 Time Limit: 20 Sec Memory Limit: 512 MB Submit: 418 Solved: 235 [ Submit][ ...
- BZOJ 3110 [Zjoi2013]K大数查询(整体二分)
3110: [Zjoi2013]K大数查询 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 11654 Solved: 3505[Submit][St ...
随机推荐
- [USACO13JAN] Seating
https://www.luogu.org/problem/show?pid=3071 题目描述 To earn some extra money, the cows have opened a re ...
- Linux - IDA - 安装 ( 带F5功能 )
Linux - IDA - 安装 ( 带F5功能 ) 0x00 前言 最近在熟悉deepin系统,想着把逆向的一些软件也迁移过去,但像ida,Ollydbg这些工具一般都是在windows下使用,所以 ...
- bootstrap的ajax提交
一般后台界面都用bootstrap框架,这是一个css框架,里面封装了ajax方法,只需要在样式中指定就行,根本自己不用写 <td> <eq name='item.status' v ...
- maven安装、配置
maven的安装和配置 1.将maven解压到自定义文件夹下.例如解压到如下目录(解压目录最好不要有中文字): 2:配置环境变量:一定要注意要用分号:与其他值隔开 3.在cmd中测试,验证是否安装成 ...
- linux下的Shell编程(6)case和select
第一个,除了if语句之外,Shell Script中也有类似C语言中多分支结构的case语句,它的语法是: case var in pattern 1 ) - ;; pattern 2 ) - ;; ...
- Linux将端口设置进防火墙的白名单
1.先检查linux服务器的端口是否被防火墙拦住 `telnet 172.168.1.101 8080后面跟端口号,如果连接上证明是防火墙白名单.如果没有配置 vi /etc/sysconfig/ip ...
- vue组件详解(五)——组件高级用法
一.递归组件 组件在它的模板内可以递归地调用自己, 只要给组件设置name 的选项就可以了. 示例如下: <div id="app19"> <my-compone ...
- 深度学习中Xavier初始化
"Xavier"初始化方法是一种很有效的神经网络初始化方法,方法来源于2010年的一篇论文<Understanding the difficulty of training ...
- PHP+MySQL分页显示示例分析
Web开发是今后分布式程式开发的主流,通常的web开发都要涉及到与数据库打交道,客户端从服务器端读取通常都是以分页的形式来显示,一页一页的阅读起来既方便又美观.所以说写分页程序是web开发的一个重要组 ...
- Lintcode375 Clone Binary Tree solution 题解
[题目描述] For the given binary tree, return a deep copy of it. 深度复制一个二叉树,给定一个二叉树,返回一个他的克隆品. [题目链接] www. ...