BZOJ.3110.[ZJOI2013]K大数查询(整体二分 树状数组/线段树)
整体二分求的是第K小(利用树状数组)。求第K大可以转为求第\(n-K+1\)小，但是这样好像得求一个\(n\)。
注意到所有数的绝对值\(\leq N\)，将所有数的大小关系反过来 第\(K\)大就是第\(K\)小了。所有数\(A[i]\)改为\(n-A[i]\)，输出的时候也改为\(n-Ans[i]\)。
区间加入一个数\(C\)可以直接用线段树区间加，也可以树状数组维护常数会小很多。
\(n*m=2.5*1e9 > MAX\_INT\) 也是没谁了。。
现学了一波树状数组区间修改区间查询...->见这里。
洛谷Rank1在BZOJ上60多。。
树状数组：
//4976KB 1220MS
#include <cstdio>
#include <cctype>
#include <algorithm>
#define lb(x) (x&-x)
//#define gc() getchar()
#define MAXIN 50000
#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
typedef long long LL;
const int N=50005;
int n,m,A[N],Ans[N];
char IN[MAXIN],*SS=IN,*TT=IN;
struct Operation
{
LL K; int l,r,pos;//K!=0:Query [l,r] K; K=0:Add [l,r] pos.
Operation() {}
Operation(LL K,int l,int r,int p):K(K),l(l),r(r),pos(p) {}
}q[N],q1[N],q2[N];
inline int read()
{
int now=0,f=1;register char c=gc();
for(;!isdigit(c);c=='-'&&(f=-1),c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now*f;
}
inline LL readll()
{
LL now=0,f=1;register char c=gc();
for(;!isdigit(c);c=='-'&&(f=-1),c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now*f;
}
namespace T
{
int n;
LL sum1[N],sum2[N];
inline void Modify(int p,int v)
{
for(int i=p; i<=n; i+=lb(i))
sum1[i]+=v, sum2[i]+=p*v;
}
inline void Modify_Range(int l,int r,int v){
Modify(l,v), Modify(r+1,-v);
}
inline LL Query(int p)
{
LL res1=0,res2=0;
for(int i=p; i; i^=lb(i))
res1+=sum1[i], res2+=sum2[i];
return res1*(p+1)-res2;
}
inline LL Query_Range(int l,int r){
return Query(r)-Query(l-1);
}
}
void Solve(int l,int r,int h,int t)
{
if(h>t) return;
if(l==r){
for(int i=h; i<=t; ++i) if(q[i].K) Ans[q[i].pos]=l;
return;
}
bool goon=0;//无询问时直接return。
for(int i=h; i<=t; ++i) if(q[i].K) {goon=1; break;}
if(!goon) return;
int mid=l+r>>1, t1=0, t2=0;
for(int i=h; i<=t; ++i)
if(q[i].K)
{
LL tmp=T::Query_Range(q[i].l,q[i].r);
if(tmp>=q[i].K) q1[t1++]=q[i];
else q[i].K-=tmp, q2[t2++]=q[i];
}
else
{
if(q[i].pos<=mid) T::Modify_Range(q[i].l,q[i].r,1), q1[t1++]=q[i];
else q2[t2++]=q[i];
}
for(int i=0; i<t1; ++i) if(!q1[i].K) T::Modify_Range(q1[i].l,q1[i].r,-1);//q1→_→
for(int i=0; i<t1; ++i) q[h+i]=q1[i];
for(int i=0; i<t2; ++i) q[h+t1+i]=q2[i];
Solve(l,mid,h,h+t1-1), Solve(mid+1,r,h+t1,t);
}
int main()
{
T::n=n=read(), m=read(); int tot=0;
for(int l,r,i=1; i<=m; ++i)
if(read()==1) l=read(), r=read(), q[i]=Operation(0,l,r,n-read());
else l=read(), r=read(), q[i]=Operation(readll(),l,r,++tot);
Solve(-N,N,1,m);
for(int i=1; i<=tot; ++i) printf("%d\n",n-Ans[i]);
return 0;
}
线段树：（常数大的一匹）
//7320kb 6064ms
#include <cstdio>
#include <cctype>
#include <algorithm>
//#define gc() getchar()
#define MAXIN 50000
#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
typedef long long LL;
const int N=50005;
int n,m,A[N],Ans[N];
char IN[MAXIN],*SS=IN,*TT=IN;
struct Operation
{
LL K; int l,r,pos;//K!=0:Query [l,r] K; K=0:Add [l,r] pos.
Operation() {}
Operation(LL K,int l,int r,int p):K(K),l(l),r(r),pos(p) {}
}q[N],q1[N],q2[N];
inline int read()
{
int now=0,f=1;register char c=gc();
for(;!isdigit(c);c=='-'&&(f=-1),c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now*f;
}
inline LL readll()
{
LL now=0,f=1;register char c=gc();
for(;!isdigit(c);c=='-'&&(f=-1),c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now*f;
}
namespace T
{
#define lson rt<<1
#define rson rt<<1|1
#define ToL l,m,rt<<1
#define ToR m+1,r,rt<<1|1
LL sum[N<<2],tag[N<<2];
inline void Update(int rt){
sum[rt]=sum[lson]+sum[rson];
}
inline void PushDown(int rt,LL m)//long long!
{
tag[lson]+=tag[rt], tag[rson]+=tag[rt];
sum[lson]+=tag[rt]*(m-(m>>1)), sum[rson]+=tag[rt]*(m>>1);
tag[rt]=0;
}
void Modify(int l,int r,int rt,int L,int R,int v)
{
if(L<=l && r<=R)
tag[rt]+=v, sum[rt]+=1ll*(r-l+1)*v;
else
{
if(tag[rt]) PushDown(rt,r-l+1);
int m=l+r>>1;
if(L<=m) Modify(ToL,L,R,v);
if(m<R) Modify(ToR,L,R,v);
Update(rt);
}
}
LL Query(int l,int r,int rt,int L,int R)
{
if(L<=l && r<=R) return sum[rt];
if(tag[rt]) PushDown(rt,r-l+1);
int m=l+r>>1;
if(L<=m)
if(m<R) return Query(ToL,L,R)+Query(ToR,L,R);
else return Query(ToL,L,R);
return Query(ToR,L,R);
}
}
void Solve(int l,int r,int h,int t)
{
if(h>t) return;
if(l==r){
for(int i=h; i<=t; ++i) if(q[i].K) Ans[q[i].pos]=l;
return;
}
bool goon=0;//无询问时直接return。
for(int i=h; i<=t; ++i) if(q[i].K) {goon=1; break;}
if(!goon) return;
int mid=l+r>>1, t1=0, t2=0;
for(int i=h; i<=t; ++i)
if(q[i].K)
{
LL tmp=T::Query(1,n,1,q[i].l,q[i].r);
if(tmp>=q[i].K) q1[t1++]=q[i];
else q[i].K-=tmp, q2[t2++]=q[i];
}
else
{
if(q[i].pos<=mid) T::Modify(1,n,1,q[i].l,q[i].r,1), q1[t1++]=q[i];
else q2[t2++]=q[i];
}
for(int i=0; i<t1; ++i) if(!q1[i].K) T::Modify(1,n,1,q1[i].l,q1[i].r,-1);//q1→_→
for(int i=0; i<t1; ++i) q[h+i]=q1[i];
for(int i=0; i<t2; ++i) q[h+t1+i]=q2[i];
Solve(l,mid,h,h+t1-1), Solve(mid+1,r,h+t1,t);
}
int main()
{
n=read(), m=read(); int tot=0;
for(int l,r,i=1; i<=m; ++i)
if(read()==1) l=read(), r=read(), q[i]=Operation(0,l,r,n-read());
else l=read(), r=read(), q[i]=Operation(readll(),l,r,++tot);
Solve(-N,N,1,m);
for(int i=1; i<=tot; ++i) printf("%d\n",n-Ans[i]);
return 0;
}
线段树 标记永久化：WA了，不知道该怎么着。。（应该没问题啊）
//标记永久化替代PushDown可以减小常数，但是怎么不对exm？
#include <cstdio>
#include <cctype>
#include <algorithm>
#define gc() getchar()
//#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
typedef long long LL;
const int N=50005;
int n,m,A[N],Ans[N];
struct Operation
{
LL K; int l,r,pos;//K!=0:Query [l,r] K; K=0:Add [l,r] pos.
Operation() {}
Operation(LL K,int l,int r,int p):K(K),l(l),r(r),pos(p) {}
}q[N],q1[N],q2[N];
inline int read()
{
int now=0,f=1;register char c=gc();
for(;!isdigit(c);c=='-'&&(f=-1),c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now*f;
}
inline LL readll()
{
LL now=0,f=1;register char c=gc();
for(;!isdigit(c);c=='-'&&(f=-1),c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now*f;
}
namespace T
{
#define lson rt<<1
#define rson rt<<1|1
#define ToL l,m,rt<<1
#define ToR m+1,r,rt<<1|1
LL sum[N<<2],all[N<<2];
void Modify(int l,int r,int rt,int L,int R,int v)
{
sum[rt]+=1ll*v*(R-L+1);
if(L<=l && r<=R) all[rt]+=v;
else
{
int m=l+r>>1;
if(L<=m) Modify(ToL,L,m,v);
if(m<R) Modify(ToR,m+1,R,v);
}
}
LL Query(int l,int r,int rt,int L,int R)
{
if(L<=l && r<=R) return sum[rt];
int m=l+r>>1;
if(L<=m)
if(m<R) return 1ll*all[rt]*(R-L+1)+Query(ToL,L,R)+Query(ToR,L,R);
else return 1ll*all[rt]*(R-L+1)+Query(ToL,L,R);
return 1ll*all[rt]*(R-L+1)+Query(ToR,L,R);
}
}
void Solve(int l,int r,int h,int t)
{
if(h>t) return;
if(l==r){
for(int i=h; i<=t; ++i) if(q[i].K) Ans[q[i].pos]=l;
return;
}
bool goon=0;//无询问时直接return。
for(int i=h; i<=t; ++i) if(q[i].K) {goon=1; break;}
if(!goon) return;
int mid=l+r>>1, t1=0, t2=0;
for(int i=h; i<=t; ++i)
if(q[i].K)
{
LL tmp=T::Query(1,n,1,q[i].l,q[i].r);
if(tmp>=q[i].K) q1[t1++]=q[i];
else q[i].K-=tmp, q2[t2++]=q[i];
}
else
{
if(q[i].pos<=mid) T::Modify(1,n,1,q[i].l,q[i].r,1), q1[t1++]=q[i];
else q2[t2++]=q[i];
}
for(int i=0; i<t1; ++i) if(!q1[i].K) T::Modify(1,n,1,q1[i].l,q1[i].r,-1);//q1
for(int i=0; i<t1; ++i) q[h+i]=q1[i];
for(int i=0; i<t2; ++i) q[h+t1+i]=q2[i];
Solve(l,mid,h,h+t1-1), Solve(mid+1,r,h+t1,t);
}
int main()
{
n=read(), m=read(); int tot=0;
for(int l,r,i=1; i<=m; ++i)
if(read()==1) l=read(), r=read(), q[i]=Operation(0,l,r,n-read());
else l=read(), r=read(), q[i]=Operation(readll(),l,r,++tot);
Solve(-N,N,1,m);
for(int i=1; i<=tot; ++i) printf("%d\n",n-Ans[i]);
return 0;
}
BZOJ.3110.[ZJOI2013]K大数查询(整体二分 树状数组/线段树)的更多相关文章
- BZOJ 3110: [Zjoi2013]K大数查询 [整体二分]
有N个位置,M个操作.操作有两种,每次操作如果是1 a b c的形式表示在第a个位置到第b个位置,每个位置加入一个数c如果是2 a b c形式,表示询问从第a个位置到第b个位置,第C大的数是多少. N ...
- BZOJ 3110 [Zjoi2013]K大数查询 ——整体二分
[题目分析] 整体二分显而易见. 自己YY了一下用树状数组区间修改,区间查询的操作. 又因为一个字母调了一下午. 貌似树状数组并不需要清空,可以用一个指针来维护,可以少一个log 懒得写了. [代码] ...
- BZOJ 3110 [ZJOI2013]K大数查询 (整体二分+线段树)
和dynamic rankings这道题的思想一样 只不过是把树状数组换成线段树区间修改,求第$K$大的而不是第$K$小的 这道题还有负数,需要离散 #include <vector> # ...
- BZOJ 3110 [Zjoi2013]K大数查询（整体二分）
3110: [Zjoi2013]K大数查询 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 11654 Solved: 3505[Submit][St ...
- BZOJ 3110: [Zjoi2013]K大数查询 [树套树]
3110: [Zjoi2013]K大数查询 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 6050 Solved: 2007[Submit][Sta ...
- bzoj 3110: [Zjoi2013]K大数查询 树状数组套线段树
3110: [Zjoi2013]K大数查询 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 1384 Solved: 629[Submit][Stat ...
- 树套树专题——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大数查询( 树状数组套主席树 )
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][ ...
随机推荐
- python——脚本和print
脚本和print 1.脚本文件 <Python 基础教程>(第二版)中 P118页,原操作为下: 1 _metaclass_ = type 2 3 class Person: 4 def ...
- python小工具之读取host文件
# -*- coding: utf-8 -*- # @Time : 2018/9/12 21:09 # @Author : cxa # @File : readhostfile.py # @Softw ...
- redis数据淘汰机制
volatile-lru:从已设置过期时间的数据集(server.db[i].expires)中挑选最近最少使用 的数据淘汰 volatile-ttl:从已设置过期时间的数据集(server.db[i ...
- 在Mac上搭建ReactNative开发环境
1.安装Homebrew, Mac系统的包管理器,用于安装NodeJS和一些其他必需的工具软件. /usr/bin/ruby -e "$(curl -fsSL https://raw.g ...
- 密码记录工具keepass保存密码
https://www.cnblogs.com/wicub/p/5753005.html
- 简单计算器的C实现-函数指针，main函数传参
/** 程序功能:简单计算器,实现加减乘除平方* 作者版本日期:2015.11.08 zhouhb OK* 源代码:李明 <新概念C语言培训>第33集 C语言Shell命令解释器的实现* ...
- No.1 selenium学习之路之浏览器操作
selenium基础,首先就是浏览器的相关操作 下面描述几种浏览器的常用操作 1.打开浏览器 webdriver后面添加想要打开的浏览器 Ie或者Chrome 2.打开指定页面(百度) 3.休眠时间 ...
- Zookeeper安装以及配置说明（三）
Zookeeper的安装和配置非常的简单,既可以配置成单机模式,也可以配置成集群模式.如下图所示: 下面将分别进行介绍: 单机模式 下载最新稳定版本zookeeper的安装包之后(看第一篇博文), 解 ...
- SQL Server日期计算
通常,你需要获得当前日期和计算一些其他的日期,例如,你的程序可能需要判断一个月的第一天或者最后一天.你们大部分人大概都知道怎样把日期进行分割(年.月.日等),然后仅仅用分割出来的年.月.日等放在几个函 ...
- Oracle约束
1.非空约束 DROP TABLE member PURGE; CREATE TABLE member( mid NUMBER, name ) NOT NULL ); 2.唯一约束 DROP TABL ...