POJ2349 prim算法

#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
double x[1001],y[1001],a[1001][1001],d[1001],p[1001],ans[10001],asd[10001];
int main(){
int i,j,k,m,n;
int t;
cin>>t;
int v=0;
for(int q=1;q<=t;q++){
int u=0;
cin>>m>>n;
for(i=1;i<=n;i++)
for(j=1;j<=n;j++){
a[i][j]=0;
}
for(i=1;i<=n;i++){
cin>>x[i]>>y[i];
p[i]=0;d[i]=10000000.0;
}
for(i=1;i<=n;i++)
for(j=1;j<=n;j++){
a[i][j]=sqrt((x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j]));
}
d[1]=0;
int min;
for(i=1;i<=n;i++){
min=1000000000.0;
for(j=1;j<=n;j++)
if(!p[j] && d[j]<min){
k=j;
min=d[j];
}
p[k]=1;
ans[++u]=d[k];
for(j=1;j<=n;j++)
if(!p[j] && a[i][j]!=0 && d[j]>a[j][k])
d[j]=a[j][k];
}
for(i=1;i<u;i++)
for(j=i+1;j<=u;j++){
if(ans[i]<ans[j]){
int tmp;
tmp=ans[i];ans[i]=ans[j];ans[j]=tmp;
}
}
asd[++v]=ans[m];
for(i=1;i<=u;i++)ans[i]=0;
}
for(i=1;i<=t;i++){
cout<<asd[i]<<endl;
}
return 0;
}