HDU - 1863 畅通工程 (最小生成树–Kruskal&prime算法)
省政府“畅通工程”的目标是使全省任何两个村庄间都可以实现公路交通(但不一定有直接的公路相连,只要能间接通过公路可达即可)。经过调查评估,得到的统计表中列出了有可能建设公路的若干条道路的成本。现请你编写程序,计算出全省畅通需要的最低成本。
Input
测试输入包含若干测试用例。每个测试用例的第1行给出评估的道路条数 N、村庄数目M ( < 100 );随后的 N
行对应村庄间道路的成本,每行给出一对正整数,分别是两个村庄的编号,以及此两村庄间道路的成本(也是正整数)。为简单起见,村庄从1到M编号。当N为0时,全部输入结束,相应的结果不要输出。
Output
对每个测试用例,在1行里输出全省畅通需要的最低成本。若统计数据不足以保证畅通,则输出“?”。
Sample Input
3 3
1 2 1
1 3 2
2 3 4
1 3
2 3 2
0 100
Sample Output
3
?
这是一个裸的最小生成树
借助这个题来简单介绍一下求最小生成树的两种方法
- Kruskal算法
又叫加边法,每一次选一条权值最小的边,看是否成环,如果成环就不加这条边。如果能加到n-1条边就说明找到了最小生成树(n是顶点的个数)。因为一个树的边数就是等于顶点个数减一。
具体实现细节:
1.用一个结构体来存边,里面存的是一条边的两个顶点和这条边的权值。
struct node{
int s;
int e;
int t;
}side[maxn];
2.把所有的边按权值从小到大排序。然后小到大开始一条条的判断能不能加入。
3.加边的话用并查集来操作,判断加不加的话,就看这条边的两个父节点是否是一个,如果是同一父节点。那就形成环了,如果不同,那就把merge()这两个点。
if(find(side[i].e)!=find(side[i].s))
{
merge(side[i].e,side[i].s);
ans+=side[i].t;
cnt++;
}
完整代码如下:
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int maxn=100010;
int f[1010];
struct node{
int s;
int e;
int t;
}side[maxn];
int n,m;
int ans=0;
void init()
{
for(int i=0;i<=n;i++)
f[i]=i;
ans=0;
return ;
}
int find(int x)
{
if(x==f[x])
return x;
else
f[x]=find(f[x]);
return f[x];
}
void merge(int x,int y)
{
int fx=find(x);
int fy=find(y);
if(fx!=fy)
{
if(fx>fy)
f[fx]=fy;
else
f[fy]=fx;
}
return ;
}
bool cmp(struct node a,struct node b)
{
return a.t<b.t;
}
int main()
{
while(cin>>m>>n)
{
if(m==0)
break;
for(int i=1;i<=m;i++)
{
int a,b,c;
cin>>a>>b>>c;
side[i].s=a;
side[i].e=b;
side[i].t=c;
}
init();
sort(side+1,side+1+m,cmp);
int cnt=0;
int flag=0;
for(int i=1;i<=m;i++)
{
if(find(side[i].e)!=find(side[i].s))
{
merge(side[i].e,side[i].s);
ans+=side[i].t;
cnt++;
}
if(cnt==n-1)
{
flag=1;
break;
}
}
if(!flag)
cout<<"?"<<endl;
else
cout<<ans<<endl;
}
return 0;
}
- prime算法
Prim算法通常以邻接矩阵作为储存结构。
它的基本思想是以顶点为主导地位,从起始顶点出发,通过选择当前可用的最小权值边把顶点加入到生成树当中来:
1.从连通网络N={V,E}中的某一顶点U0出发,选择与它关联的具有最小权值的边(U0,V),将其顶点加入到生成树的顶点集合U中。
2.以后每一步从一个顶点在U中,而另一个顶点不在U中的各条边中选择权值最小的边(U,V),把它的顶点加入到集合U中。如此继续下去,直到网络中的所有顶点都加入到生成树顶点集合U中为止。
引自https://blog.csdn.net/Feynman1999/article/details/73385463
代码如下:
#include <iostream>
#include <cstring>
#include <cstdio>
#include <cmath>
#define inf 0x3f3f3f3f
using namespace std;
const int maxx=110;
int dis[maxx];
int mp[maxx][maxx];
int vis[maxx];
int ans=0;
int n,m;
void init()
{
memset(vis,0,sizeof(vis));
memset(mp,inf,sizeof(mp));
ans=0;
return ;
}
int main()
{
while(cin>>m>>n)
{
if(m==0)
break;
init();
for(int i=0;i<m;i++)
{
int a,b,c;
cin>>a>>b>>c;
mp[a][b]=mp[b][a]=c;
}
vis[1]=1;
dis[1]=0;
int flag=1;
for(int i=1;i<=n;i++)
dis[i]=mp[1][i];
for(int i=1;i<n;i++)///循环找n-1次 找n-1条边 能找到说明找到了最小生成树
{
int minn=inf;
int minid=inf;
for(int j=1;j<=n;j++)
{
if(vis[j]==0&&dis[j]<minn)
{
minn=dis[j];
minid=j;
}
}
//cout<<minn<<endl;
if(minid==inf)
{
flag=0;
break;
}
vis[minid]=1;
ans+=minn;
///更新最小距离
for(int j=1;j<=n;j++)
{
if(vis[j]==1)
continue;
dis[j]=min(dis[j],mp[minid][j]);
}
}
if(flag==0)
cout<<"?"<<endl;
else
cout<<ans<<endl;
}
return 0;
}