图论算法----最小生成树Prim算法详解

一、题目描述

最短网络

题目描述

农民约翰被选为他们镇的镇长！他其中一个竞选承诺就是在镇上建立起互联网，并连接到所有的农场。当然，他需要你的帮助。

约翰已经给他的农场安排了一条高速的网络线路，他想把这条线路共享给其他农场。为了用最小的消费，他想铺设最短的光纤去连接所有的农场。

你将得到一份各农场之间连接费用的列表，你必须找出能连接所有农场并所用光纤最短的方案。每两个农场间的距离不会超过100000

输入

第一行： 农场的个数，N（3<=N<=100）。

第二行..结尾: 后来的行包含了一个N*N的矩阵,表示每个农场之间的距离。理论上，他们是N行，每行由N个用空格分隔的数组成，实际上，他们限制在80个字符，因此，某些行会紧接着另一些行。当然，对角线将会是0，因为不会有线路从第i个农场到它本身。

输出

只有一个输出，其中包含连接到每个农场的光纤的最小长度。

样例输入

4
0 4 9 21
4 0 8 17
9 8 0 16
21 17 16 0

样例输出

28

二、算法分析

Prim算法是一个经典的最小生成树算法，有点像Dijkstra算法，都用的是贪心的策略，

下面是Dijkstra的算法主体和Prim的算法主体的对比：

Dijkstra:

for(i=1;i<=n;i++)
dis[i]=a[q][i];
b[q]=1;
dis[q]=0;
for(i=1;i<=n-1;i++)

{
xi=10000000;
k=0;
for(j=1;j<=n;j++)
if(b[j]==0&&dis[j]<xi){

xi=dis[j];
k=j;  
}
b[k]=1;
for(j=1;j<=n;j++)

if(dis[j]>dis[k]+a[k][j])
dis[j]=dis[k]+a[k][j];
}

Prim:

dis[1]=0;
for(i=1;i<=n;i++)

{
k=0;
for(j=1;j<=n;j++)
if(b[j]==0&&dis[j]<dis[k])
k=j;
b[k]=1;
for(j=1;j<=n;j++){

if(b[j]==0&&a[k][j]<dis[j])
dis[j]=a[k][j];
}
}

for(i=1;i<=n;i++)
sum+=dis[i];

可以看出Prim的初始化比Dijkstra的初始化简单得多，只是第二个循环和第三个循环的判断条件不同了，总体结构差不多。

AC代码：

#include<cstdio>
#include<cstring>
int a[105][105],dis[105],sum;
bool b[105];
int main()
{
    int n,i,j,k;
    scanf("%d",&n);
    for(i=1;i<=n;i++)
        for(j=1;j<=n;j++)
            scanf("%d",&a[i][j]);
    memset(dis,0x7f,sizeof(dis));
    dis[1]=0;
    for(i=1;i<=n;i++){
        k=0;
        for(j=1;j<=n;j++)
            if(b[j]==0&&dis[j]<dis[k])
                k=j;
        b[k]=1;
        for(j=1;j<=n;j++){
            if(b[j]==0&&a[k][j]<dis[j])
                dis[j]=a[k][j];
        }
    }
    for(i=1;i<=n;i++)
        sum+=dis[i];
    printf("%d",sum);
}