HEU Tunnels(常规最小生成树)

注意数据类型,使用double!

#include <cstdio>
#include <cmath>
#include <cstring>

const int N = 110;
struct Node {
    double x, y;
}point[N];
int n;
double g[N][N];
const double INF = 10000000000.0;
const double eps = 1e-9;

double prim() {
    double ans = 0.0, mi;
    bool vis[N]; 
    memset( vis, 0, sizeof(vis));
    vis[0] = true;
    int v;
    for ( int u = 0; u < n-1; ++u ) {
        mi = INF;
        for ( int i = 0; i < n; ++i ) if ( !vis[i] && mi - g[0][i] > eps ) mi = g[0][i], v = i;
        ans += sqrt(mi);
        vis[v] = true;
        for ( int i = 0; i < n; ++i ) if ( !vis[i] && g[0][i] - g[v][i] > eps ) g[0][i] = g[v][i];
    }
    return ans;
}
double dis( int i, int j ) {
    return (point[i].x - point[j].x)*(point[i].x - point[j].x) + (point[i].y - point[j].y)*(point[i].y - point[j].y);
}
int main()
{
    while ( scanf("%d", &n) != EOF && n ) {
        for ( int i = 0; i < n; ++i ) 
            scanf("%lf%lf", &point[i].x, &point[i].y);
        for ( int i = 0; i < n; ++i )
            for ( int j = i; j < n; ++j ) 
                if ( i == j ) g[i][j] = INF;
                else g[i][j] = g[j][i] = dis( i, j );
        printf("%.2lf\n", prim());
    }
}



  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值