最小生成樹

用來解決工程中的代價問題。

?

一：普里姆算法

?

?

具體代碼用C語言實(shí)現(xiàn)如下：

typedef int VRType;typedef char InfoType;#define MAX_NAME 3 /* 頂點(diǎn)字符串的最大長度+1 */#define MAX_INFO 20 /* 相關(guān)信息字符串的最大長度+1 */typedef char VertexType[MAX_NAME]; #include<string.h>#include<ctype.h>#include<malloc.h> /* malloc()等 */#include<limits.h> /* INT_MAX等 */#include<stdio.h> /* EOF(=^Z或F6),NULL */#include<stdlib.h> /* atoi() */#include<io.h> /* eof() */#include<math.h> /* floor(),ceil(),abs() */#include<process.h> /* exit() *//* 函數(shù)結(jié)果狀態(tài)代碼 */#define TRUE 1#define FALSE 0#define OK 1#define ERROR 0#define INFEASIBLE -1typedef int Status; /* Status是函數(shù)的類型,其值是函數(shù)結(jié)果狀態(tài)代碼，如OK等 */typedef int Boolean; /* Boolean是布爾類型,其值是TRUE或FALSE */#define MAX_NAME 5 /* 頂點(diǎn)字符串的最大長度 */typedef int InfoType;typedef char VertexType[MAX_NAME]; /* 字符串類型 *//* 圖的鄰接表存儲表示 */#define MAX_VERTEX_NUM 20typedef enum{DG,DN,AG,AN}GraphKind; /* {有向圖,有向網(wǎng),無向圖,無向網(wǎng)} */typedef struct ArcNode{int adjvex; /* 該弧所指向的頂點(diǎn)的位置 */struct ArcNode *nextarc; /* 指向下一條弧的指針 */InfoType *info; /* 網(wǎng)的權(quán)值指針） */}ArcNode; /* 表結(jié)點(diǎn) */typedef struct{VertexType data; /* 頂點(diǎn)信息 */ArcNode *firstarc; /* 第一個表結(jié)點(diǎn)的地址,指向第一條依附該頂點(diǎn)的弧的指針 */}VNode,AdjList[MAX_VERTEX_NUM]; /* 頭結(jié)點(diǎn) */typedef struct{AdjList vertices;int vexnum,arcnum; /* 圖的當(dāng)前頂點(diǎn)數(shù)和弧數(shù) */int kind; /* 圖的種類標(biāo)志 */}ALGraph;/* 圖的鄰接表存儲的基本操作(15個) */int LocateVex(ALGraph G,VertexType u){ /* 初始條件: 圖G存在,u和G中頂點(diǎn)有相同特征 *//* 操作結(jié)果: 若G中存在頂點(diǎn)u,則返回該頂點(diǎn)在圖中位置;否則返回-1 */int i;for(i=0;i<G.vexnum;++i)if(strcmp(u,G.vertices[i].data)==0) //比較頂點(diǎn)是否相等return i;return -1;}Status CreateGraph(ALGraph *G){ /* 采用鄰接表存儲結(jié)構(gòu),構(gòu)造沒有相關(guān)信息的圖G(用一個函數(shù)構(gòu)造4種圖) */int i,j,k;int w; /* 權(quán)值 */VertexType va,vb;ArcNode *p;printf("請輸入圖的類型(有向圖:0,有向網(wǎng):1,無向圖:2,無向網(wǎng):3): ");scanf("%d",&(*G).kind);printf("請輸入圖的頂點(diǎn)數(shù),邊數(shù): ");scanf("%d,%d",&(*G).vexnum,&(*G).arcnum);printf("請輸入%d個頂點(diǎn)的值(<%d個字符):\n",(*G).vexnum,MAX_NAME);for(i=0;i<(*G).vexnum;++i) /* 構(gòu)造頂點(diǎn)向量 */{scanf("%s",(*G).vertices[i].data);(*G).vertices[i].firstarc=NULL;}if((*G).kind==1||(*G).kind==3) /* 網(wǎng) */printf("請順序輸入每條弧(邊)的權(quán)值、弧尾和弧頭(以空格作為間隔):\n");else /* 圖 */printf("請順序輸入每條弧(邊)的弧尾和弧頭(以空格作為間隔):\n");for(k=0;k<(*G).arcnum;++k) /* 構(gòu)造表結(jié)點(diǎn)鏈表 */{if((*G).kind==1||(*G).kind==3) /* 網(wǎng) */scanf("%d%s%s",&w,va,vb);else /* 圖 */scanf("%s%s",va,vb);i=LocateVex(*G,va); /* 弧尾 */j=LocateVex(*G,vb); /* 弧頭 */p=(ArcNode*)malloc(sizeof(ArcNode));p->adjvex=j;if((*G).kind==1||(*G).kind==3) /* 網(wǎng) */{p->info=(int *)malloc(sizeof(int));*(p->info)=w; //網(wǎng)的權(quán)值指針}elsep->info=NULL; /* 圖 */p->nextarc=(*G).vertices[i].firstarc; /* 插在表頭 */(*G).vertices[i].firstarc=p;if((*G).kind>=2) /* 無向圖或網(wǎng),產(chǎn)生第二個表結(jié)點(diǎn) */{p=(ArcNode*)malloc(sizeof(ArcNode));p->adjvex=i;if((*G).kind==3) /* 無向網(wǎng) */{p->info=(int*)malloc(sizeof(int));*(p->info)=w;}elsep->info=NULL; /* 無向圖 */p->nextarc=(*G).vertices[j].firstarc; /* 插在表頭 */(*G).vertices[j].firstarc=p;}}return OK;}void DestroyGraph(ALGraph *G){ /* 初始條件: 圖G存在。操作結(jié)果: 銷毀圖G */int i;ArcNode *p,*q;(*G).vexnum=0;(*G).arcnum=0;for(i=0;i<(*G).vexnum;++i){p=(*G).vertices[i].firstarc; //第一個表結(jié)點(diǎn)的地址，頭結(jié)點(diǎn)while(p){q=p->nextarc;if((*G).kind%2) /* 網(wǎng) */free(p->info);free(p);p=q;}}}VertexType* GetVex(ALGraph G,int v){ /* 初始條件: 圖G存在,v是G中某個頂點(diǎn)的序號。操作結(jié)果: 返回v的值 */if(v>=G.vexnum||v<0)exit(ERROR);return &G.vertices[v].data;}Status PutVex(ALGraph *G,VertexType v,VertexType value){ /* 初始條件: 圖G存在,v是G中某個頂點(diǎn) *//* 操作結(jié)果: 對v賦新值value */int i;i=LocateVex(*G,v); //先找到該頂點(diǎn)if(i>-1) /* v是G的頂點(diǎn) */{strcpy((*G).vertices[i].data,value);return OK;}return ERROR;}int FirstAdjVex(ALGraph G,VertexType v){ /* 初始條件: 圖G存在,v是G中某個頂點(diǎn) *//* 操作結(jié)果: 返回v的第一個鄰接頂點(diǎn)的序號。若頂點(diǎn)在G中沒有鄰接頂點(diǎn),則返回-1 */ArcNode *p;int v1;v1=LocateVex(G,v); /* v1為頂點(diǎn)v在圖G中的序號 */p=G.vertices[v1].firstarc;if(p)return p->adjvex;elsereturn -1;}int NextAdjVex(ALGraph G,VertexType v,VertexType w){ /* 初始條件: 圖G存在,v是G中某個頂點(diǎn),w是v的鄰接頂點(diǎn) *//* 操作結(jié)果: 返回v的(相對于w的)下一個鄰接頂點(diǎn)的序號。 *//* 若w是v的最后一個鄰接點(diǎn),則返回-1 */ArcNode *p;int v1,w1;v1=LocateVex(G,v); /* v1為頂點(diǎn)v在圖G中的序號 */w1=LocateVex(G,w); /* w1為頂點(diǎn)w在圖G中的序號 */p=G.vertices[v1].firstarc;while(p&&p->adjvex!=w1) /* 指針p不空且所指表結(jié)點(diǎn)不是w */p=p->nextarc;if(!p||!p->nextarc) /* 沒找到w或w是最后一個鄰接點(diǎn) */return -1;else /* p->adjvex==w */return p->nextarc->adjvex; /* 返回v的(相對于w的)下一個鄰接頂點(diǎn)的序號 */}void InsertVex(ALGraph *G,VertexType v){ /* 初始條件: 圖G存在,v和圖中頂點(diǎn)有相同特征 *//* 操作結(jié)果: 在圖G中增添新頂點(diǎn)v(不增添與頂點(diǎn)相關(guān)的弧,留待InsertArc()去做) */strcpy((*G).vertices[(*G).vexnum].data,v); /* 構(gòu)造新頂點(diǎn)向量 */(*G).vertices[(*G).vexnum].firstarc=NULL;(*G).vexnum++; /* 圖G的頂點(diǎn)數(shù)加1 */}Status DeleteVex(ALGraph *G,VertexType v){ /* 初始條件: 圖G存在,v是G中某個頂點(diǎn) *//* 操作結(jié)果: 刪除G中頂點(diǎn)v及其相關(guān)的弧 */int i,j;ArcNode *p,*q;j=LocateVex(*G,v); /* j是頂點(diǎn)v的序號 */if(j<0) /* v不是圖G的頂點(diǎn) */return ERROR;p=(*G).vertices[j].firstarc; /* 刪除以v為出度的弧或邊 */while(p){q=p;p=p->nextarc;if((*G).kind%2) /* 網(wǎng) */free(q->info);free(q);(*G).arcnum--; /* 弧或邊數(shù)減1 */}(*G).vexnum--; /* 頂點(diǎn)數(shù)減1 */for(i=j;i<(*G).vexnum;i++) /* 頂點(diǎn)v后面的頂點(diǎn)前移 */(*G).vertices[i]=(*G).vertices[i+1];for(i=0;i<(*G).vexnum;i++) /* 刪除以v為入度的弧或邊且必要時修改表結(jié)點(diǎn)的頂點(diǎn)位置值 */{p=(*G).vertices[i].firstarc; /* 指向第1條弧或邊 */while(p) /* 有弧 */{if(p->adjvex==j){if(p==(*G).vertices[i].firstarc) /* 待刪結(jié)點(diǎn)是第1個結(jié)點(diǎn) */{(*G).vertices[i].firstarc=p->nextarc;if((*G).kind%2) /* 網(wǎng) */free(p->info);free(p);p=(*G).vertices[i].firstarc;if((*G).kind<2) /* 有向 */(*G).arcnum--; /* 弧或邊數(shù)減1 */}else{q->nextarc=p->nextarc;if((*G).kind%2) /* 網(wǎng) */free(p->info);free(p);p=q->nextarc;if((*G).kind<2) /* 有向 */(*G).arcnum--; /* 弧或邊數(shù)減1 */}}else{if(p->adjvex>j)p->adjvex--; /* 修改表結(jié)點(diǎn)的頂點(diǎn)位置值(序號) */q=p;p=p->nextarc;}}}return OK;}Status InsertArc(ALGraph *G,VertexType v,VertexType w){ /* 初始條件: 圖G存在,v和w是G中兩個頂點(diǎn) *//* 操作結(jié)果: 在G中增添弧<v,w>,若G是無向的,則還增添對稱弧<w,v> */ArcNode *p;int w1,i,j;i=LocateVex(*G,v); /* 弧尾或邊的序號 */j=LocateVex(*G,w); /* 弧頭或邊的序號 */if(i<0||j<0)return ERROR;(*G).arcnum++; /* 圖G的弧或邊的數(shù)目加1 */if((*G).kind%2) /* 網(wǎng) */{printf("請輸入弧(邊)%s→%s的權(quán)值: ",v,w);scanf("%d",&w1);}p=(ArcNode*)malloc(sizeof(ArcNode));p->adjvex=j;if((*G).kind%2) /* 網(wǎng) */{p->info=(int*)malloc(sizeof(int));*(p->info)=w1;}elsep->info=NULL;p->nextarc=(*G).vertices[i].firstarc; /* 插在表頭 */(*G).vertices[i].firstarc=p;if((*G).kind>=2) /* 無向,生成另一個表結(jié)點(diǎn) */{p=(ArcNode*)malloc(sizeof(ArcNode));p->adjvex=i;if((*G).kind==3) /* 無向網(wǎng) */{p->info=(int*)malloc(sizeof(int));*(p->info)=w1;}elsep->info=NULL;p->nextarc=(*G).vertices[j].firstarc; /* 插在表頭 */(*G).vertices[j].firstarc=p;}return OK;}Status DeleteArc(ALGraph *G,VertexType v,VertexType w){ /* 初始條件: 圖G存在,v和w是G中兩個頂點(diǎn) *//* 操作結(jié)果: 在G中刪除弧<v,w>,若G是無向的,則還刪除對稱弧<w,v> */ArcNode *p,*q;int i,j;i=LocateVex(*G,v); /* i是頂點(diǎn)v(弧尾)的序號 */j=LocateVex(*G,w); /* j是頂點(diǎn)w(弧頭)的序號 */if(i<0||j<0||i==j)return ERROR;p=(*G).vertices[i].firstarc; /* p指向頂點(diǎn)v的第一條出弧 */while(p&&p->adjvex!=j) /* p不空且所指之弧不是待刪除弧<v,w> */{ /* p指向下一條弧 */q=p;p=p->nextarc;}if(p&&p->adjvex==j) /* 找到弧<v,w> */{if(p==(*G).vertices[i].firstarc) /* p所指是第1條弧 */(*G).vertices[i].firstarc=p->nextarc; /* 指向下一條弧 */elseq->nextarc=p->nextarc; /* 指向下一條弧 */if((*G).kind%2) /* 網(wǎng) */free(p->info);free(p); /* 釋放此結(jié)點(diǎn) */(*G).arcnum--; /* 弧或邊數(shù)減1 */}if((*G).kind>=2) /* 無向,刪除對稱弧<w,v> */{p=(*G).vertices[j].firstarc; /* p指隙サ鉾的第一條出弧 */while(p&&p->adjvex!=i) /* p不空且所指之弧不是待刪除弧<w,v> */{ /* p指向下一條弧 */q=p;p=p->nextarc;}if(p&&p->adjvex==i) /* 找到弧<w,v> */{if(p==(*G).vertices[j].firstarc) /* p所指是第1條弧 */(*G).vertices[j].firstarc=p->nextarc; /* 指向下一條弧 */elseq->nextarc=p->nextarc; /* 指向下一條弧 */if((*G).kind==3) /* 無向網(wǎng) */free(p->info);free(p); /* 釋放此結(jié)點(diǎn) */}}return OK;}typedef struct{ /* 記錄從頂點(diǎn)集U到V-U的代價最小的邊的輔助數(shù)組定義 */VertexType adjvex;VRType lowcost;}minside[MAX_VERTEX_NUM];int minimum(minside SZ,MGraph G){ /* 求closedge.lowcost的最小正值 */int i=0,j,k,min;while(!SZ[i].lowcost)i++;min=SZ[i].lowcost; /* 第一個不為0的值 */k=i;for(j=i+1;j<G.vexnum;j++)if(SZ[j].lowcost>0)if(min>SZ[j].lowcost){min=SZ[j].lowcost;k=j;}return k;}void MiniSpanTree_PRIM(MGraph G,VertexType u){ /* 用普里姆算法從第u個頂點(diǎn)出發(fā)構(gòu)造網(wǎng)G的最小生成樹T,輸出T的各條邊 算法7.9 */int i,j,k;minside closedge;k=LocateVex(G,u);for(j=0;j<G.vexnum;++j) /* 輔助數(shù)組初始化 */{if(j!=k){strcpy(closedge[j].adjvex,u);closedge[j].lowcost=G.arcs[k][j].adj;}}closedge[k].lowcost=0; /* 初始,U={u} */printf("最小代價生成樹的各條邊為:\n");for(i=1;i<G.vexnum;++i){ /* 選擇其余G.vexnum-1個頂點(diǎn) */k=minimum(closedge,G); /* 求出T的下一個結(jié)點(diǎn)：第K頂點(diǎn) */printf("(%s-%s)\n",closedge[k].adjvex,G.vexs[k]); /* 輸出生成樹的邊 */closedge[k].lowcost=0; /* 第K頂點(diǎn)并入U集 */for(j=0;j<G.vexnum;++j)if(G.arcs[k][j].adj<closedge[j].lowcost){ /* 新頂點(diǎn)并入U集后重新選擇最小邊 */strcpy(closedge[j].adjvex,G.vexs[k]);closedge[j].lowcost=G.arcs[k][j].adj;}}}void main(){MGraph G;CreateAN(&G);MiniSpanTree_PRIM(G,G.vexs[0]);}

注意：普里姆算法的時間復(fù)雜度為O()，與網(wǎng)中的邊數(shù)無關(guān)，因此適用于求邊稠密的網(wǎng)的最小生成樹。

?

?

二：克魯斯卡爾算法

?

分析過程：在網(wǎng)中先選取最小邊，繼續(xù)找下一條次小邊，如果和已選的邊構(gòu)成回路則舍棄，繼續(xù)尋找直到所有頂點(diǎn)。

?

克魯斯卡爾算法的時間復(fù)雜度為O(e*)。
.

創(chuàng)作挑戰(zhàn)賽新人創(chuàng)作獎勵來咯，堅持創(chuàng)作打卡瓜分現(xiàn)金大獎

總結(jié)

以上是生活随笔為你收集整理的最小生成树——普里姆算法和克鲁斯卡尔算法的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。