还是畅通工程
Problem Description
某省调查乡村交通状况,得到的统计表中列出了任意两村庄间的距离。省政府“畅通工程”的目标是使全省任何两个村庄间都可以实现公路交通(但不一定有直接的公路相连,只要能间接通过公路可达即可),并要求铺设的公路总长度为最小。请计算最小的公路总长度。
Input
测试输入包含若干测试用例。每个测试用例的第1行给出村庄数目N ( < 100 );随后的N(N-1)/2行对应村庄间的距离,每行给出一对正整数,分别是两个村庄的编号,以及此两村庄间的距离。为简单起见,村庄从1到N编号。
当N为0时,输入结束,该用例不被处理。
Output
对每个测试用例,在1行里输出最小的公路总长度。
Sample Input
3
1 2 1
1 3 2
2 3 4
4
1 2 1
1 3 4
1 4 1
2 3 3
2 4 2
3 4 5
0
Sample Output
3
5
Hint
Huge input, scanf is recommended.
有n个城镇,给出了一些城镇之间的距离,要求把所有城镇连接起来的最短距离总和;
这题就是最小生成树的板子题。
最小生成树有两种算法可以解决(Prim算法和Kruskal算法)。
Kruskal算法就是把所有边拿出来比较,需要用到并查集,把所有边先进行从小到大排序,然后从最小的边开始判断,如果这条边连接的两个城镇在之前没有联系,这条边就要被选用,反之,则不选,一直选下去,最后得到的就是最短的距离和。
Prim算法强调的则是所有的点,先选一个点作为起点,然后把这个点连接的所有边都加到队列里,选权值最小的那条边,然后判断此时这条边连接的另一个点是否被访问过,如果没有,则继续把这个点的所有边也加入队列,一直这么判断。
Kruskal算法代码:
#include<cstdio>
#include<cstring>
#include<cmath>
#include<set>
#include<map>
#include<iostream>
#include<algorithm>
using namespace std;
#define ll long long
const int inf=0x3f3f3f3f;
const int maxn=5010;
#define eps 1e-8
struct node
{
int u,v,w;
}a[maxn];
int pre[maxn];
int n,m;
bool cmp(node x,node y)
{
return x.w<y.w;
}
int Find(int x)
{
if(x==pre[x]) return x;
return pre[x]=Find(pre[x]);
}
void init()
{
for(int i=0;i<=n;i++)
pre[i]=i;
}
void join(int x,int y)
{
int fx=Find(x);
int fy=Find(y);
if(fx!=fy)
pre[fx]=fy;
}
int main()
{
while(~scanf("%d",&n))
{
if(n==0) break;
init();
m=n*(n-1)/2;
for(int i=0;i<m;i++)
scanf("%d %d %d",&a[i].u,&a[i].v,&a[i].w);
sort(a,a+m,cmp);
int ans=0;
for(int i=0;i<m;i++)
{
if(Find(a[i].u)!=Find(a[i].v))
{
ans+=a[i].w;
join(a[i].u,a[i].v);
}
}
printf("%d\n",ans);
}
return 0;
}
Prim算法代码:
#include<cstdio>
#include<cstring>
#include<cmath>
#include<set>
#include<map>
#include<iostream>
#include<algorithm>
using namespace std;
#define ll long long
const int inf=0x3f3f3f3f;
const int maxn=110;
#define eps 1e-8
int n,m;
int mp[maxn][maxn];
int dis[maxn],vis[maxn];
void init()
{
for(int i=1; i<=n; i++)
for(int j=1; j<=n; j++)
if(i==j) mp[i][j]=0;
else mp[i][j]=inf;
}
int prim()
{
int sum=0;
for(int i=1; i<=n; i++)
{
vis[i]=0;
dis[i]=mp[1][i];
}
for(int i=1; i<=n; i++)
{
int minn=inf,v;
for(int j=1; j<=n; j++)
{
if(!vis[j]&&dis[j]<minn)
{
minn=dis[j];
v=j;
}
}
vis[v]=1;
sum+=dis[v];
for(int j=1; j<=n; j++)
{
if(!vis[j]&&mp[v][j]<dis[j])
{
dis[j]=mp[v][j];
}
}
}
return sum;
}
int main()
{
while(~scanf("%d",&n))
{
init();
if(n==0) break;
m=n*(n-1)/2;
int u,v,w;
for(int i=0; i<m; i++)
{
scanf("%d %d %d",&u,&v,&w);
mp[u][v]=min(mp[u][v],w);
mp[v][u]=mp[u][v];
}
printf("%d\n",prim());
}
return 0;
}