还是畅通工程
某省调查乡村交通状况,得到的统计表中列出了任意两村庄间的距离。省政府“畅通工程”的目标是使全省任何两个村庄间都可以实现公路交通(但不一定有直接的公路相连,只要能间接通过公路可达即可),并要求铺设的公路总长度为最小。请计算最小的公路总长度。
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
Huge input, scanf is recommended.
思路:最小生成树模板题
Kruskal算法,
1.把图中的所有边按权值从小到大排序。 2.把图中的n个顶点看成独立的n棵独立的树。 3.所选的边连接的两个顶点a,b,应属于两颗不同的树,则成为最小生成树的一条边,并将这两颗树合并作为一颗树。 4. 重复,直到所有顶点都在一颗树内或者树上有n-1条边为止。
代码:
#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
struct node
{
int a,b,s;
}p[10005];
int f[105];
bool cmp(node x,node y) //按权值对边排序
{
return x.s<y.s;
}
int find(int x)
{
if(f[x]==x)
return x;
return f[x]=find(f[x]);
}
int join(int x,int y) //判断两点是否在一棵树上
{
int t1,t2;
t1=find(x);
t2=find(y);
if(t1!=t2)//如果不是就把边加入树中并把两个点连到一棵树上
{
f[t2]=t1;
return 1;
}
return 0;
}
int main()
{
int n,m;
while(~scanf("%d",&n)&&n)
{
memset(f,0,sizeof(0));
int sum=0,count=0;
m=(n-1)*n/2;
for(int i=0;i<m;i++)
scanf("%d%d%d",&p[i].a,&p[i].b,&p[i].s);
sort(p,p+m,cmp);
for(int i=1;i<=n;i++)
f[i]=i;
for(int i=0;i<m;i++)
{
if(join(p[i].a,p[i].b))
{
count++;
sum+=p[i].s;
}
if(count==n-1) //n个点之间最多n-1条线
break;
}
printf("%d\n",sum);
}
return 0;
}
Prim算法:
“加点法”,每次迭代选择代价最小的边对应的点,加入到最小生成树中。从某一个顶点s开始,逐渐长大覆盖整个连通网的所有顶点。
1.图的所有顶点集合为VV;初始令集合u={s},v=V−uu={s},v=V−u;
2.选择一条权值最小的边(u0,v0),把点加入到最小生成树中
3.重复上述步骤,直到最小生成树有n-1条边或者n个顶点为止。
由于不断向集合u中加点,所以最小代价边必须同步更新;需要建立一个辅助数组,用来维护集合中每个顶点与最小代价边信息。
代码:
#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
const int N=10005;
const int inf=0x3f3f3f;
int n,vis[N],dis[N]; //dis数组用来记录点i到生成树的最小距离
int mp[105][105];
struct node
{
int u,v,w;
}p[N];
void Prim()
{
for(int i=1;i<=n;i++)
dis[i]=mp[1][i];
vis[1]=1;
int sum=0;
for(int i=1;i<n;i++)
{
int minn=inf,k;
for(int j=1;j<=n;j++)
if(vis[j]==0&&dis[j]<minn)
minn=dis[j],k=j;
vis[k]=1;
sum+=dis[k];
for(int j=1;j<=n;j++)
if(vis[j]==0&&mp[k][j]<dis[j])
dis[j]=mp[k][j];
}
printf("%d\n",sum);
}
int main()
{
while(~scanf("%d",&n)&&n)
{
memset(mp,0,sizeof(mp));
memset(vis,0,sizeof(vis));
memset(dis,0,sizeof(dis));
int k=n*(n-1)/2;
for(int i=0;i<k;i++)
{
scanf("%d%d%d",&p[i].u,&p[i].v,&p[i].w);
mp[p[i].u][p[i].v]=mp[p[i].v][p[i].u]=p[i].w;
}
Prim();
}
return 0;
}