最小生成树

最小生成树

概念：
一个有 n 个结点的连通图的生成树是原图的极小连通子图，且包含原图中的所有 n 个结点，并且有保持图连通的最少的边。
简单来说：最小生成树就是在一个连通图（每个点都相连的无向图）中使得权值和最小的树，保证每个点都在里面。
最小生成树其实是最小权重生成树的简称。
最小生成树又叫“MST”。
应用：
例如：要在n个城市之间铺设光缆，主要目标是要使这 n 个城市的任意两个之间都可以通信，但铺设光缆的费用很高，且各个城市之间铺设光缆的费用不同，因此另一个目标是要使铺设光缆的总费用最低。这就需要找到带权的最小生成树。
求最小生成树的方法：
最小生成树可以用kruskal（克鲁斯卡尔）算法或prim（普里姆）算法求出。
prim算法的时间复杂度不依赖于排序算法，并且主要与点的个数有关，适用于密集图。
kruskal算法需要排序，但使用幷查集可节省判断时间，主要与边的条数有关，适用于稀疏图。
简单来说：Prim算法，适用于点少的图。Kruskal算法，适用于边少的图。

Prim算法

思想：
贪心思想：每次选取最小边。
算法描述：
（1）输入：一个加权连通图，其中顶点集合为V，边集合为E；
（2）初始化：Vnew = {x}，其中x为集合V中的任一节点（起始点），Enew = {},为空；
（3）重复下列操作，直到Vnew = V：
a.在集合E中选取权值最小的边< u, v >，其中u为集合Vnew中的元素，而v不在Vnew集合当中，并且v∈V（如果存在有多条满足前述条件即具有相同权值的边，则可任意选取其中之一）；
b.将v加入集合Vnew中，将< u, v >边加入集合Enew中；
（4）输出：使用集合Vnew和Enew来描述所得到的最小生成树。
图例描述：



时间复杂度：
顶点数v，边数e
邻接矩阵:O(v)
邻接表:O(elog2v)
代码：

#include<iostream>using namespace std;const int maxn=1000000;int n,k,tmp,ans,map[1001][1001],dis[maxn];bool flag[maxn];void prim(){    for(int i=1;i<=n;i++)    dis[i]=maxn;//初始化    dis[1]=0;    for(int i=1;i<=n;i++)    {        tmp=maxn;        for(int j=1;j<=n;j++)        if(!flag[j]&&tmp>dis[j])        {            tmp=dis[j];            k=j;        }//找出最小距离的节点         flag[k]=1;//把访问的节点做标记        for(int j=1;j<=n;j++)        if(!flag[j]&&dis[j]>map[k][j])        dis[j]=map[k][j];//更新最短距离    }}int main(){    cin>>n;    for(int i=1;i<=n;i++)//邻接矩阵储存权值       for(int j=1;j<=n;j++)      cin>>map[i][j];    prim();    for(int i=1;i<=n;i++)    ans+=dis[i];    cout<<ans;    return 0;}

Kruskal算法

基本思路：
克鲁斯卡尔算法是在剩下的所有未选取的边中，找最小边，如果和已选取的边构成回路，则放弃，选取次小边。

#include<iostream>#include<algorithm>using namespace std;int n,e,fa[101],sum=0;struct node{    int o;    int u;    int t;}a[101];int find(int x){    if(father[x]!=x)    father[x]=find(father[x])    return father[x];}int cmp(node a,node b)        {    return a.t < b.t;}void unionn(int x,int y){    int f1=find(x);    int f2=find(y);    if(f1!=f2)fa[f2]=f1;}int main(){    int k=0;    cin>>n>>e;    for(int i=1;i<=e;i++)    cin>>a[i].o>>a[i].u>>a[i].t;    for(int i=1;i<=n;i++)    fa[i]=i;    sort(a+1,a+1+e,cmp);    for(int i=1;i<=e;i++)    {        if(find(a[i].o)!=find(a[i].u))        {          unionn(a[i].o,a[i].u);          sum=sum+a[i].t;          k++;        }        if(k==n-1)break;    }    cout<<sum;}