连通问题的几种算法.c

连通问题(connectivity problem)是一种图论问题.在可以把问题归结为建立两元素间联系的数学问题中，思维图的运算可以归结为安插新的结点和边，最终完成包含表示这两元素结点的一个连通图。

程序1.1 连通问题的快速查找算法

这个程序从标准输入读取小于N的非负整数对序列（对p-q表示“把对象B连接到q”），并且输出还未联通的输入对。程序中使用数组id，每个元素代表一个对象，且具有以下性质，当且仅当p和q是连通的，id[p]和id[q]相等。为简化起见，定义N为编译时的常数。另一方面，也可以从输入得到它，并动态地为它分配id数组。

//   Date:2020/4/11
//   Author:xiezhg5
#include <stdio.h>
#define N 10000
int main(void)
{
	int i,p,q,t,id[N];
	for(i=0;i<N;i++)
		id[i]=i;
	while(scanf("%d %d\n",&p,&q)==2)
	{
		if(id[p]==id[q])
			continue;
		for(t=id[p],i=0;i<N;i++)
		{
			if(id[i]==t)
				id[i]=id[q];
		}
		printf("%d %d\n",p,q);
	}
	return 0;
}

程序1.2 连通问题的快速合并算法

如果用这段代码替换程序1.1中的while循环体，我们就得到一个同样满足程序1.1中说明的程序，但其合并操作的计算时间较短，查找操作的计算时间更长。这个代码中的for循环以及后面的if语句指明了在数组id中p和q是连通的充分必要条件。赋值语句id[i]=j实现合并操作。

//   Date:2020/4/11
//   Author:xiezhg5
#include <stdio.h>
#define N 10000
int main(void)
{
	int i,p,q,t,id[N];
	for(i=0;i<N;i++)
		id[i]=i;
	while(scanf("%d %d\n",&p,&q)==2)
	{
		for(i=p;i!=id[i];i=id[i]);
		for(t=q;t!=id[t];t=id[t]);
		if(i==t)
		continue;
		id[i]=t;
		printf("%d %d\n",p,q);
	}
	return 0;
}

程序1.3 加权快速合并算法

这个程序是快速合并算法的一个改进。用另一个数组sz记录每个id[i]==i的对象所在树中的节点数，使得合并操作能够吧较小的树连接到较大的树上，以防止树中长路径的增长。

//   Date:2020/4/11
//   Author:xiezhg5
#include <stdio.h>
#define N 10000
int main(void)
{
	int i,p,q,t,id[N],sz[N];
	for(i=0;i<N;i++)
	{
		id[i]=i;
		sz[i]=1;
	}
	while(scanf("%d %d\n",&p,&q)==2)
	{
		for(i=p;i!=id[i];i=id[i]);
		for(t=q;t!=id[t];t=id[t]);
		if(i==t)
		continue;
		if(sz[i]<sz[t])
		{
			id[i]=t;
			sz[t]+=sz[i];
		}
		else
		{
			id[t]=i;
			sz[i]+=sz[t];
		}
		printf("%d %d\n",p,q);
	}
	return 0;
}

程序1.4 等分路径压缩

如果我们用这段代码替换程序1.3中的for循环，就等分了遍历的任何路径的长度。改变的结果是经过的长序列的合并操作之后，树几乎完全是扁平的。

//   Date:2020/4/11
//   Author:xiezhg5
#include <stdio.h>
#define N 10000
int main(void)
{
	int i,p,q,t,id[N],sz[N];
	for(i=0;i<N;i++)
	{
		id[i]=i;
		sz[i]=1;
	}
	while(scanf("%d %d\n",&p,&q)==2)
	{
		for(i=p;i!=id[i];i=id[i])
		id[i]=id[id[i]];
		for(t=q;t!=id[t];t=id[t])
		id[t]=id[id[t]];
		printf("%d %d\n",p,q);
	}
	return 0;
}