并查集

前言

一直以来,并查集学浅了……除了最基本的路径压缩什么都不会


路径压缩

基本操作

详细理论见蓝书。

并查集的初始化

for(int i=1;i<=n;i++) f[i]=i;

并查集的查询

int find(int x){return f[x]==x?x:f[x]=find(f[x]);}

并查集的合并

void merge(int x,int y){
  int f1=find(x),f2=find(y);
  if(f1!=f2) 不在一个集合中
  f[f1]=f2;
}

[NOI2015] 程序自动分析

题目

描述

在实现程序自动分析的过程中,常常需要判定一些约束条件是否能被同时满足。

考虑一个约束满足问题的简化版本:假设 \(x_1,x_2,x_3...\) 代表程序中出现的变量,给定 \(n\) 个形如 \(x_i=x_j\)\(x_i≠x_j\) 的变量相等/不相等的约束条件,请判定是否可以分别为每一个变量赋予恰当的值,使得上述所有约束条件同时被满足。例如,一个问题中的约束条件为:\(x_1=x_2,x_2=x_3 ,x_3=x_4 , x4≠x_1\) ,这些约束条件显然是不可能同时被满足的,因此这个问题应判定为不可被满足。

现在给出一些约束满足问题,请分别对它们进行判定。

输入格式

输入的第一行包含一个正整数 \(t\) 表示需要判定的问题个数。注意这些问题之间是相互独立的。

对于每个问题,包含若干行:

第一行包含一个正整数 \(n\) ,表示该问题中需要被满足的约束条件个数。接下来 \(n\) 行,每行包括三个整数 \(i,j,e\) ,描述一个相等/不等的约束条件,相邻整数之间用单个空格隔开。若 \(e=1\) ,则该约束条件为 \(x_i=x_j\)\(e=0\) ,则该约束条件为 \(x_i≠x_j\)

输出格式

输出包括 \(t\) 行。

输出文件的第 \(k\) 行输出一个字符串 \(YES\) 或者 \(NO\)(字母全部大写),\(YES\) 表示输入中的第 \(k\) 个问题判定为可以被满足,\(NO\) 表示不可被满足。

思路

首先,把有相等关系的变量合并,每次出现不等关系时,查询 \(i,j\) 是否在同一个集合中,若在,则矛盾。

由于题目的数据范围是 $i,j≤ 10^9 $ 所以需要离散化。

code

#include<bits/stdc++.h>
using namespace std;
const int N=1e6+105;
int f[N],u[N<<2],T,n,cnt;
bool flag;
struct qwq{int x,y,e;}q[N];
bool cmp(qwq a,qwq b){return a.e>b.e;}
inline void clean(){
	flag=1;cnt=-1;
	memset(u,0,sizeof(u));
	memset(q,0,sizeof(q));
	memset(f,0,sizeof(f));
}
void pre(int k){for(int i=1;i<=k;i++)f[i]=i;}
int find(int x){return x==f[x]?x:f[x]=find(f[x]);}
void check(int x,int y,int e){
	int f1=find(x),f2=find(y);
	if(e==1) f[f1]=f2;
	else if(f1==f2){flag=0;printf("NO\n");}
}
inline int read(){
	int x=0,f=1;char ch=getchar();
	while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
	while(isdigit(ch)){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
	return x*f;
}
int main(){
	T=read();
	while(T--){
		n=read();clean();
		for(int i=1;i<=n;i++){
			q[i].x=read();q[i].y=read();q[i].e=read();
			u[++cnt]=q[i].x;u[++cnt]=q[i].y;
		}
		sort(u,u+cnt);int pwp=unique(u,u+cnt)-u;
		for(int i=1;i<=n;i++){
			q[i].x=lower_bound(u,u+pwp,q[i].x)-u;
			q[i].y=lower_bound(u,u+pwp,q[i].y)-u;
		}
		pre(pwp);sort(q+1,q+n+1,cmp);
		for(int i=1;i<=n;i++){
			check(q[i].x,q[i].y,q[i].e);
			if(!flag) break;
		}
		if(flag) printf("YES\n");
	}
	return 0;
}

超市

题目

思路

code


扩展域和边带权

建立一个数组 \(d\) ,\(d[x]\) 记录战舰 \(x\)\(f[x]\) 之间的边的权值。在路径压缩把 \(x\) 指向树根的同事,把 \(d[x]\) 更新为从 \(x\) 到树根的路径上的所有边权之和。

建立一个数组 \(siz\) 在每个树根上记录集合大小。

查询+路径压缩

int find(int x){
	if(x==f[x]) return x;
	int root=find(f[x]);
	d[x]+=d[f[x]];
	return f[x]=root;
}

合并

void merge(int x,int y){
	x=find(x);y=find(y);
	f[x]=y;d[x]=siz[y];
	siz[y]+=siz[x];
}

[NOI2002] 银河英雄传说

题目

描述

公元 \(5801\) 年,地球居民迁至金牛座 \(\alpha\) 第二行星,在那里发表银河联邦创立宣言,同年改元为宇宙历元年,并开始向银河系深处拓展。

宇宙历 \(799\) 年,银河系的两大军事集团在巴米利恩星域爆发战争。泰山压顶集团派宇宙舰队司令莱因哈特率领十万余艘战舰出征,气吞山河集团点名将杨威利组织麾下三万艘战舰迎敌。

杨威利擅长排兵布阵,巧妙运用各种战术屡次以少胜多,难免恣生骄气。在这次决战中,他将巴米利恩星域战场划分成 \(30000\) 列,每列依次编号为 \(1, 2,\ldots ,30000\) 。之后,他把自己的战舰也依次编号为 \(1, 2, \ldots , 30000\) ,让第 \(i\) 号战舰处于第 \(i\) 列,形成“一字长蛇阵”,诱敌深入。这是初始阵形。当进犯之敌到达时,杨威利会多次发布合并指令,将大部分战舰集中在某几列上,实施密集攻击。合并指令为 M i j ,含义为第 \(i\) 号战舰所在的整个战舰队列,作为一个整体(头在前尾在后)接至第 \(j\) 号战舰所在的战舰队列的尾部。显然战舰队列是由处于同一列的一个或多个战舰组成的。合并指令的执行结果会使队列增大。

然而,老谋深算的莱因哈特早已在战略上取得了主动。在交战中,他可以通过庞大的情报网络随时监听杨威利的舰队调动指令。

在杨威利发布指令调动舰队的同时,莱因哈特为了及时了解当前杨威利的战舰分布情况,也会发出一些询问指令: C i j 。该指令意思是,询问电脑,杨威利的第 \(i\) 号战舰与第 \(j\) 号战舰当前是否在同一列中,如果在同一列中,那么它们之间布置有多少战舰。

作为一个资深的高级程序设计员,你被要求编写程序分析杨威利的指令,以及回答莱因哈特的询问。

输入格式

第一行有一个整数 \(T\) \((1 \le T \le 5 \times 10^5 )\) ,表示总共有 \(T\) 条指令。

以下有 \(T\) 行,每行有一条指令。指令有两种格式:

M i j\(i\)\(j\) 是两个整数\((1 \le i,j \le 30000)\),表示指令涉及的战舰编号。该指令是莱因哈特窃听到的杨威利发布的舰队调动指令,并且保证第 \(i\) 号战舰与第 \(j\) 号战舰不在同一列。

C i j\(i\)\(j\) 是两个整数 \((1 \le i,j \le 30000 )\) ,表示指令涉及的战舰编号。该指令是莱因哈特发布的询问指令。

输出格式
依次对输入的每一条指令进行分析和处理:

  • 如果是杨威利发布的舰队调动指令,则表示舰队排列发生了变化,你的程序要注意到这一点,但是不要输出任何信息。

  • 如果是莱因哈特发布的询问指令,你的程序要输出一行,仅包含一个整数,表示在同一列上,第 \(i\) 号战舰与第 \(j\) 号战舰之间布置的战舰数目。如果第 \(i\) 号战舰与第 \(j\) 号战舰当前不在同一列上,则输出 \(-1\)

code

#include<bits/stdc++.h>
using namespace std;
const int N=1e5+105;
int f[N],d[N],siz[N];
int find(int x){
	if(x==f[x]) return x;
	int root=find(f[x]);
	d[x]+=d[f[x]];
	return f[x]=root;
}
void merge(int x,int y){
	x=find(x);y=find(y);
	f[x]=y;d[x]=siz[y];
	siz[y]+=siz[x];
}
inline int read(){
	int x=0,f=1;char ch=getchar();
	while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
	while(isdigit(ch)){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
	return x*f;
}
int T,x,y;
char c;
int main(){
	T=read();
	for(int i=1;i<=N-105;i++) f[i]=i,siz[i]=1;
	while(T--){
		cin>>c;x=read();y=read();
		if(c=='M')merge(x,y);
		else if(c=='C'){
			int X=find(x),Y=find(y);
			if(X!=Y)printf("-1\n");
			if(X==Y) printf("%d\n",abs(d[x]-d[y])-1);			
		}
	}
	return 0;
}

关于爆空间导致的结果是 \(WA\) 所以瞎调了好长时间这档事

预处理循环一开始写的是 i<=N 然后因为数组是 int f[N] 导致 To short on line 1


奇偶游戏

或者这里: [CEOI1999]Parity Game

题目

思路

边带权

\(sum\) 表示前缀和。

  1. \(S[l\) ~ \(r]\) 有偶数个 \(1\),等价于 \(sum[l-1]\)\(sum[r]\) 的奇偶性相同

  2. \(S[l\) ~ \(r]\) 有奇数个 \(1\),等价于 \(sum[l-1]\)\(sum[r]\) 的奇偶性不同

那么可以转换为:

  • \(x_1\)\(x_2\) 的奇偶性相同, \(x_2\)\(x_3\) 的奇偶性也相同,那么 \(x_1\)\(x_3\) 的奇偶性相同

  • \(x_1\)\(x_2\) 的奇偶性相同, \(x_2\)\(x_3\) 的奇偶性不同,那么 \(x_1\)\(x_3\) 的奇偶性不同

(若 \(x_3\)\(x_2\) 的奇偶性不同,\(x_2\)\(x_1\) 的奇偶性相同,那么 \(x_1\)\(x_3\) 的奇偶性不同)

  • \(x_1\)\(x_2\) 的奇偶性不同, \(x_2\)\(x_3\) 的奇偶性也不同,那么 \(x_1\)\(x_3\) 的奇偶性相同

\(d[x]=0\) 表示 \(x\)\(f[x]\) 的奇偶性相同。

\(d[x]=1\) 表示 \(x\)\(f[x]\) 的奇偶性不同。

检查每一个询问 \((x,y,v)\)\(x\)\(y\)是否在同一个集合内

code

边带权做法:

#include<bits/stdc++.h>
using namespace std;
const int N=5e4+105;
struct qwq{int l,r,v;}ask[N];
int a[N],f[N],d[N],n,m,tot;
inline int read(){
	int x=0,f=1;char ch=getchar();
	while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
	while(isdigit(ch)){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
	return x*f;
}
inline void lsh(){
	n=read();m=read();
	for(int i=1;i<=m;i++){
		char s[5];
		ask[i].l=read();ask[i].r=read();	
		scanf("%s",s);ask[i].v=(s[0]=='o'?1:0);
		a[++tot]=ask[i].l-1;a[++tot]=ask[i].r;
	}
	sort(a+1,a+1+tot);n=unique(a+1,a+tot+1)-a-1;
}
int find(int x){
	if(x==f[x]) return x;
	int root=find(f[x]);d[x]^=d[f[x]];
	return f[x]=root;
}
int main(){
	lsh();
	for(int i=1;i<=n;i++) f[i]=i;
	for(int i=1;i<=m;i++){
		int x=lower_bound(a+1,a+1+n,ask[i].l-1)-a;
		int y=lower_bound(a+1,a+1+n,ask[i].r)-a;
		int X=find(x),Y=find(y);
		if(X==Y){
			if((d[x]^d[y])!=ask[i].v)
				{printf("%d\n",i-1);return 0;}
		}
		else f[X]=Y,d[X]=d[x]^d[y]^ask[i].v;					
	}
	printf("%d",m);
}

[NOI2001] 食物链

题目

思路

code

[NOIP2010 提高组] 关押罪犯

题目

思路

code

posted @ 2021-09-12 08:51  Nickle-NI  阅读(39)  评论(0编辑  收藏  举报