基本介紹
- 中文名:普里姆算法
- 外文名:Prim Algorithm
- 別稱:最小生成樹算法
- 提出者:沃伊捷赫·亞爾尼克(Vojtěch Jarník)
- 提出時間:1930年
- 套用學科:計算機,數據結構,數學(圖論)
- 適用領域範圍:套用圖論知識的實際問題
- 算法:貪心
算法描述
時間複雜度
最小邊、權的數據結構 | 時間複雜度(總計) |
---|---|
鄰接矩陣、搜尋 | O(V^2) |
O((V + E) log(V)) = O(E log(V)) | |
O(E + V log(V)) |
圖例描述
圖例 | 說明 | 不可選 | 可選 | 已選(Vnew) |
此為原始的加權連通圖。每條邊一側的數字代表其權值。 | - | - | - | |
頂點D被任意選為起始點。頂點A、B、E和F通過單條邊與D相連。A是距離D最近的頂點,因此將A及對應邊AD以高亮表示。 | C, G | A, B, E, F | D | |
下一個頂點為距離D或A最近的頂點。B距D為9,距A為7,E為15,F為6。因此,F距D或A最近,因此將頂點F與相應邊DF以高亮表示。 | C, G | B, E, F | A, D | |
算法繼續重複上面的步驟。距離A為7的頂點B被高亮表示。 | C | B, E, G | A, D, F | |
在當前情況下,可以在C、E與G間進行選擇。C距B為8,E距B為7,G距F為11。點E最近,因此將頂點E與相應邊BE高亮表示。 | 無 | C, E, G | A, D, F, B | |
這裡,可供選擇的頂點只有C和G。C距E為5,G距E為9,故選取C,並與邊EC一同高亮表示。 | 無 | C, G | A, D, F, B, E | |
頂點G是唯一剩下的頂點,它距F為11,距E為9,E最近,故高亮表示G及相應邊EG。 | 無 | G | A, D, F, B, E, C | |
現在,所有頂點均已被選取,圖中綠色部分即為連通圖的最小生成樹。在此例中,最小生成樹的權值之和為39。 | 無 | 無 | A, D, F, B, E, C, G |
代碼
PASCAL代碼
programprim;constmap:array[1..6,1..6]ofinteger=((0,6,1,5,0,0),(6,0,5,0,3,0),(1,5,0,5,6,4),(5,0,5,0,0,2),(0,3,6,0,0,6),(0,0,4,2,6,0));//樣例輸入vari,j,l:integer;min,minn:longint;f,d:array[1..6]ofinteger;s:array[1..6,1..3]ofinteger;v,p:setof1..6;beginl:=1;p:=[];v:=[];fori:=2to6dov:=v+[i];p:=p+[1];fori:=1to6dof[i]:=1000;//還可以寫成filldword(f,sizeof(f)div2,1000)f[1]:=0;fori:=1to6dod[i]:=0;//還可以寫成fillchar(d,sizeof(d),0);s[1,3]:=0;fori:=1to5dobeginmin:=1000;forj:=1to6dobeginif(f[j]>map[l,j])and(jinv)and(map[l,j]<>0)thenbeginf[j]:=map[l,j];d[j]:=l;end;if(f[j]<min)and(f[j]>0)thenbeginmin:=f[j];minn:=j;end;writeln(d[j]);end;f[minn]:=0;v:=v-[minn];p:=p+[minn];s[i,1]:=d[minn];l:=minn;s[i,2]:=minn;s[i,3]:=min;end;fori:=1to5dowrite(s[i,1],'to',s[i,2],'=',s[i,3],'-->');readln;end.
c代碼
#include<stdio.h>#include<stdlib.h>#definemax1000000000;inta[1001][1001],d[1001],p[1001];intmain(){inti,j,k,m,n,min,ans,t;intx,y,z;scanf("%d%d",&n,&m);for(i=1;i<=m;i++){scanf("%d%d%d",&x,&y,&z);a[x][y]=z;a[y][x]=z;}for(i=1;i<=n;i++)d[i]=1000000000;d[1]=0;for(i=2;i<=n;i++){min=max;for(j=1;j<=n;j++)if(!p[j]&&min>d[j])min=d[j];t=j;}p[t]=j;for(j=1;j<=n;j++)if(a[t][j]=0&&d[j]>a[t][j]){d[j]=a[t][j];ans+=min;}printf("%d",ans);return0;}
C++代碼
#defineMAXN1000#defineINF1<<30intclosest[MAXN],lowcost[MAXN],m;//m為節點的個數intG[MAXN][MAXN];//鄰接矩陣intprim(){for(inti=0;i<m;i++){lowcost[i]=INF;}for(inti=0;i<m;i++){closest[i]=0;}closest[0]=-1;//加入第一個點,-1表示該點在集合U中,否則在集合V中intnum=0,ans=0,e=0;//e為最新加入集合的點while(num<m-1)//加入m-1條邊{intmicost=INF,miedge=-1;for(inti=0;i<m;i++)if(closest[i]!=-1){inttemp=G[e][i];if(temp<lowcost[i]){lowcost[i]=temp;closest[i]=e;}if(lowcost[i]<micost)micost=lowcost[miedge=i];}ans+=micost;closest[e=miedge]=-1;num++;}returnans;}
時間複雜度
鄰接矩陣:O(v) 鄰接表:O(elog2v)