BZOJ3224 洛谷3369 Tyvj 1728 普通平衡树 splay

本文介绍了一种使用伸展树(Splay Tree)解决的数据结构问题,提供了详细的代码实现及优化思路,包括插入、删除、查找等操作。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

欢迎访问~原文出处——博客园-zhouzhendong

去博客园看该题解


题目传送门 - BZOJ3224


题意概括

  您需要写一种数据结构(可参考题目标题),来维护一些数,其中需要提供以下操作:
1. 插入x数
2. 删除x数(若有多个相同的数,因只删除一个)
3. 查询x数的排名(若有多个相同的数,因输出最小的排名)
4. 查询排名为x的数
5. 求x的前驱(前驱定义为小于x,且最大的数)
6. 求x的后继(后继定义为大于x,且最小的数)


题解

  splay模板题。

  


代码

 UPD(2018-03-20):省选一试前夜,抱佛脚打板子,发现之前的代码真丑。注意,在splay中预先放一个非常大的数可以有效的排除平衡树空的特殊情况。我在细节处理上也增加了大量改进,代码量爆缩,运行效率略高了一些,老的那份代码在后面。

#include <bits/stdc++.h>
using namespace std;
const int N=100005;
int n,root=1,size=1,val[N],cnt[N],son[N][2],fa[N],tot[N];
int wson(int x){
	return son[fa[x]][1]==x;
}
void pushup(int x){
	tot[x]=cnt[x]+tot[son[x][0]]+tot[son[x][1]];
}
void rotate(int x){
	if (!x)
		return;
	int y=fa[x],z=fa[y],L=wson(x),R=L^1;
	if (z)
		son[z][wson(y)]=x;
	fa[x]=z,fa[y]=x,fa[son[x][R]]=y;
	son[y][L]=son[x][R],son[x][R]=y;
	pushup(y),pushup(x);
}
void splay(int x,int k){
	if (!x)
		return;
	if (!k)
		root=x;
	for (int y=fa[x];fa[x]!=k;rotate(x),y=fa[x])
		if (fa[y]!=k)
			rotate(wson(x)==wson(y)?y:x);
}
int find(int x,int v){
	return val[x]==v?x:find(son[x][v>val[x]],v);
}
int findkth(int x,int k){
	if (k<=tot[son[x][0]])
		return findkth(son[x][0],k);
	k-=tot[son[x][0]];
	if (k<=cnt[x])
		return x;
	k-=cnt[x];
	return findkth(son[x][1],k);
}
int findnxt(int x,int v){
	if (!x)
		return 0;
	if (val[x]<=v)
		return findnxt(son[x][1],v);
	else {
		int res=findnxt(son[x][0],v);
		return res?res:x;
	}
}
int findpre(int x,int v){
	if (!x)
		return 0;
	if (val[x]>=v)
		return findpre(son[x][0],v);
	else {
		int res=findpre(son[x][1],v);
		return res?res:x;
	}
}
void insert(int &x,int pre,int v){
	if (!x){
		x=++size;
		val[x]=v,cnt[x]=tot[x]=1,fa[x]=pre;
		splay(x,0);
		return;
	}
	tot[x]++;
	if (val[x]==v){
		cnt[x]++;
		return;
	}
	insert(son[x][v>val[x]],x,v);
}
void Insert(int v){insert(root,0,v);}
void Delete(int v){
	int x;
	splay(x=find(root,v),0);
	if (--cnt[x])
		return;
	splay(findnxt(root,v),root);
	root=son[x][1];
	son[root][0]=son[x][0];
	fa[son[x][0]]=root;
	fa[root]=son[x][0]=son[x][1]=0;
	pushup(root);
}
int Rank(int v){
	splay(find(root,v),0);
	return tot[son[root][0]]+1;
}
int main(){
	val[1]=2147483647;
	cnt[1]=tot[1]=1;
	scanf("%d",&n);
	while (n--){
		int opt,x;
		scanf("%d%d",&opt,&x);
		if (opt==1) Insert(x);
		if (opt==2) Delete(x);
		if (opt==3) printf("%d\n",Rank(x));
		if (opt==4) printf("%d\n",val[findkth(root,x)]);
		if (opt==5) printf("%d\n",val[findpre(root,x)]);
		if (opt==6) printf("%d\n",val[findnxt(root,x)]);
	}
	return 0;
}

  

 

#include <cstring>
#include <algorithm>
#include <cstdio>
#include <cmath>
#include <cstdlib>
using namespace std;
const int N=100005;
struct Splay{
	int fa[N],lc[N],rc[N],size[N],cnt[N],val[N],root,tot;
	void clear(){//splay³õʼ»¯ 
		root=tot=0,memset(size,0,sizeof size),memset(cnt,0,sizeof cnt);
	}
	void update(int x){//¼ÆËã(¸üÐÂ)ijһ¸ö½ÚµãµÄsize 
		if (x){
			size[x]=cnt[x];
			size[x]+=lc[x]?size[lc[x]]:0;
			size[x]+=rc[x]?size[rc[x]]:0;
		}
	}
	void zig(int x){//ÓÒÐý 
		if (!fa[x])
			return;
		int y=fa[x],z=fa[y];
		if (z)
			if (lc[z]==y)
				lc[z]=x;
			else
				rc[z]=x;
		fa[x]=z;
		fa[y]=x;
		fa[rc[x]]=y;
		lc[y]=rc[x];
		rc[x]=y;
		update(y);
		update(x);
	}
	void zag(int x){//×óÐý 
		if (!fa[x])
			return;
		int y=fa[x],z=fa[y];
		if (z)
			if (lc[z]==y)
				lc[z]=x;
			else
				rc[z]=x;
		fa[x]=z;
		fa[y]=x;
		fa[lc[x]]=y;
		rc[y]=lc[x];
		lc[x]=y;
		update(y);
		update(x);
	}
	void splay(int rt,int x){//splayÖ÷¹ý³Ì(ÔÚÒÔrtΪ¸ùµÄ×ÓÊ÷ÖÐsplay£¬×îÖÕx³ÉΪµ±Ç°×ÓÊ÷µÄÐÂ׿×Ú)
		rt=fa[rt];
		int y,z;
		while (fa[x]!=rt){
			y=fa[x];
			z=fa[y];
			if (z&&z!=rt){
				if (lc[z]==y&&lc[y]==x)
					zig(y),zig(x);
				else if (lc[z]==y&&rc[y]==x)
					zag(x),zig(x);
				else if (rc[z]==y&&lc[y]==x)
					zig(x),zag(x);
				else
					zag(y),zag(x);
			}
			else if (lc[y]==x)
				zig(x);
			else
				zag(x);
		}
		if (!fa[x])
			root=x;
	}
	void splay(int x){
		splay(root,x);
	}
	void ins(int &k,int key,int p){//ÔÚkÕâ¸ö×ÓÊ÷ÖвåÈëֵΪkeyµÄÊý£¬µ±Ç°½ÚµãµÄ¸¸Ç×Ϊp 
		if (!k){
			size[k=++tot]=cnt[k]=1;
			lc[k]=rc[k]=0;
			val[k]=key;
			fa[k]=p;
			splay(k);
			return;
		}
		size[k]++;
		if (key==val[k]){
			cnt[k]++;
			return;
		}
		if (key<val[k])
			ins(lc[k],key,k);
		else
			ins(rc[k],key,k);
	}
	void ins(int key){//ÔÚÕû¸ösplayÀïÃæ²åÈëָΪkeyµÄÊý 
		if (!root){
			size[root=++tot]=cnt[root]=1;
			lc[root]=rc[root]=0;
			val[root]=key;
			fa[root]=0;
			return;
		}
		ins(root,key,0);
	}
	int find(int k,int key){//ÔÚÒÔkΪ¸ùµÄsplayÀïÃæ²éÕÒֵΪkeyµÄÊýµÄ´æ´¢±àºÅ  
		int tmp=k;
		while (k)
			if (val[k]==key)
				break;
			else
				k=key<val[k]?lc[k]:rc[k];
		if (k)
			splay(k);
		return k;
	}
	int find(int key){//splayÀïÃæ²éÕÒֵΪkeyµÄÊýµÄ´æ´¢±àºÅ  
		return find(root,key);
	}
	void del(int key){//ɾ³ýֵΪkeyµÄÊý£¨µ«ÊÇtot²»¼õ£¬Òâζ×ÅÓÐ1¸ö¿Õ¼ä»áÀË·Ñ£© 
		int x=find(key),ls=lc[x],rs=rc[x];
		if (!x)
			return;
		if (--cnt[x])
			return;
		if (!ls&&!rs){
			clear();
			return;
		}
		if (!ls)
			root=rs,fa[rs]=0;
		else if (!rs)
			root=ls,fa[ls]=0;
		else {
			int lson=getmax(ls);
			swap(lson,ls);
			fa[lson]=0;
			splay(ls);
			rc[ls]=rs;
			fa[rs]=ls;
			update(ls);
		}
	}
	int getmin(int k){//ÔÚ¸ùΪkµÄsplayÖÐÕÒ×îСֵµÄ´æ´¢±àºÅ 
		return lc[k]?getmin(lc[k]):k;
	}
	int getmax(int k){//ÔÚ¸ùΪkµÄsplayÖÐÕÒ×î´óÖµµÄ´æ´¢±àºÅ 
		return rc[k]?getmax(rc[k]):k;
	}
	int findkth(int k){//ѰÕÒµÚkСµÄ 
		int t=root;
		while (t){
			if (size[lc[t]]<k&&k<=size[lc[t]]+cnt[t])
				break;
			if (k<=size[lc[t]])
				t=lc[t];
			else
				k-=size[lc[t]]+cnt[t],t=rc[t];
		}
		splay(t);
		return t;
	}
	int findpre(int k,int key){//ÔÚ¸ùΪkµÄsplayÖÐÕÒkeyµÄǰÇý£¨¼´Öµ<keyµÄ×î´óÖµ£©µÄ´æ´¢±àºÅ 
		if (!k)
			return 0;
		if (key<=val[k])
			return findpre(lc[k],key);
		else {
			int tmp=findpre(rc[k],key);
			return tmp?tmp:k;
		}
	}
	int findpre(int key){return findpre(root,key);}//ÕÒkeyµÄǰÇý£¨¼´Öµ<keyµÄ×î´óÖµ£©µÄ´æ´¢±àºÅ 
	int findnxt(int k,int key){//ÔÚ¸ùΪkµÄsplayÖÐÕÒkeyµÄºó¼Ì£¨¼´Öµ>keyµÄ×îСֵ£©µÄ´æ´¢±àºÅ 
		if (!k)
			return 0;
		if (key>=val[k])
			return findnxt(rc[k],key);
		else {
			int tmp=findnxt(lc[k],key);
			return tmp?tmp:k;
		}
	}
	int findnxt(int key){return findnxt(root,key);}//ÕÒkeyµÄºó¼Ì£¨¼´Öµ>keyµÄ×îСֵ£©µÄ´æ´¢±àºÅ 
}s;
int n,opt,x;
int main(){
	scanf("%d",&n);
	s.clear();
	for (int i=1;i<=n;i++){
		scanf("%d%d",&opt,&x);
		if (opt==1)	s.ins(x);
		if (opt==2)	s.del(x);
		if (opt==3){//ÇóֵΪxµÄÃû´Î 
			s.splay(s.find(x));
			printf("%d\n",s.size[s.lc[s.root]]+1);
		}
		if (opt==4)	printf("%d\n",s.val[s.findkth(x)]);
		if (opt==5)	printf("%d\n",s.val[s.findpre(x)]);
		if (opt==6)	printf("%d\n",s.val[s.findnxt(x)]);
	}
	return 0;
}

 

  

 

转载于:https://www.cnblogs.com/zhouzhendong/p/BZOJ3224.html

### 关于 BZOJ1728 Two-Headed Cows (双头牛) 的算法解析 此问题的核心在于如何通过有效的图论方法解决给定约束下的最大独立集问题。以下是详细的分析和解答。 #### 问题描述 题目要求在一个无向图中找到最大的一组节点集合,使得这些节点之间满足特定的颜色匹配条件。具体来说,每条边连接两个节点,并附带一种颜色标记(A 或 B)。对于任意一条边 \(u-v\) 和其对应的颜色 \(c\),如果这条边属于最终选取的子集中,则必须有至少一个端点未被选入该子集或者两端点均符合指定颜色关系。 #### 解决方案概述 本题可以通过 **二分枚举 + 图染色验证** 来实现高效求解。核心思想如下: 1. 假设当前最优解大小为 \(k\),即尝试寻找是否存在一个大小为 \(k\) 的合法子集。 2. 枚举每一个可能作为起点的节点并将其加入候选子集。 3. 对剩余部分执行基于 BFS/DFS 的图遍历操作,在过程中动态调整其他节点的状态以确保整体合法性。 4. 如果某次试探能够成功构建符合条件的大规模子集,则更新答案;反之则降低目标值重新测试直至收敛至最佳结果。 这种方法利用了贪心策略配合回溯机制来逐步逼近全局最优点[^1]。 #### 实现细节说明 ##### 数据结构设计 定义三个主要数组用于记录状态信息: - `color[]` : 存储每个顶点所分配到的具体色彩编号; - `used[]`: 表示某个定点是否已经被处理过; - `adjList[][]`: 记录邻接表形式表示的原始输入数据结构便于后续访问关联元素。 ##### 主要逻辑流程 ```python from collections import deque def check(k, n): def bfs(start_node): queue = deque([start_node]) used[start_node] = True while queue: u = queue.popleft() for v, c in adjList[u]: if not used[v]: # Assign opposite color based on edge constraint 'c' target_color = ('B' if c == 'A' else 'A') if color[u]==c else c if color[v]!=target_color and color[v]!='?': return False elif color[v]=='?': color[v]=target_color queue.append(v) used[v] =True elif ((color[u]==c)==(color[v]==('B'if c=='A'else'A'))): continue return True count=0 success=True for i in range(n): if not used[i]: temp_count=count+int(color[i]=='?' or color[i]=='A') if k<=temp_count: color_copy=color[:] if bfs(i): count=temp_count break else : success=False return success n,m=list(map(int,input().split())) colors=[['?']*m]*n for _ in range(m): a,b,c=input().strip().split() colors[int(a)-1].append((int(b),c)) low ,high,res=0,n,-1 while low<=high: mid=(low+high)//2 color=['?']*n used=[False]*n if check(mid,n): res=mid low=mid+1 else : high=mid-1 print(res) ``` 上述代码片段展示了完整的程序框架以及关键函数 `check()` 的内部运作方式。它接受参数 \(k\) 并返回布尔值指示是否有可行配置支持如此规模的选择[^2]。 #### 复杂度分析 由于采用了二分查找技术缩小搜索空间范围再加上单轮 DFS/BFS 时间复杂度 O(V+E),总体性能表现良好适合大规模实例运行需求。 ---
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值