题目链接 BZOJ

洛谷

整体二分求的是第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多。。

树状数组：

  1. //4976KB	1220MS
  2. #include <cstdio>
  3. #include <cctype>
  4. #include <algorithm>
  5. #define lb(x) (x&-x)
  6. //#define gc() getchar()
  7. #define MAXIN 50000
  8. #define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
  9. typedef long long LL;
  10. const int N=50005;
  12. int n,m,A[N],Ans[N];
  13. char IN[MAXIN],*SS=IN,*TT=IN;
  14. struct Operation
  15. {
  16. 	LL K; int l,r,pos;//K!=0:Query [l,r] K;	K=0:Add [l,r] pos.
  17. 	Operation() {}
  18. 	Operation(LL K,int l,int r,int p):K(K),l(l),r(r),pos(p) {}
  19. }q[N],q1[N],q2[N];
  21. inline int read()
  22. {
  23. 	int now=0,f=1;register char c=gc();
  24. 	for(;!isdigit(c);c=='-'&&(f=-1),c=gc());
  25. 	for(;isdigit(c);now=now*10+c-'0',c=gc());
  26. 	return now*f;
  27. }
  28. inline LL readll()
  29. {
  30. 	LL now=0,f=1;register char c=gc();
  31. 	for(;!isdigit(c);c=='-'&&(f=-1),c=gc());
  32. 	for(;isdigit(c);now=now*10+c-'0',c=gc());
  33. 	return now*f;
  34. }
  35. namespace T
  36. {
  37. 	int n;
  38. 	LL sum1[N],sum2[N];
  40. 	inline void Modify(int p,int v)
  41. 	{
  42. 		for(int i=p; i<=n; i+=lb(i))
  43. 			sum1[i]+=v, sum2[i]+=p*v;
  44. 	}
  45. 	inline void Modify_Range(int l,int r,int v){
  46. 		Modify(l,v), Modify(r+1,-v);
  47. 	}
  48. 	inline LL Query(int p)
  49. 	{
  50. 		LL res1=0,res2=0;
  51. 		for(int i=p; i; i^=lb(i))
  52. 			res1+=sum1[i], res2+=sum2[i];
  53. 		return res1*(p+1)-res2;
  54. 	}
  55. 	inline LL Query_Range(int l,int r){
  56. 		return Query(r)-Query(l-1);
  57. 	}
  58. }
  59. void Solve(int l,int r,int h,int t)
  60. {
  61. 	if(h>t) return;
  62. 	if(l==r){
  63. 		for(int i=h; i<=t; ++i) if(q[i].K) Ans[q[i].pos]=l;
  64. 		return;
  65. 	}
  66. 	bool goon=0;//无询问时直接return。
  67. 	for(int i=h; i<=t; ++i) if(q[i].K) {goon=1; break;}
  68. 	if(!goon) return;
  70. 	int mid=l+r>>1, t1=0, t2=0;
  71. 	for(int i=h; i<=t; ++i)
  72. 		if(q[i].K)
  73. 		{
  74. 			LL tmp=T::Query_Range(q[i].l,q[i].r);
  75. 			if(tmp>=q[i].K) q1[t1++]=q[i];
  76. 			else q[i].K-=tmp, q2[t2++]=q[i];
  77. 		}
  78. 		else
  79. 		{
  80. 			if(q[i].pos<=mid) T::Modify_Range(q[i].l,q[i].r,1), q1[t1++]=q[i];
  81. 			else q2[t2++]=q[i];
  82. 		}
  83. 	for(int i=0; i<t1; ++i) if(!q1[i].K) T::Modify_Range(q1[i].l,q1[i].r,-1);//q1→_→
  84. 	for(int i=0; i<t1; ++i) q[h+i]=q1[i];
  85. 	for(int i=0; i<t2; ++i) q[h+t1+i]=q2[i];
  86. 	Solve(l,mid,h,h+t1-1), Solve(mid+1,r,h+t1,t);
  87. }
  89. int main()
  90. {
  91. 	T::n=n=read(), m=read(); int tot=0;
  92. 	for(int l,r,i=1; i<=m; ++i)
  93. 		if(read()==1) l=read(), r=read(), q[i]=Operation(0,l,r,n-read());
  94. 		else l=read(), r=read(), q[i]=Operation(readll(),l,r,++tot);
  95. 	Solve(-N,N,1,m);
  96. 	for(int i=1; i<=tot; ++i) printf("%d\n",n-Ans[i]);
  98. 	return 0;
  99. }

线段树：（常数大的一匹）

  1. //7320kb	6064ms
  2. #include <cstdio>
  3. #include <cctype>
  4. #include <algorithm>
  5. //#define gc() getchar()
  6. #define MAXIN 50000
  7. #define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
  8. typedef long long LL;
  9. const int N=50005;
  11. int n,m,A[N],Ans[N];
  12. char IN[MAXIN],*SS=IN,*TT=IN;
  13. struct Operation
  14. {
  15. 	LL K; int l,r,pos;//K!=0:Query [l,r] K;	K=0:Add [l,r] pos.
  16. 	Operation() {}
  17. 	Operation(LL K,int l,int r,int p):K(K),l(l),r(r),pos(p) {}
  18. }q[N],q1[N],q2[N];
  20. inline int read()
  21. {
  22. 	int now=0,f=1;register char c=gc();
  23. 	for(;!isdigit(c);c=='-'&&(f=-1),c=gc());
  24. 	for(;isdigit(c);now=now*10+c-'0',c=gc());
  25. 	return now*f;
  26. }
  27. inline LL readll()
  28. {
  29. 	LL now=0,f=1;register char c=gc();
  30. 	for(;!isdigit(c);c=='-'&&(f=-1),c=gc());
  31. 	for(;isdigit(c);now=now*10+c-'0',c=gc());
  32. 	return now*f;
  33. }
  34. namespace T
  35. {
  36. 	#define lson rt<<1
  37. 	#define rson rt<<1|1
  38. 	#define ToL l,m,rt<<1
  39. 	#define ToR m+1,r,rt<<1|1
  41. 	LL sum[N<<2],tag[N<<2];
  43. 	inline void Update(int rt){
  44. 		sum[rt]=sum[lson]+sum[rson];
  45. 	}
  46. 	inline void PushDown(int rt,LL m)//long long!
  47. 	{
  48. 		tag[lson]+=tag[rt], tag[rson]+=tag[rt];
  49. 		sum[lson]+=tag[rt]*(m-(m>>1)), sum[rson]+=tag[rt]*(m>>1);
  50. 		tag[rt]=0;
  51. 	}
  52. 	void Modify(int l,int r,int rt,int L,int R,int v)
  53. 	{
  54. 		if(L<=l && r<=R)
  55. 			tag[rt]+=v, sum[rt]+=1ll*(r-l+1)*v;
  56. 		else
  57. 		{
  58. 			if(tag[rt]) PushDown(rt,r-l+1);
  59. 			int m=l+r>>1;
  60. 			if(L<=m) Modify(ToL,L,R,v);
  61. 			if(m<R) Modify(ToR,L,R,v);
  62. 			Update(rt);
  63. 		}
  64. 	}
  65. 	LL Query(int l,int r,int rt,int L,int R)
  66. 	{
  67. 		if(L<=l && r<=R) return sum[rt];
  68. 		if(tag[rt]) PushDown(rt,r-l+1);
  69. 		int m=l+r>>1;
  70. 		if(L<=m)
  71. 			if(m<R) return Query(ToL,L,R)+Query(ToR,L,R);
  72. 			else return Query(ToL,L,R);
  73. 		return Query(ToR,L,R);
  74. 	}
  75. }
  76. void Solve(int l,int r,int h,int t)
  77. {
  78. 	if(h>t) return;
  79. 	if(l==r){
  80. 		for(int i=h; i<=t; ++i) if(q[i].K) Ans[q[i].pos]=l;
  81. 		return;
  82. 	}
  83. 	bool goon=0;//无询问时直接return。
  84. 	for(int i=h; i<=t; ++i) if(q[i].K) {goon=1; break;}
  85. 	if(!goon) return;
  87. 	int mid=l+r>>1, t1=0, t2=0;
  88. 	for(int i=h; i<=t; ++i)
  89. 		if(q[i].K)
  90. 		{
  91. 			LL tmp=T::Query(1,n,1,q[i].l,q[i].r);
  92. 			if(tmp>=q[i].K) q1[t1++]=q[i];
  93. 			else q[i].K-=tmp, q2[t2++]=q[i];
  94. 		}
  95. 		else
  96. 		{
  97. 			if(q[i].pos<=mid) T::Modify(1,n,1,q[i].l,q[i].r,1), q1[t1++]=q[i];
  98. 			else q2[t2++]=q[i];
  99. 		}
  100. 	for(int i=0; i<t1; ++i) if(!q1[i].K) T::Modify(1,n,1,q1[i].l,q1[i].r,-1);//q1→_→
  101. 	for(int i=0; i<t1; ++i) q[h+i]=q1[i];
  102. 	for(int i=0; i<t2; ++i) q[h+t1+i]=q2[i];
  103. 	Solve(l,mid,h,h+t1-1), Solve(mid+1,r,h+t1,t);
  104. }
  106. int main()
  107. {
  108. 	n=read(), m=read(); int tot=0;
  109. 	for(int l,r,i=1; i<=m; ++i)
  110. 		if(read()==1) l=read(), r=read(), q[i]=Operation(0,l,r,n-read());
  111. 		else l=read(), r=read(), q[i]=Operation(readll(),l,r,++tot);
  112. 	Solve(-N,N,1,m);
  113. 	for(int i=1; i<=tot; ++i) printf("%d\n",n-Ans[i]);
  115. 	return 0;
  116. }

线段树 标记永久化：WA了，不知道该怎么着。。（应该没问题啊）

  1. //标记永久化替代PushDown可以减小常数，但是怎么不对exm？
  2. #include <cstdio>
  3. #include <cctype>
  4. #include <algorithm>
  5. #define gc() getchar()
  6. //#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
  7. typedef long long LL;
  8. const int N=50005;
  10. int n,m,A[N],Ans[N];
  11. struct Operation
  12. {
  13. 	LL K; int l,r,pos;//K!=0:Query [l,r] K;	K=0:Add [l,r] pos.
  14. 	Operation() {}
  15. 	Operation(LL K,int l,int r,int p):K(K),l(l),r(r),pos(p) {}
  16. }q[N],q1[N],q2[N];
  18. inline int read()
  19. {
  20. 	int now=0,f=1;register char c=gc();
  21. 	for(;!isdigit(c);c=='-'&&(f=-1),c=gc());
  22. 	for(;isdigit(c);now=now*10+c-'0',c=gc());
  23. 	return now*f;
  24. }
  25. inline LL readll()
  26. {
  27. 	LL now=0,f=1;register char c=gc();
  28. 	for(;!isdigit(c);c=='-'&&(f=-1),c=gc());
  29. 	for(;isdigit(c);now=now*10+c-'0',c=gc());
  30. 	return now*f;
  31. }
  32. namespace T
  33. {
  34. 	#define lson rt<<1
  35. 	#define rson rt<<1|1
  36. 	#define ToL l,m,rt<<1
  37. 	#define ToR m+1,r,rt<<1|1
  39. 	LL sum[N<<2],all[N<<2];
  41. 	void Modify(int l,int r,int rt,int L,int R,int v)
  42. 	{
  43. 		sum[rt]+=1ll*v*(R-L+1);
  44. 		if(L<=l && r<=R) all[rt]+=v;
  45. 		else
  46. 		{
  47. 			int m=l+r>>1;
  48. 			if(L<=m) Modify(ToL,L,m,v);
  49. 			if(m<R) Modify(ToR,m+1,R,v);
  50. 		}
  51. 	}
  52. 	LL Query(int l,int r,int rt,int L,int R)
  53. 	{
  54. 		if(L<=l && r<=R) return sum[rt];
  55. 		int m=l+r>>1;
  56. 		if(L<=m)
  57. 			if(m<R) return 1ll*all[rt]*(R-L+1)+Query(ToL,L,R)+Query(ToR,L,R);
  58. 			else return 1ll*all[rt]*(R-L+1)+Query(ToL,L,R);
  59. 		return 1ll*all[rt]*(R-L+1)+Query(ToR,L,R);
  60. 	}
  61. }
  62. void Solve(int l,int r,int h,int t)
  63. {
  64. 	if(h>t) return;
  65. 	if(l==r){
  66. 		for(int i=h; i<=t; ++i) if(q[i].K) Ans[q[i].pos]=l;
  67. 		return;
  68. 	}
  69. 	bool goon=0;//无询问时直接return。
  70. 	for(int i=h; i<=t; ++i) if(q[i].K) {goon=1; break;}
  71. 	if(!goon) return;
  73. 	int mid=l+r>>1, t1=0, t2=0;
  74. 	for(int i=h; i<=t; ++i)
  75. 		if(q[i].K)
  76. 		{
  77. 			LL tmp=T::Query(1,n,1,q[i].l,q[i].r);
  78. 			if(tmp>=q[i].K) q1[t1++]=q[i];
  79. 			else q[i].K-=tmp, q2[t2++]=q[i];
  80. 		}
  81. 		else
  82. 		{
  83. 			if(q[i].pos<=mid) T::Modify(1,n,1,q[i].l,q[i].r,1), q1[t1++]=q[i];
  84. 			else q2[t2++]=q[i];
  85. 		}
  86. 	for(int i=0; i<t1; ++i) if(!q1[i].K) T::Modify(1,n,1,q1[i].l,q1[i].r,-1);//q1
  87. 	for(int i=0; i<t1; ++i) q[h+i]=q1[i];
  88. 	for(int i=0; i<t2; ++i) q[h+t1+i]=q2[i];
  89. 	Solve(l,mid,h,h+t1-1), Solve(mid+1,r,h+t1,t);
  90. }
  92. int main()
  93. {
  94. 	n=read(), m=read(); int tot=0;
  95. 	for(int l,r,i=1; i<=m; ++i)
  96. 		if(read()==1) l=read(), r=read(), q[i]=Operation(0,l,r,n-read());
  97. 		else l=read(), r=read(), q[i]=Operation(readll(),l,r,++tot);
  98. 	Solve(-N,N,1,m);
  99. 	for(int i=1; i<=tot; ++i) printf("%d\n",n-Ans[i]);
  101. 	return 0;
  102. }

BZOJ.3110.[ZJOI2013]K大数查询(整体二分 树状数组/线段树)的更多相关文章

  1. BZOJ 3110: [Zjoi2013]K大数查询 [整体二分]

    有N个位置,M个操作.操作有两种,每次操作如果是1 a b c的形式表示在第a个位置到第b个位置,每个位置加入一个数c如果是2 a b c形式,表示询问从第a个位置到第b个位置,第C大的数是多少. N ...

  2. BZOJ 3110 [Zjoi2013]K大数查询 ——整体二分

    [题目分析] 整体二分显而易见. 自己YY了一下用树状数组区间修改,区间查询的操作. 又因为一个字母调了一下午. 貌似树状数组并不需要清空,可以用一个指针来维护,可以少一个log 懒得写了. [代码] ...

  3. BZOJ 3110 [ZJOI2013]K大数查询 (整体二分+线段树)

    和dynamic rankings这道题的思想一样 只不过是把树状数组换成线段树区间修改,求第$K$大的而不是第$K$小的 这道题还有负数,需要离散 #include <vector> # ...

  4. BZOJ 3110 [Zjoi2013]K大数查询（整体二分）

    3110: [Zjoi2013]K大数查询 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 11654  Solved: 3505[Submit][St ...

  5. BZOJ 3110: [Zjoi2013]K大数查询 [树套树]

    3110: [Zjoi2013]K大数查询 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 6050  Solved: 2007[Submit][Sta ...

  6. bzoj 3110: [Zjoi2013]K大数查询 树状数组套线段树

    3110: [Zjoi2013]K大数查询 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 1384  Solved: 629[Submit][Stat ...

  7. 树套树专题——bzoj 3110: [Zjoi2013] K大数查询 &amp; 3236 [Ahoi2013] 作业 题解

    [原题1] 3110: [Zjoi2013]K大数查询 Time Limit: 20 Sec  Memory Limit: 512 MB Submit: 978  Solved: 476 Descri ...

  8. BZOJ 3110: [Zjoi2013]K大数查询( 树状数组套主席树 )

    BIT+(可持久化)权值线段树, 用到了BIT的差分技巧. 时间复杂度O(Nlog^2(N)) ---------------------------------------------------- ...

  9. BZOJ 3110([Zjoi2013]K大数查询-区间第k大[段修改,在线]-树状数组套函数式线段树)

    3110: [Zjoi2013]K大数查询 Time Limit: 20 Sec   Memory Limit: 512 MB Submit: 418   Solved: 235 [ Submit][ ...

随机推荐

  1. webSQL 增删改查

    webSQL 增删改查  转载于:https://www.cnblogs.com/liuhao-web/p/7866032.html /** *数据库操作辅助类,定义对象.数据操作方法都在这里定义 * ...

  2. imperva 更改web界面的密码

    通过SSH作为用户根登录到MX(或通过另一个用户并提升) 运行命令“su oracle” //首先切换到oracle用户 sqlplus secure/(密码)   /用此命令登录到数据库     s ...

  3. Linux学习笔记-Linux系统简介

    Linux学习笔记-Linux系统简介 UNIX与Linux发展史 UNIX是父亲,Linux是儿子. UNIX发行版本 操作系统 公司 硬件平台 AIX IBM PowerPC HP-UX HP P ...

  4. Flask：静态文件&模板（0.1）

    Windows 10家庭中文版,Python 3.6.4,Flask 1.0.2 前面看了Flask的Quickstart文档,可是,一直没有练习里面的内容,这不,刚刚练习完毕,来写篇博文记录一下! ...

  5. 分别使用docx4j,jacob将文字与图片插入word中书签位置

    项目中需要将一段文字,与人员的签名(图片)插入到上传的word中,上网查询了一下,有许多种方式可以向word中插入文字,发现docx4j与jacob都为比较常见的解决方案,于是就先使用的docx4j进 ...

  6. numpy数学计算

    1.求范数 np.linalg.norm norm(x, ord=None, axis=None, keepdims=False)  范数理论的一个小推论告诉我们:ℓ1≥ℓ2≥ℓ∞

  7. 2016-2017-2 20155309南皓芯《java程序设计》第七周学习总结

    教材学习内容总结 Lambda 一种匿名方法 表达式构成 括号以及括号里用逗号分隔的参数列表 仅有一个参数的可以省略括号 ->符号 花括号以及花括号里的语句 仅有一条语句时可以省略花括号,并且这 ...

  8. MP3 Fuzz学习

    这篇文章主要是学习一波MP3格式fuzz的知识.目录如下 0x0.MP3格式的构成 0x0.MP3格式的构成 MP3是一种通俗叫法,学名叫MPEG1 Layer-3.MP3是三段式的结构,依次由ID3 ...

  9. win7下scheme环境配置

    运行lisp方言--scheme,在windows下,用eclipse. 1.eclipse装好 2.eclipse安装插件scheme48 development took, http://www. ...

  10. Spark介绍及安装部署

    一.Spark介绍 1.1 Apache Spark Apache Spark是一个围绕速度.易用性和复杂分析构建的大数据处理框架(没有数据存储).最初在2009年由加州大学伯克利分校的AMPLab开 ...