***************************************转载请注明出处:http://blog.csdn.net/lttree***************************************
畅通project
Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 15572 Accepted Submission(s): 6462
Problem Description
省*“畅通project”的目标是使全省不论什么两个村庄间都能够实现公路交通(但不一定有直接的公路相连。仅仅要能间接通过公路可达就可以)。经过调查评估。得到的统计表中列出了有可能建设公路的若干条道路的成本。
现请你编敲代码,计算出全省畅通须要的最低成本。
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
?
Source
最小生成树的基础题目。畅通project。
赤裸裸的求最小生成树。
额外加了一点要推断 能否构成最小生成树。
这次,我用的Kruskal算法。
Kruskal 构建最小生成树:
大体。就是先依照边长进行排序(由小到大),
然后再向外加边。
加边的时候推断是否能构成回路,假设能构成回路,就不能加边。
为什么这么做是对的呢?
首先。要知道,最小生成树,一定不会出现回路!
Why?自己算算。。
o(╯□╰)o。
。
。
然后,我们已经将边依照小到大排序了,所以这样加边,得到的肯定是最小生成树啦~
Kruskal算法重要的就是推断回路。
这个是用 并查集 来实现的,(并查集相关可戳:http://blog.csdn.net/lttree/article/details/23820679)
然后。最后再用并查集Find函数来找找。是否全部的点都在同一个集合,假设不在,输出?
恩,OK~
/****************************************
*****************************************
* Author:Tree *
*From :http://blog.csdn.net/lttree *
* Title : 畅通project *
*Source: hdu 1863 *
* Hint : 最小生成树(Kruskal) *
*****************************************
****************************************/ #include <stdio.h>
#include <algorithm>
using namespace std;
struct EDGE
{
int u,v,cost;
}eg[100001];
int n,m,father[100001]; bool cmp(EDGE e1,EDGE e2)
{
return e1.cost<e2.cost;
} // 并查集 初始化函数
void Init( int m )
{
int i;
for(i=1;i<=m;i++)
father[i]=i;
}
// 并查集 查找函数
int Find(int x)
{
while(father[x]!=x)
x=father[x];
return x;
}
// 并查集 合并函数
void Combine(int a,int b)
{
int temp_a,temp_b;
temp_a=Find(a);
temp_b=Find(b); if(temp_a!=temp_b)
father[temp_a]=temp_b;
} // 最小生成树 Kruskal 算法
int Kruskal( void )
{
EDGE e;
int i,res;
sort(eg,eg+n,cmp);
// 并查集 初始化
Init(m); // 构建最小生成树
res=0;
for( i=0;i<n;++i )
{
e=eg[i];
if( Find(e.u)!=Find(e.v) )
{
Combine(e.u,e.v);
res+=e.cost;
}
}
return res;
} int main()
{
int i,ans;
bool bl;
while( scanf("%d%d",&n,&m) && n )
{
for( i=0;i<n;++i )
scanf("%d%d%d",&eg[i].u,&eg[i].v,&eg[i].cost);
ans=Kruskal(); // 是否全部的点都在同一个集合
bl=true;
for(i=2;i<=m;++i)
if( Find(1)!=Find(i) )
{
bl=false;
break;
} if( bl ) printf("%d\n",ans);
else printf("?\n");
}
return 0;
}
又搞了搞最小生成树的Prim算法。
。。
Prim算法就是 从一个点慢慢扩展到全图。
原理:
就是从一个点出发,然后从全部与这个点直接相连的点中。找权值最小的那条边。进行扩展。
可是。不用每次都寻找,仅仅须要在增加一个点后,
更新这个集合到其它各个点的距离,就可以。
Prim和Kruskal差别:
我的理解是:
Prim是一个集合战斗,慢慢扩展,一个个吞并,最后构成一个大树。
而Kruskal 是多个集合(≥1,也可能是一个集合) 分别作战,最后合并成一个大树。
Prim和Kruskal优劣性:
Prim须要每次都维护mincost数组(距离各个点的最短距离),所以须要O(n^2)
可是同最短路的Dijkstra一样,假设用 堆 来维护,则复杂度可降到 O(n log n)
Kruskal算法仅仅是在排序上最费时,算法复杂度可看做 O( n log n )
本题的Prim算法:
/****************************************
*****************************************
* Author:Tree *
*From :http://blog.csdn.net/lttree *
* Title : 畅通project *
*Source: hdu 1863 *
* Hint : 最小生成树(Prim ) *
*****************************************
****************************************/
#include <stdio.h>
#include <string.h> #define RANGE 101
#define MAX 0x7fffffff
int cost[RANGE][RANGE];
int mincost[RANGE];
bool used[RANGE]; // n个点。m条边
int n,m; int Min(int a,int b)
{
return a<b?a:b;
} void prim( void )
{
// sum 记录最小生成树权值
int i,v,u,sum;
// 从1到各个点距离,初始化used数组
for( i=1;i<=n;++i )
{
used[i]=false;
mincost[i]=cost[1][i];
}
sum=0; while( true )
{
v=-1; // 从没有连接到的点中,找近期的点
for( u=1;u<=n;++u )
if( !used[u] && (v==-1 || mincost[u]<mincost[v]) )
v=u; if( v==-1 ) break;
if( mincost[v]==MAX ) break; used[v]=true;
sum+=mincost[v]; // 更新到各个点的距离
for( u=1;u<=n;++u )
mincost[u]=Min( mincost[u],cost[v][u] );
} // 推断是否能构成最小生成树
for( i=1;i<=n;++i )
{
if( used[i]==false )
{
printf("?\n");
return;
}
}
printf("%d\n",sum);
} int main()
{
int i,j;
int u,v,c; while( scanf("%d%d",&m,&n) && m )
{
// init cost by MAX
for( i=1;i<=n;++i )
for( j=1;j<=i;++j )
{
if( i==j ) cost[i][j]=0;
else cost[i][j]=cost[j][i]=MAX;
}
for( i=0;i<m;++i )
{
scanf("%d%d%d",&u,&v,&c);
cost[u][v]=cost[v][u]=c;
} prim();
}
return 0;
}
版权声明:本文博客原创文章,博客,未经同意,不得转载。