【SSL】1612最优布线问题(最小生成树)

【SSL】1612最优布线问题(最小生成树)

Time Limit:1000MS
Memory Limit:65536K

Description

学校有n台计算机，为了方便数据传输，现要将它们用数据线连接起来。两台计算机被连接是指它们之间有数据线连接。由于计算机所处的位置不同，因此不同的两台计算机的连接费用往往是不同的。
当然，如果将任意两台计算机都用数据线连接，费用将是相当庞大的。为了节省费用，我们采用数据的间接传输手段，即一台计算机可以间接的通过若干台计算机（作为中转）来实现与另一台计算机的连接。
现在由你负责连接这些计算机，你的任务是使任意两台计算机都连通（不管是直接的或间接的）。

Input

输入文件wire.in，第一行为整数n（2<=n<=100），表示计算机的数目。此后的n行，每行n个整数。第x+1行y列的整数表示直接连接第x台计算机和第y台计算机的费用。

Output

输出文件wire.out，一个整数，表示最小的连接费用。

Sample Input

3
0 1 2
1 0 1
2 1 0

Sample Output

2

（注：表示连接1和2，2和3，费用为2）

Hint

Fillchar(f,sizeof(f),$7f)把所有值赋值为最大

思路

prim

Prim算法采用与Dijkstra、Bellman-Ford算法一样的“蓝白点”思想：白点代表已经进入最小生成树的点，蓝点代表未进入最小生成树的点。
算法描述：
以1为起点生成最小生成树，min[v]表示蓝点v与白点相连的最小边权。
MST表示最小生成树的权值之和。
a）初始化：min[v]= ∞(v≠1); min[1]=0;MST=0;
b）for (i = 1; i<= n; i++)
1.寻找min[u]最小的蓝点u。
2.将u标记为白点
3.MST+=min[u]
4.for 与白点u相连的所有蓝点v
if (w[u][v]<min[v])
min[v]=w[u][v];
c）算法结束： MST即为最小生成树的权值之和

算法分析&思想讲解：
Prim算法每次循环都将一个蓝点u变为白点，并且此蓝点u与白点相连的最小边权min[u]还是当前所有蓝点中最小的。这样相当于向生成树中添加了n-1次最小的边，最后得到的一定是最小生成树。
算法时间复杂度：O (N²)。
可以用邻接矩阵和邻接表。

kruskal

假设连通网G＝（V，E），则令最小生成树的初始状态为只有n个顶点而无边的非连通网，图中每个顶点自成一个连通分量。在E中选择代价最小的边，若该边依附的顶点落在T中不同的连通分量上，则将此边加入到T中，否则舍去此边而选择下一条代价最小的边。依此类推，直至T中所有顶点都在同一连通分量上为止。
Kruskal算法将一个连通块当做一个集合。Kruskal首先将所有的边按从小到大顺序排序（一般使用快排），并认为每一个点都是孤立的，分属于n个独立的集合。然后按顺序枚举每一条边。如果这条边连接着两个不同的集合，那么就把这条边加入最小生成树，这两个不同的集合就合并成了一个集合；如果这条边连接的两个点属于同一集合，就跳过。直到选取了n-1条边为止。

算法描述:
初始化并查集。father[x]=x。
tot=0
将所有边用快排从小到大排序。
计数器 k=0;
for (i=1; i<=M; i++) //循环所有已从小到大排序的边
if 这是一条u,v不属于同一集合的边(u,v)(因为已经排序，所以必为最小)，
begin
　①合并u,v所在的集合，相当于把边(u,v)加入最小生成树。
　 ②tot=tot+W(u,v)
③k++
④如果k=n-1,说明最小生成树已经生成，则break;
end;
结束，tot即为最小生成树的总权值之和。

Kruskal每次都选择一条最小的边，而且这条边的两个顶点分属于两个不同的集合。将选取的这条边加入最小生成树，并且合并集合。

Kruskal算法的时间复杂度为O(E*logE)，E为边数。

代码

prim邻接矩阵

#include<iostream>
#include<cstdio>
#include<queue>
#include<cstring>
#include<algorithm>
using namespace std;
int n,a[110][110],ans=0,tot=0,mn[110];
bool d[110];
void input()
{
	int i,j;
	memset(mn,100,sizeof(mn));
	scanf("%d",&n);
	for(i=1;i<=n;i++)
		for(j=1;j<=n;j++)
			scanf("%d",&a[i][j]);
	return;
}
void prim()
{
	int i,j,w;
	memset(d,true,sizeof(d));
	mn[1]=0;
	for(i=1;i<n;i++)
	{
		for(w=0,j=1;j<=n;j++)//找最近蓝点
			if(d[j]&&mn[j]<mn[w])
				w=j;
		d[w]=0;
		for(j=1;j<=n;j++)//修改最短距离
			if(d[j]&&a[w][j]<mn[j])
				mn[j]=a[w][j];
	}
	for(i=1;i<=n;i++)//累计距离和
		ans+=mn[i];
	return;
}
int main()
{
	input();
	prim();
	printf("%d",ans);
	return 0;
}

prim邻接表

#include<iostream>
#include<cstdio>
#include<queue>
#include<cstring>
#include<algorithm>
using namespace std;
int n,a[110][110],ans=0,tot=0,mn[110],head[110];
bool d[110];
struct jgt
{
	int x,y,w,nst;
}f[10010];
void input()
{
	int i,j;
	memset(mn,100,sizeof(mn));
	memset(head,0,sizeof(head));
	scanf("%d",&n);
	for(i=1;i<=n;i++)
		for(j=1;j<=n;j++)
		{
			scanf("%d",&a[i][j]);
			if(a[i][j])	//建边
			{
				tot++;
				f[tot].x=i;
				f[tot].y=j;
				f[tot].w=a[i][j];
				f[tot].nst=head[i];
				head[i]=tot;
			}
		}
	return;
}
void prim()
{
	int i,j,w;
	memset(d,true,sizeof(d));
	mn[1]=0;
	for(i=1;i<n;i++)
	{
		for(w=0,j=1;j<=n;j++)//找最近蓝点
			if(d[j]&&mn[j]<mn[w])
				w=j;
		d[w]=0;
		for(j=head[w];j;j=f[j].nst)//修改最短距离
			if(d[f[j].y]&&a[w][f[j].y]<mn[f[j].y])
				mn[f[j].y]=a[w][f[j].y];
	}
	for(i=1;i<=n;i++)//累计距离和
		ans+=mn[i];
	return;
}
int main()
{
	input();
	prim();
	printf("%d",ans);
	return 0;
}

kruskal

#include<iostream>
#include<cstdio>
#include<queue>
#include<cstring>
#include<algorithm>
using namespace std;
int n,a[110][110],ans=0,tot=0,b[110];
bool d[110];
struct jgt
{
	int x,y,w;
}f[10010];
int bcj(int x)//并查集
{
	if(b[x]==x)
		return x;
	b[x]=bcj(b[x]);
	return b[x];
}
bool cmp(jgt t1,jgt t2)
{
	return t1.w<t2.w;
}
void input()
{
	int i,j;
	scanf("%d",&n);
	for(i=1;i<=n;b[i]=i,i++)
		for(j=1;j<=n;j++)
		{
			scanf("%d",&a[i][j]);
			if(i<j)
			{
				f[tot].x=i;
				f[tot].y=j;
				f[tot++].w=a[i][j];
			}
		}
	return;
}
void kruskal()
{
	int i,x;
	sort(f,f+tot,cmp);//排序
	for(i=0;n;i++)
	{
		if(bcj(f[i].x)!=bcj(f[i].y))//若不连通
		{
			n--;
			ans+=f[i].w;
			b[bcj(f[i].x)]=bcj(f[i].y);//合并
		}
	}
	return;
}
int main()
{
	input();
	kruskal();
	printf("%d",ans);
	return 0;
}