最小生成樹(shù)有兩個(gè)經(jīng)典算法：Prim算法和Kruskal算法，Prim適合于點(diǎn)較少的圖，對(duì)于一個(gè)節(jié)點(diǎn)數(shù)為N的連通圖來(lái)說(shuō)，其時(shí)間復(fù)雜度為O(N^2)；Kruskal適合于邊較少的圖，對(duì)一個(gè)邊為E的連通圖來(lái)說(shuō)，其時(shí)間復(fù)雜度為O(ElogE)，因此要根據(jù)不同情況選擇合適的算法。
這里說(shuō)一下Prim算法。
Prim的具體步驟為把所有點(diǎn)分為兩個(gè)部分：屬于集合S，或不屬于S，當(dāng)所有點(diǎn)都屬于S時(shí)，算法結(jié)束。
1.初始條件先將第一個(gè)點(diǎn)p0劃到S中，然后利用p0關(guān)聯(lián)的所有邊更新cost[]（sost[i]表示pi與S中點(diǎn)相連的最短的那條邊長(zhǎng)）
2.每次從sost[]中選出最小的那一個(gè)cost[i]（i不能屬于S），將i加入到S中，并利用與i相關(guān)的邊更新cost[]（已加入到S中的點(diǎn)不用再更新）
3.反復(fù)執(zhí)行第二步，直到圖連通。（我們知道一個(gè)有n個(gè)節(jié)點(diǎn)的圖，最少只需要n-1條邊就可以連通了，所以第二步會(huì)執(zhí)行n-1次，每次都會(huì)在圖中加入一條邊）
關(guān)于Kruskal請(qǐng)參閱：
http://www.shnenglu.com/hoolee/archive/2012/08/04/186253.html下面是zoj1203的Prim算法代碼：
#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#define MAX 1000000
#define LENN 110
typedef struct
{
double x;
double y;
}Point;
int main()
{
int i, j, k;
int N, n;
Point ps[LENN];
double mp[LENN][LENN];
double cost[LENN];
int bl[LENN];
scanf("%d", &N);
n = 0;
int gard = 0;
while(N != 0)
{
for(i = 0; i < N; i++)
scanf("%lf%lf", &ps[i].x, &ps[i].y);
for(i = 0; i < N; i++)// make map[][]
for(j = i + 1; j < N; j++)
{
double dx = ps[i].x - ps[j].x;
double dy = ps[i].y - ps[j].y;
double lent = sqrt(dx * dx + dy * dy);
mp[i][j] = mp[j][i] = lent;
}
for(i = 0; i < N; i++)
mp[i][i] = MAX;
double len = 0;
bl[0] = 1;
for(i = 1; i < N; i++)
{
cost[i] = mp[0][i];
bl[i] = 0;
}
for(i = 0; i < N - 1; i++)//prim
{
double min = MAX;
int t;
for(k = 0; k < N; k++)
if(bl[k] == 0 && cost[k] < min)
{
min = cost[k];
t = k;
}
bl[t] = 1;
len += min;
for(j = 0; j < N; j++)// update
if(bl[j] == 0 && mp[t][j] < cost[j])
cost[j] = mp[t][j];
}
if(gard++ != 0)
putchar(10);
printf("Case #%d:\n", ++n);
printf("The minimal distance is: %.2lf\n", len);
scanf("%d", &N);
}
//system("pause");
}
posted on 2012-08-06 17:46
小鼠標(biāo) 閱讀(3142)
評(píng)論(0) 編輯 收藏 引用 所屬分類(lèi):
圖論