HDU 1233はやはりスムーズエンジニアリング(最小生成ツリー)


やはりスムーズな工事です
Time Limit: 4000/2000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others) Total Submission(s): 37570    Accepted Submission(s): 16923
Problem Description
ある省は田舎の交通状況を調査し、得られた統計表には任意の2つの村間の距離がリストされている.省政府の「円滑な工事」の目標は、全省のどの2つの村の間でも道路交通を実現できるようにすることである(しかし、直接的な道路がつながっているとは限らず、間接的に道路を通過すればよい)、敷設された道路の総長さを最小限に抑えることを要求している.最小の道路の全長を計算してください.
 
Input
テスト入力には、いくつかのテスト例が含まれます.各試験例の1行目は、村の数N(<100)を与える.その後のN(N−1)/2行は村間の距離に対応し,各行には2つの村の番号とこの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
Hint
Huge input, scanf is recommended.

 
Source
浙大コンピュータ大学院生の再試験の上機試験-2006年
 
題意:各村を結ぶ最短道路、典型的な最小生成木を要求する図である.
TLEコードを...
#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<string>
#include<cstdlib>
#include<iomanip>
#include<algorithm>
#include<time.h>
typedef long long LL;
using namespace std;
#define MIN 65535
#define MAX_Point 120   //      
#define MAX_Edge 14400  //      
int flag1 =0;
double sum;
double arr_list[MAX_Point][MAX_Point]; 
struct Edge
{
	int point;
	double lowcost;
	int flag;
};
Edge edge[MAX_Edge];
double prim(int n)
{
	int i,j,k=1,flag;
	double min,sum2=0;
	j=1;
	for(i=1;i<=n;i++)
	{
		if(i!=j)
		{
			edge[i].point=i;
			edge[i].lowcost=arr_list[j][i];
			edge[i].flag=0;
		}
	}
	edge[j].lowcost=0;
	edge[j].flag=1;
	for(i=2;i<=n;i++)
	{
		k=1;
		min=MIN;
		flag=0;
		for(j=2;j<=n;j++)
		{
			if(edge[j].flag==0 && edge[j].lowcost<min)
			{
				k=j;
				min=edge[j].lowcost;
				flag=1;
			}
		}
		if(!flag) return -1;
		sum2+=min;
		edge[k].flag=1;
		for(j=2;j<=n;j++)
		{
			if(edge[j].flag==0 && arr_list[k][j]<edge[j].lowcost)
			{
				edge[j].point=k;
				edge[j].lowcost=arr_list[k][j];
			}
		}
	}
	return sum2;
}	
int main()
{
	double t,time_data;
	int n,m;
	int a,b;
	while(cin>>m)
	{
		if(m==0) break;
		sum=0;
		for(int i=1;i<=m;i++)
		  for(int j=1;j<=m;j++)
		   arr_list[i][j]=65535;
		   for(int i=1;i<=m*(m-1)/2;i++)
		   {
		   	cin>>a>>b>>time_data;
		   	if(time_data<arr_list[a][b]) //     ,         
             {
             	arr_list[a][b]=time_data;
				arr_list[b][a]=time_data; 
             }
		   }
		   sum=prim(m);
		   printf("%.lf
",sum); } return 0; }

マルチポーズACコード:
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define N 110
#define INF 0x7ffffff
using namespace std;
int n;
int dist[N],map[N][N],isvisited[N];
int prim(){
         int i,j,min,pos;
         int sum=0;
         memset(isvisited,false,sizeof(isvisited));
         //   
         for(i=1;i<=n;i++){
                   dist[i]=map[1][i];
         } 
        // 1  
         isvisited[1]=true;
         dist[1]=65535;
         //             
         for(i=1;i<n;i++){
                   min=65535;
                   for(j=1;j<=n;j++){
                            if(!isvisited[j] && dist[j]<min){
                                     min=dist[j];
                                     pos=j;
                            }
                   }        
                   sum+=dist[pos]; //    
                   isvisited[pos]=true; 
                   //    
                  for(j=1;j<=n;j++){
                            if(!isvisited[j] && dist[j]>map[pos][j]){

                                     dist[j]=map[pos][j];
                            }
                  }
         }        
         return sum;
}
int main()
{
    while(~scanf("%d",&n)&&n)
    {
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                map[i][j]=INF;
        int m=n*(n-1)/2;
        for(int i=0;i<m;i++)
        {
            int u,v,t;
            scanf("%d%d%d",&u,&v,&t);
            if(map[u][v]>t) 
			 map[u][v]=map[v][u]=t;
        }
        cout<<prim()<<endl;
    }
}

ACコード:
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define N 110
#define INF 0x7ffffff
using namespace std;
int v[N],d[N],mp[N][N],n;
int prim()
{
    for(int i=1;i<=n;i++) 
   {
    	 d[i]=INF,v[i]=0;
    } 
    d[1]=0;
    for(int i=0;i<n;i++)
    {
        int mmin=INF,t;
        for(int j=1;j<=n;j++){
        	if(!v[j]&&mmin>d[j])
                mmin=d[t=j];
        }   
        v[t]=1;
        for(int j=1;j<=n;j++)
        {
        	if(!v[j])
                d[j]=min(d[j],mp[t][j]);
        }          
    }
    int ans=0;
    for(int i=1;i<=n;i++)    
	   ans+=d[i];
    return ans;
}

int main()
{
    while(~scanf("%d",&n)&&n)
    {
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                mp[i][j]=INF;
        int m=n*(n-1)/2;
        for(int i=0;i<m;i++)
        {
            int u,v,t;
            scanf("%d%d%d",&u,&v,&t);
            if(mp[u][v]>t) 
			 mp[u][v]=mp[v][u]=t;
        }
        cout<<prim()<<endl;
    }
}