HDU-1875(畅通工程再续，最小生成树)

畅通工程再续 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 37053    Accepted Submission(s): 12398   Problem Description 相信大家都听说一个“百岛湖”的地方吧，百岛湖的居民生活在不同的小岛中，当他们想去其他的小岛时都要通过划小船来实现。现在政府决定大力发展百岛湖，发展首先要解决的问题当然是交通问题，政府决定实现百岛湖的全畅通！经过考察小组RPRush对百岛湖的情况充分了解后，决定在符合条件的小岛间建上桥，所谓符合条件，就是2个小岛之间的距离不能小于10米，也不能大于1000米。当然，为了节省资金，只要求实现任意2个小岛之间有路通即可。其中桥的价格为 100元/米。     Input 输入包括多组数据。输入首先包括一个整数T(T <= 200)，代表有T组数据。 每组数据首先是一个整数C(C <= 100),代表小岛的个数，接下来是C组坐标，代表每个小岛的坐标，这些坐标都是 0 <= x, y <= 1000的整数。     Output 每组输入数据输出一行，代表建桥的最小花费，结果保留一位小数。如果无法实现工程以达到全部畅通，输出”oh!”.     Sample Input   2 2 10 10 20 20 3 1 1 2 2 1000 1000     Sample Output   1414.2 oh!     Author 8600     Source 2008浙大研究生复试热身赛（2）——全真模拟

Prime算法：

#include<iostream>
#include<cmath>
#include<algorithm>
#include<climits>
#include<cstring>
using namespace std;                                 
const int MAX=2000+10;
const double INF=100000007.0;
struct point{
	double x;
	double y;
};
double diss(point a,point b){
	return sqrt((a.x -b.x)*(a.x -b.x)+(a.y -b.y)*(a.y -b.y));
}
double Map[MAX][MAX];              			
double Prim(int from,int to,int n){
	double dis[MAX];                          	
	bool vis[MAX];                		
	memset(vis,false,sizeof(vis));
 
	    for(int i=1;i<=n;i++)
	    dis[i]=Map[1][i];            		   
	
		vis[from]=true;
		double sum=0.0,Min;
		int f;
	    for(int i=1;i<=n-1;i++){
	       Min=INF;
		
		for(int k=1;k<=n;k++){
			 if(!vis[k] && dis[k]<Min)   
			  Min=dis[k],f=k;              
		    }
			if(Min==INF) return INF; 
		   vis[f]=true;
		   sum+=dis[f];
		for(int k=1;k<=n;k++){
	        if(!vis[k]&&Map[f][k]<dis[k])
	        dis[k]=Map[f][k];
		}		
	}
	 return sum;               
}
int main(){
	int n,t;
	cin>>t;
	while(t--){
	 cin>>n;
	 point p[1200];
     for(int i=1;i<=n;i++){
    	for(int j=1;j<=n;j++){
    		if(i==j)
    		Map[i][j]=0.0;
    		else 
    		Map[i][j]=INF;
		}
	}

	for(int i=1;i<=n;i++){
		cin>>p[i].x>>p[i].y;
	}
	for(int i=1;i<=n-1;i++){
	    for(int j=i+1;j<=n;j++){
	    	double temp=diss(p[i],p[j]);
	    	if(temp<10.0||temp>1000.0)
	    	Map[i][j]=Map[j][i]=INF;
	    	else 
			Map[i][j]=Map[j][i]=temp;
		}
	}
	
	    double ans=Prim(1,n,n);
	    if(ans==INF) cout<<"oh!"<<endl;
    	else{
    		printf("%.1lf\n",ans*100.0);
		}
	}
	return 0;
}

再来个看起来清爽一点的版本：

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int MAX_N=2050;
const double INF=100000007.0;

struct point{
	double x;
	double y;
};
double dist(point a,point b){
	return sqrt((a.x -b.x)*(a.x -b.x)+(a.y -b.y)*(a.y -b.y));
}
double mincost[MAX_N];
bool used[MAX_N];
double cost[MAX_N][MAX_N];
int V;
double Prim(){
	for(int i=1;i<=V;i++){
		mincost[i]=cost[1][i];   //初始化的时候注意一下方式，这个不能直接全部INF 
		used[i]=false;
	}
	mincost[0]=0.0;
	double res=0.0;
	while(true){
		int v=-1;
		for(int u=1;u<=V;u++){
			if(!used[u] && (v == -1||mincost[u]<mincost[v]))
			v=u;
		}
		if(v==-1) break;
		
		used[v]=true;
		res+=mincost[v];
		for(int u=1;u<=V;u++){
			mincost[u]=min(mincost[u],cost[v][u]);
		}
	}
	return res;
}
int main(){
	int t;
	cin>>t;
	while(t--){
		cin>>V;
		point p[1200];
		for(int i=1;i<=V;i++){
			for(int j=1;j<=V;j++){
				cost[i][j]=(i==j? 0.0:INF);
			}
		}
		for(int i=1;i<=V;i++){
			cin>>p[i].x>>p[i].y;
		}
		
		for(int i=1;i<=V-1;i++){
			for(int j=i+1;j<=V;j++){
				double temp=dist(p[i],p[j]);
				cost[i][j]=cost[j][i]=((temp<10.0||temp>1000.0)? INF:temp);
			}
		}
		double ans=Prim();
		if(ans>=INF)cout<<"oh!"<<endl;
		else
		printf("%.1lf\n",ans*100);
	}
	return 0;
}