【fhqtreap】航空管制

最新推荐文章于 2024-01-05 23:05:51 发布
最新推荐文章于 2024-01-05 23:05:51 发布
阅读量1.9k 收藏
点赞数
分类专栏： poj 数据结构
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/huyuncong/article/details/22963863
版权

飞机因为航空管制晚点了一个多小时，于是正好学习一下fhq介绍过的treap...

poj2352

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
using namespace std;
int l[200000],r[200000],s[2000000],w[200000],x[200000],y[200000];
int ss,root,n,ans[200000];
bool random(double p)
{
	return (double)rand()/RAND_MAX<p;
}
void print(int x)
{
	if (!x) return ;
	cout<<x<<' '<<l[x]<<' '<<r[x]<<endl;
	print(l[x]),print(r[x]);
}
int ori(int x)
{
	++ss,l[ss]=r[ss]=0;
	w[ss]=x,s[ss]=1;
	return ss;
}
void updata(int x)
{
	s[x]=s[l[x]]+1+s[r[x]];
}
int merge(int x,int y)
{
	if (!x) return y;
	if (!y) return x;
	if (random((double)s[x]/(s[x]+s[y]))) {
		r[x]=merge(r[x],y),updata(x);
		return x;
	}
	l[y]=merge(x,l[y]),updata(y);
	return y;
}
void split(int root,int k,int &a,int &b)
{
	if (!root) {
		a=b=0;
		return ;
	}
	if (w[root]>k) {								 
		int p,q;
		split(l[root],k,p,q);
		l[root]=0,updata(root);
		a=p,b=merge(q,root);
		return ;
	}
	int p,q;													//保证根节点被选 
	split(r[root],k,p,q);
	r[root]=0,updata(root);
	a=merge(root,p),b=q;
	return ;
}
int ask(int x)
{
	int a,b;
	split(root,x,a,b);
	int ans=s[a];
	root=merge(a,b);
	return ans;
}
int main()
{
	freopen("input.txt","r",stdin);
	scanf("%d",&n);
	for (int i=1;i<=n;i++) scanf("%d%d",&x[i],&y[i]);
	root=ss=0;
	for (int i=1;i<=n;i++) {
		int sum=ask(x[i]);
		ans[sum]++;
		int a,b;
		split(root,x[i],a,b);
		root=ori(x[i]);
		root=merge(a,root),root=merge(root,b);
	}
	for (int i=0;i<n;i++) printf("%d\n",ans[i]);
	return 0;
}

poj3580

上飞机前没调出来...

用叉姐的rand函数tle了，原因是值域只有32767,但是即便扩大范围也用了1900+ms,而黄哥的程序不但各种指针，stl，rand函数值域也只有32767还比我要快，难道真的是数组寻址太慢了不成...不过数据范围应该没那么大

值得一提的是我以前同样是数组实现的splay也比这个treap快的多，虽然说用这个treap处理某些操作非常方便...

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <ctime>
const int oo=1073741819;	
const int Rand_Max=2147483647;
using namespace std;
int l[250000],r[250000],w[250000],s[250000],Min[250000];
int Ad[250000],Sw[250000];
int root,ss,n,m;
void print(int x)
{
	if (!x) return ;
	cout<<x<<' '<<l[x]<<' '<<r[x]<<endl;
	print(l[x]);
	print(r[x]);
}
inline int ran()
{
	static int x = 1;
	x += (x << 2) + 1;
	return x & 2147483647;
}

bool random(double p)
{
	return (double) (ran())/Rand_Max < p;
}
int ori(int x)
{
	++ss,l[ss]=r[ss]=0;
	w[ss]=x,s[ss]=1,Min[ss]=x;
	Ad[ss]=0,Sw[ss]=0;
	return ss;
}
void add(int x,int W)
{
	if (!x) return ;
	w[x]+=W,Ad[x]+=W,Min[x]+=W;
}
void swap(int x)
{
	if (!x) return ;
	int e;
	e=l[x],l[x]=r[x],r[x]=e;
	Sw[x]^=1;
}
void updata(int x)
{
	if (!x) return ;
	s[x]=s[l[x]]+s[r[x]]+1;
	Min[x]=min(Min[l[x]],min(Min[r[x]],w[x]));
}
void pushdown(int x)
{
	if (!x) return ;
	if (Ad[x]) {
		add(l[x],Ad[x]),add(r[x],Ad[x]);
		Ad[x]=0;
	}
	if (Sw[x]) {
		swap(l[x]),swap(r[x]);
		Sw[x]=0;
	}
	updata(x);
}
int merge(int x,int y)
{
	if (!x) return y;
	if (!y) return x;
	//if (random((double)s[x]/(s[x]+s[y]))) {
	//if (w[x]<w[y]) {
	if (ran()%(s[x]+s[y])<s[x]) {
		pushdown(x);
		r[x]=merge(r[x],y);
		updata(x);
		return x;
	}
	pushdown(y);
	l[y]=merge(x,l[y]);
	updata(y);
	return y;
}
void split(int x,int k,int &a,int &b)
{
	if (!x) {
		a=b=0;
		return ;
	}
	pushdown(x);
	if (s[l[x]]>=k) {
		int p,q;
		split(l[x],k,p,q);
		l[x]=0,updata(x);
		a=p,b=merge(q,x);
		return ;
	}
	int p,q;
	split(r[x],k-(s[l[x]]+1),p,q);
	r[x]=0,updata(x);
	a=merge(x,p),b=q;
	return ;
}
int main()
{
	freopen("input.txt","r",stdin);
	freopen("output.txt","w",stdout);
	scanf("%d",&n);
	root=ss=0,Min[0]=oo,s[0]=0,w[0]=oo;
	for (int i=1;i<=n;i++) {
		int x;
		scanf("%d",&x);
		int y=ori(x);
		root=merge(root,y);
	}
	scanf("%d",&m);
	char ch[200];
	//print(root);
	for (int i=1;i<=m;i++) {
		scanf("%s",ch+1);
		if ('A'==ch[1]) {
			int x,y,z;
			scanf("%d%d%d",&x,&y,&z);
			int a,b,c,d;
			split(root,x-1,a,b);
			split(b,y-x+1,c,d);
			add(c,z);
			root=merge(a,merge(c,d));
		}
		else if ('R'==ch[1] && 'E'==ch[4]) {
			int x,y;
			scanf("%d%d",&x,&y);
			int a,b,c,d;
			split(root,x-1,a,b);
			split(b,y-x+1,c,d);
			swap(c);
			root=merge(a,merge(c,d));
		}
		else if ('R'==ch[1] && 'O'==ch[4]) {
			int x,y,z;
			scanf("%d%d%d",&x,&y,&z);
			//if ((z>(y-x+1))) continue;
			z%=(y-x+1);
			int a,b,c,d,e,f;
			split(root,x-1,a,b);
			split(b,y-x+1,c,d);
			split(c,(y-x+1)-z,e,f);
			root=merge(a,merge(f,merge(e,d)));
		}
		else if ('I'==ch[1]) {
			int x,y;
			scanf("%d%d",&x,&y);
			int p=ori(y);
			int a,b;
			split(root,x,a,b);
			root=merge(a,merge(p,b));
		}
		else if ('D'==ch[1]) {
			int x;
			scanf("%d",&x);
			int a,b,c,d;
			split(root,x-1,a,b);
			split(b,1,c,d);
			root=merge(a,d);
		}
		else if ('M'==ch[1]) {
			int x,y;
			scanf("%d%d",&x,&y);
			int a,b,c,d;
			split(root,x-1,a,b);
			split(b,y-x+1,c,d);
			int ans=Min[c];
			root=merge(a,merge(c,d));
			printf("%d\n",ans);
		}
	}
	cout<<clock()<<endl;
	return 0;
}


update 2014.4.14

去掉了一个无用的updata，时间减少了100ms，又去掉了split的两个merge，时间瞬间减少了400+ms

果然细节对常数的影响巨大，不过split的两个merge确实是没有必要的，但是叉姐的板子里面有...

update 2014.5.11 

其实我也不知道去掉split的merge会不会导致树高不是log= =

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <ctime>
const int oo=1073741819;	
const int Rand_Max=1073741819;
using namespace std;
int l[250000],r[250000],w[250000],s[250000],Min[250000];
int Ad[250000],Sw[250000];
int root,ss,n,m;
void print(int x)
{
	if (!x) return ;
	cout<<x<<' '<<l[x]<<' '<<r[x]<<endl;
	print(l[x]);
	print(r[x]);
}
inline int ran()
{
	static int x = 1;
	x += (x << 2) + 1;
	return x & 1073741819;
}

bool random(double p)
{
	return (double) (ran())/Rand_Max < p;
}
int ori(int x)
{
	++ss,l[ss]=r[ss]=0;
	w[ss]=x,s[ss]=1,Min[ss]=x;
	Ad[ss]=0,Sw[ss]=0;
	return ss;
}
void add(int x,int W)
{
	if (!x) return ;
	w[x]+=W,Ad[x]+=W,Min[x]+=W;
}
void swap(int x)
{
	if (!x) return ;
	int e;
	e=l[x],l[x]=r[x],r[x]=e;
	Sw[x]^=1;
}
void updata(int x)
{
	if (!x) return ;
	s[x]=s[l[x]]+s[r[x]]+1;
	Min[x]=min(Min[l[x]],min(Min[r[x]],w[x]));
}
void pushdown(int x)
{
	if (!x) return ;
	if (Ad[x]) {
		add(l[x],Ad[x]),add(r[x],Ad[x]);
		Ad[x]=0;
	}
	if (Sw[x]) {
		swap(l[x]),swap(r[x]);
		Sw[x]=0;
	}
	//updata(x);
}
int merge(int x,int y)
{
	if (!x) return y;
	if (!y) return x;
	//if (random((double)s[x]/(s[x]+s[y]))) {
	//if (w[x]<w[y]) {
	if (ran()%(s[x]+s[y])<s[x]) {
		pushdown(x);
		r[x]=merge(r[x],y);
		updata(x);
		return x;
	}
	pushdown(y);
	l[y]=merge(x,l[y]);
	updata(y);
	return y;
}
void split(int x,int k,int &a,int &b)
{
	if (!x) {
		a=b=0;
		return ;
	}
	pushdown(x);
	if (s[l[x]]>=k) {
		int p,q;
		split(l[x],k,p,q);
		l[x]=q,updata(x);
		a=p,b=x;
		return ;
	}
	int p,q;
	split(r[x],k-(s[l[x]]+1),p,q);
	r[x]=p,updata(x);
	a=x,b=q;
	return ;
}
int main()
{
	freopen("input.txt","r",stdin);
	freopen("output.txt","w",stdout);
	scanf("%d",&n);
	root=ss=0,Min[0]=oo,s[0]=0,w[0]=oo;
	for (int i=1;i<=n;i++) {
		int x;
		scanf("%d",&x);
		int y=ori(x);
		root=merge(root,y);
	}
	scanf("%d",&m);
	char ch[200];
	//print(root);
	for (int i=1;i<=m;i++) {
		scanf("%s",ch+1);
		if ('A'==ch[1]) {
			int x,y,z;
			scanf("%d%d%d",&x,&y,&z);
			int a,b,c,d;
			split(root,x-1,a,b);
			split(b,y-x+1,c,d);
			add(c,z);
			root=merge(a,merge(c,d));
		}
		else if ('R'==ch[1] && 'E'==ch[4]) {
			int x,y;
			scanf("%d%d",&x,&y);
			int a,b,c,d;
			split(root,x-1,a,b);
			split(b,y-x+1,c,d);
			swap(c);
			root=merge(a,merge(c,d));
		}
		else if ('R'==ch[1] && 'O'==ch[4]) {
			int x,y,z;
			scanf("%d%d%d",&x,&y,&z);
			//if ((z>(y-x+1))) continue;
			z%=(y-x+1);
			int a,b,c,d,e,f;
			split(root,x-1,a,b);
			split(b,y-x+1,c,d);
			split(c,(y-x+1)-z,e,f);
			root=merge(a,merge(f,merge(e,d)));
		}
		else if ('I'==ch[1]) {
			int x,y;
			scanf("%d%d",&x,&y);
			int p=ori(y);
			int a,b;
			split(root,x,a,b);
			root=merge(a,merge(p,b));
		}
		else if ('D'==ch[1]) {
			int x;
			scanf("%d",&x);
			int a,b,c,d;
			split(root,x-1,a,b);
			split(b,1,c,d);
			root=merge(a,d);
		}
		else if ('M'==ch[1]) {
			int x,y;
			scanf("%d%d",&x,&y);
			int a,b,c,d;
			split(root,x-1,a,b);
			split(b,y-x+1,c,d);
			int ans=Min[c];
			root=merge(a,merge(c,d));
			printf("%d\n",ans);
		}
	}
	//cout<<clock()<<endl;
	return 0;
}


huyuncong
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
【数据结构】FHQ Treap详解
聚娃科技开发者博客
01-27 1262
FHQ Treap是什么？ FHQ Treap，又名无旋Treap，是一种不需要旋转的平衡树，是范浩强基于Treap发明的。FHQ Treap具有代码短，易理解，速度快的优点。（当然跟红黑树比一下就是……）至少它在OI中算是很优秀的数据结构了。 前置知识： C++ 二叉搜索树的基本性质，下面会讲 二叉堆 二叉搜索树的基本性质 很简单，就这几个。 在二叉搜索树中，每个结点都满足左...
FHQ-Treap(非旋treap/平衡树)——从入门到入坟
yyh_getAC的博客
07-11 4865
平衡树有很多种，其中fhq-treap 可以说是最强的一种平衡树之一，它可以维护值域，也可以维护下标，还能维护区间修改，更难能可贵的是，它可以完成splay都完成不了的可持久化操作。其唯一弱于splay的就是在LCT上，splay维护LCT的复杂度是O(nlog(n)),而fhq-treap的复杂度为O(nlog^2(n))，稍微大了一点，但是比splay好写多了。............
FHQ Treap学习记录（详解）
CE_Automata的博客
03-01 700
FHQ Treap，其中 FHQ 指范浩强神犇，是依赖于分裂合并操作实现的 Treap，这种操作方式使得它天生支持维护序列、可持久化等特性……
FHQ Treap 详解
AKAPinkman的博客
06-21 5288
0.简单介绍 FHQ Treap,以下简写为fhq,是一种treap(树堆)的变体,功能比treap强大,代码比splay好写,易于理解,常数稍大. fhq不需要通过一般平衡树的左右旋转来保持平衡,而是通过分裂(split)和合并(merge)来实现操作. 本文力求简明易懂,在看本文前最好先了解过一种平衡树的基本写法,否则代码部分有可能看不懂. 1. 基本操作 结构 以结构体作为树的每一个节点,存储左子树和右子树的位置,权值,堆权值,和子树大小. struct Node { int l, r;
fhq Treap
weixin_51606521的博客
02-10 1697
普通Treap，既Tree加Heap。它让平衡树上的每一个结点存放两个信息：值和一个随机的索引。其中值满足二叉搜索树的性质，索引满足堆的性质，结合二叉搜索树和二叉堆的性质来使树平衡。这也是Treap的性质。 Treap用二叉堆来维护随机索引，其实就是相当于把插入次序随机化。插入一数值后你必然要让索引去满足二叉堆的特性，但又因为索引是随机的，那就会导致插入的数后不知道搞到哪里去了，相当于插入次序随机了。 （随机大法好，吊打一切毒瘤数据） fhq Treap的奇怪操作并且是核心操作只有两个：分裂和合并
2020年9月最新民用航空管制频率表.pdf
01-22
2020年9月最新民用航空管制频率表.pdf
2018年上半年最新航空管制频率表_2018年上半年最新航空管制频率表_
10-04
航空管制频率表是航空交通管理中的关键参考资料，对于飞行员、空中交通管制员以及航空爱好者来说都具有重要意义。2018年上半年的最新航空管制频率表详细列出了各个机场、航路和飞行情报区(FIR)的通信频率，确保了...
最新航空管制频.doc
06-06
"航空管制频率" 本文档主要介绍了航空管制频率的相关知识点，涵盖了航空管制的各个阶段，从飞机起降过程到航路航线的分类。 航空管制过程 航空管制过程可以分为五个阶段：地面控制（Ground Control）、塔台控制...
2021年3月最新民用航空管制频率表.pdf
03-29
资料收集自网上开放的2021年03版航图,包含各大机场,及中小机场航空管制频率表,同时旧版频率表废止!严禁用于非法用途，仅供学习参考研究使用。
2020年3月民用航空管制频率表.pdf
03-06
资料收集自网上开放的2020年01版航图,包含各大机场,及中小机场航空管制频率表,同时旧版频率表废止!严禁用于非法用途，仅供学习参考研究使用。
非旋treap线段树
03-05
指针构建非旋treap,代替线段树,可支持持久化,luogu3372可过.
[C++ 学习笔记] FHQ-Treap (无旋 Treap)
最新发布
lzl20220405的博客
01-05 1229
C++ FHQ-Treap - 无旋 Treap - 树堆 学习笔记
FHQ-Treap 简介
zhaoweiming2019的博客
08-11 754
FHQ-treap 即非旋Treap，是一种短小精悍，功能丰富的平衡树。据说它的效率介于 Treap 和 Splay 之间（可能是我的FHQ常数比较小，跑得比我的Treap还快）。它可以实现 Splay 可以实现的所有功能，包括平衡树的基本操作和区间翻转操作。它的实现难度比 Splay 要简单很多，没有 Splay 那么多转来转去的操作，不会令人头晕，而且FHQ代码的对称性良好，易于调试。有两种 FHQ-Treap，一种是按值为关键字的，一种是按下标为关键字的（适用于区间翻转等操作时）。...
fhq treap入门
qq_35577488的博客
03-23 655
目录一.前言二.原理三.具体实现模板题目一.洛谷P3369:普通平衡树二.洛谷P3391 文艺平衡树 一.前言 这个数据结构早有耳闻.只是一直没学.跟着AgOH学了一波.发现很简单. 二.原理 二叉搜索树之所以不够优秀,是因为存在一些数据能够将二叉搜索树退化成一条链.从而导致平方级的复杂度.那么对于这个缺陷,大佬们想出了很多方法去解决,衍生出了①替罪羊树②有旋/无旋Treap③splay.今天我入门学习的就是无旋treap. Treap之所以能够克服二叉搜索树的缺陷,主要因为它以随机的方式构造二叉搜索树.使
【算法学习】FHQ Treap （无旋Treap）
Pengwill's Blog
09-28 5072
FHQ Treap 简介 FHQ Treap和普通的Treap都是一个二叉搜索堆，其同时满足二叉树的性质（左子树的权值小于等于当前节点权值，右子树权值大于当前节点权值）和堆的性质（对于小根堆，当前节点的优先级是堆中最小的)。 FHQ Treap与一般的Treap的不同之处主要在于： 不用旋转，用split和merge来为维护堆的优先级。 能够可持久化。 操作 FHQ Treap的核心操作是sp...
从零开始的数据结构：FHQ-treap
m0_51791233的博客
10-01 342
这里是引用
无旋Treap——FHQ-Treap详解
Andysun06的博客
01-29 720
FHQ-treap——无旋treap 前置知识： 二叉排序树 treap 了解什么是函数式编程 引入 “平衡树好烦啊，转来转去的，而且treap遇到区间序列问题就gg。” “有一种不用旋转的treap，叫FHQ-treap，可以解决区间序列问题。” 正文 约定 我们做出以下约定： int ch[MAXN][3];//0 左孩子，1右孩子 int val[MAXN];//每个点的权值 int rnd[MAXN];//每个点的随机权值 int size[MAXN];//以每个点为根的树的大小
非旋转Treap-总结
linkfqy
04-18 2563
【前言】 非旋转Treap，又称FHQ_Treap，函数式Treap。由神犇FHQ首先提出（OrzOrzOrz）。由于使用了函数式编程的思想，非旋转Treap可以实现可持久化。 此外，非旋转Treap的大多数操作都基于分裂与合并，所以代码简短精炼，非常适合OIer食用。 【操作】 就讲讲分裂与合并好了……其他都是可以用脚趾头想出来的。分裂(split) 对于一棵Treap，我们把它分为ke
平衡树详解（fhq Treap）
Wu_while的博客
01-25 1126
你并不需要先学习Treap！

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值