最小生成树

MST（Minimum Spanning Tree，最小生成树）问题有两种通用的解法，Prim算法就是其中之一，它是从点的方面考虑构建一颗MST，大致思想是：设图G顶点集合为U，首先任意选择图G中的一点作为起始点a，将该点加入集合V，再从集合U-V中找到另一点b使得点b到V中任意一点的权值最小，此时将b点也加入集合V；以此类推，现在的集合V={a，b}，再从集合U-V中找到另一点c使得点c到V中任意一点的权值最小，此时将c点加入集合V，直至所有顶点全部被加入V，此时就构建出了一颗MST。因为有N个顶点，所以该MST就有N-1条边，每一次向集合V中加入一个点，就意味着找到一条MST的边。

用图示和代码说明：

初始状态：

设置2个数据结构：

lowcost[i]:表示以i为终点的边的最小权值,当lowcost[i]=0说明以i为终点的边的最小权值=0,也就是表示i点加入了MST

mst[i]:表示对应lowcost[i]的起点，即说明边<mst[i],i>是MST的一条边，当mst[i]=0表示起点i加入MST   //起点存索引方便一些

我们假设V1是起始点，进行初始化（*代表无限大，即无通路）：

lowcost[2]=6，lowcost[3]=1，lowcost[4]=5，lowcost[5]=*，lowcost[6]=*

mst[2]=1，mst[3]=1，mst[4]=1，mst[5]=1，mst[6]=1，（所有点默认起点是V1）

明显看出，以V3为终点的边的权值最小=1，所以边<mst[3],3>=1加入MST

此时，因为点V3的加入，需要更新lowcost数组和mst数组：

lowcost[2]=5，lowcost[3]=0，lowcost[4]=5，lowcost[5]=6，lowcost[6]=4

mst[2]=3，mst[3]=0，mst[4]=1，mst[5]=3，mst[6]=3

明显看出，以V6为终点的边的权值最小=4，所以边<mst[6],6>=4加入MST

此时，因为点V6的加入，需要更新lowcost数组和mst数组：

lowcost[2]=5，lowcost[3]=0，lowcost[4]=2，lowcost[5]=6，lowcost[6]=0

mst[2]=3，mst[3]=0，mst[4]=6，mst[5]=3，mst[6]=0

明显看出，以V4为终点的边的权值最小=2，所以边<mst[4],4>=4加入MST

此时，因为点V4的加入，需要更新lowcost数组和mst数组：

lowcost[2]=5，lowcost[3]=0，lowcost[4]=0，lowcost[5]=6，lowcost[6]=0

mst[2]=3，mst[3]=0，mst[4]=0，mst[5]=3，mst[6]=0

明显看出，以V2为终点的边的权值最小=5，所以边<mst[2],2>=5加入MST

此时，因为点V2的加入，需要更新lowcost数组和mst数组：

lowcost[2]=0，lowcost[3]=0，lowcost[4]=0，lowcost[5]=3，lowcost[6]=0

mst[2]=0，mst[3]=0，mst[4]=0，mst[5]=2，mst[6]=0

很明显，以V5为终点的边的权值最小=3，所以边<mst[5],5>=3加入MST

lowcost[2]=0，lowcost[3]=0，lowcost[4]=0，lowcost[5]=0，lowcost[6]=0

mst[2]=0，mst[3]=0，mst[4]=0，mst[5]=0，mst[6]=0

至此，MST构建成功，如图所示：

根据以上得到以下代码(望各位看下,发表一下意见，有什么不足之处，还望指出)

#include <iostream>
using namespace std;
#define VerMax 100
typedef char VerType;
typedef int ArcType;
typedef struct{
VerType Vertex[VerMax];
ArcType arc[VerMax][VerMax];
int vexnum,arcnum;
}MGraph;
 struct Edge//辅助结构体,记录u--v-u的最小代价边的情况
{
 VerType vex;//顶点名称 ,起点
 ArcType lowcost;//以j为终点的最小边,这两个确定具体边,贪心算法;
};
Edge e[VerMax];
int min(MGraph *G)//求出生成树的下一个顶点
{
int i,min=INT_MAX,index=0;
    for (i=1;i<=G->vexnum;i++)
if (e[i].lowcost<min&&e[i].lowcost!=0)
{
index=i;
            min=e[i].lowcost;
}
return index;
}
void prim(MGraph *G,VerType c)
{
  int k,i,j;
  for (i=1;i<=G->vexnum;i++)
 if (G->Vertex[i]==c)
 {
 k=i;
 //将每个元素都以这个元素为起点
           for (j=1;j<=G->vexnum;j++)
  e[j].vex=c;
 }
for (i=1;i<=G->vexnum;i++)
if (i!=k)
e[i].lowcost=G->arc[k][i];
e[k].lowcost=0;
     for (i=1;i<G->vexnum;i++)
{
                k=min(G);//我每个最小边都保存了起点和终点;
cout<<"("<<e[k].vex<<","<< G->Vertex[k]<<")"<<endl;
e[k].lowcost=0;
for (j=1;j<=G->vexnum;j++)
{
if (G->arc[k][j]<e[j].lowcost&&e[j].lowcost!=0)//这里就将u和u-v的以j为终点的边做了比较
{
e[j].vex=G->Vertex[k];//以k为起点;
e[j].lowcost=G->arc[k][j];//在添加了一个顶点进u中,更新以j为终点的边.
}
}
}
     
 


}
int main()
{
MGraph G;
cin>>G.vexnum>>G.arcnum;
//memset(G.arc,INT_MAX,sizeof(G.arcnum));//这个没赋成功，这能赋字符类型的,它是按一个哟个字节赋的;
//memset(visited,0,sizeof(visited));
int i,j,v;
for (i=0;i<=G.vexnum;i++)//初始化边,求小赋大
for (j=0;j<=G.vexnum;j++)
G.arc[i][j]=INT_MAX;


for (i=0;i<VerMax;i++)
    e[i].lowcost=INT_MAX;


for (j=1;j<=G.vexnum;j++)
{
cin>>e[j].vex;//初始化顶点
G.Vertex[j]=e[j].vex;
e[j].lowcost=INT_MAX;
}
while (G.arcnum--)
{
cin>>i>>j>>v;//初始化边；
G.arc[i][j]=v;
G.arc[j][i]=v;
}
char v1;
cin>>v1;
    prim(&G,v1);
return 0;

}