bzoj2816.。。。。zjoi2012网络

本文深入探讨了Splay树森林的应用及其实现细节,包括旋转操作、路径压缩等关键概念,并介绍了如何通过Splay树解决复杂的数据结构问题。

水水的splay  准确说是splay森林


其实不止颜色数的splay树  (还有一些单点也可以存在)


不过都可以用 一个 [C][N]  存下来


其实splay之间的边  都不是实际上连起来的边

每个点其实实际上是和它的前驱后继连起来的


不然zigzag就会破坏结果  原因自行脑补


还有要注意有时候需要一个reverse  


要连接的两条链   需要 一个的   节点 的d 儿子 

另一个节点的 d ^ 1 儿子

均为空  即等于0


才可以连接  否则为非法情况


剩下的就靠码力了


拒绝眼残 拒绝手残。。


就好了








#include <cstdio>
#include <map>
#define lp(i,j,k) for(int i = j;i <= k;++i)
#define CH(i,j,k) ch[i][j][k]
#define FA(i,j) fa[i][j]
#define D 10010
using namespace std;
map<pair<int,int>,int>MAP;
int a[D],n,m,c,k,u,v,w,fa[11][D],x[11][D],ch[11][D][2];
bool lazy[11][D];
void reverse (int col,int X) {
	int temp = CH(col,X,1);
	CH(col,X,1) = CH(col,X,0);
	CH(col,X,0) = temp;
	if(X)
		lazy[col][X] ^= 1;
}
void push_down (int col,int X) {
	if(lazy[col][X]) {
		reverse(col,CH(col,X,0));
		reverse(col,CH(col,X,1));
	}
	lazy[col][X] = 0;
}
void push_up (int col,int X) {
	x[col][X] = max(x[col][CH(col,X,0)],x[col][CH(col,X,1)]);
	x[col][X] = max(x[col][X],a[X]);
}
void rot (int col,int X,int d) {
	x[col][X] = x[col][FA(col,X)];
	x[col][FA(col,X)] = max(x[col][CH(col,X,d)],a[FA(col,X)]);
	x[col][FA(col,X)] = max(x[col][FA(col,X)],x[col][CH(col,FA(col,X),d)]);
	CH(col,FA(col,X),d ^ 1) = CH(col,X,d);
	if(CH(col,X,d))
		FA(col,CH(col,X,d)) = FA(col,X);
	CH(col,X,d) = FA(col,X);
	int temp = FA(col,FA(col,X));
	FA(col,FA(col,X)) = X;
	FA(col,X) = temp;
	if(!temp)
		return;
	if(CH(col,temp,1) == CH(col,X,d))
		CH(col,temp,1) = X;
	else CH(col,temp,0) = X;
}
void splay (int col,int X) {
	while(FA(col,X)) {
		if(FA(col,FA(col,X))) {
			push_down(col,FA(col,FA(col,X)));
			push_down(col,FA(col,X));
			push_down(col,X);
			int d =  CH(col,FA(col,FA(col,X)),0) == FA(col,X);
			if(CH(col,FA(col,X),d) == X)
				rot(col,X,d ^ 1),rot(col,X,d);
			else rot(col,X,d),rot(col,X,d);
		}
		else {
			push_down(col,FA(col,X));
			push_down(col,X);
			if(CH(col,FA(col,X),1) == X)
				rot(col,X,0);
			else rot(col,X,1);
		}
	}
}
int main () {
	scanf("%d%d%d%d",&n,&m,&c,&k);
	lp(i,1,n) {
		scanf("%d",&a[i]);
		lp(j,1,c)
			x[j][i] = a[i];
	}
	lp(i,1,m) {
		scanf("%d%d%d",&u,&v,&w);
		++w;
		MAP[make_pair(u,v)] = w;
		MAP[make_pair(v,u)] = w;
		splay(w,v);
		splay(w,u);
		int d;
		if(CH(w,u,0))
			d = 1;
		else d = 0;
		if(CH(w,v,d ^ 1))
			reverse(w,v);
		CH(w,u,d) = v;
		FA(w,v) = u;
		push_up(w,u);
		splay(w,v);
	}
	lp(i,1,k) {
		scanf("%d",&w);
		if(w == 0) {
			scanf("%d%d",&u,&v);
			a[u] = v;
			lp(j,1,c) {
				splay(j,u);
				x[j][u] = max(x[j][u],v);
			}
		}
		else {
			if(w == 1) {
				scanf("%d%d%d",&u,&v,&w);
				++w;
				int col = MAP[make_pair(u,v)];
				if(!col)
					printf("No such edge.\n");
				else {
					if(col == w) {
						printf("Success.\n");
						continue;
					}
					else {
						splay(w,u);
						splay(w,v);
						if((CH(w,u,0) && CH(w,u,1)) || (CH(w,v,0) && CH(w,v,1))) {
							printf("Error 1.\n");
							continue;
						}
						bool f = 0;
						lp(hh,0,1) {							
							if(CH(w,v,hh) == u && CH(w,v,hh ^ 1) == 0) {
								f = printf("Error 2.\n"),1;
								break;
							}
						}
						if(f)
							continue;
						splay(col,u);
						splay(col,v);
						int d;
						if(CH(col,v,0) == u)
							CH(col,v,0) = 0;
						else CH(col,v,1) = 0;
						FA(col,u) = 0;
						push_up(col,u);
						push_up(col,v);
						if(CH(w,u,0))
							d = 1;
						else d = 0;
						if(CH(w,v,d ^ 1))
							reverse(w,v);
						CH(w,u,d) = v;
						FA(w,v) = u;
						push_up(w,u);
						MAP[make_pair(u,v)] = w;
						MAP[make_pair(v,u)] = w;
						printf("Success.\n");
					}
				}
			}
			else {
				scanf("%d%d%d",&w,&u,&v);
				++w;
				splay(w,u);
				splay(w,v);
				int d = CH(w,v,1) == u;
				if(v != u) {
					if(d) {
						int pr = max(x[w][CH(w,u,0)],a[v]);
						pr = max(pr,a[u]);
						printf("%d\n",pr);
					}
					else {
						if(CH(w,v,0) == u) {
							int pr = max(x[w][CH(w,u,1)],a[v]);
							pr = max(pr,a[u]);
							printf("%d\n",pr);
						}
						else printf("-1\n");
					}
				}
				else printf("%d\n",a[u]);
			}
		}
	}
}













### 关于 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、付费专栏及课程。

余额充值