Install Air Conditioning HDU - 4756（最小生成树+树形dp）

#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
using namespace std;

const int maxn = 1e3 + 10;
const double inf = 0x3f3f3f3f3f3f;
double d[maxn][maxn], lowc[maxn], dp[maxn][maxn];
bool vis[maxn], is_tree[maxn][maxn];
int pre[maxn], head[maxn];
int cnt, n, m;
double sum, ans;

struct Point {
    double x, y;
} a[maxn];
struct Edge {
    int to, next;
} edge[maxn<<1];

inline double distan(const Point& lhs, const Point& rhs) {
    return sqrt((lhs.x - rhs.x) * (lhs.x - rhs.x) + (lhs.y - rhs.y) * (lhs.y - rhs.y));
}

inline void addedge(int u, int v) {
    edge[cnt].to = v;
    edge[cnt].next = head[u];
    head[u] = cnt++;
}

void Prim() {
    sum = 0.0;
    memset(vis, 0, sizeof(vis));
    memset(pre, 0, sizeof(pre));
    for (int i = 1; i < n; i++) lowc[i] = d[0][i];
    vis[0] = true;
    for (int i = 1; i < n; i++) {
        double minc = inf;
        int p = -1;
        for (int j = 0; j < n; j++) {
            if (!vis[j] && minc > lowc[j]) {
                minc = lowc[j];
                p = j;
            }
        }
        sum += minc;
        vis[p] = true;
        addedge(p, pre[p]);
        addedge(pre[p], p);
        for (int j = 0; j < n; j++) {
            if (!vis[j] && lowc[j] > d[p][j]) {
                lowc[j] = d[p][j];
                pre[j] = p;
            }
        }
    }
}

double dfs(int u, int fa, int root) {
    double ans = fa == root ? inf : d[root][u];
    for (int i = head[u]; ~i; i = edge[i].next) {
        int v = edge[i].to;
        if (v == fa) continue;
        double temp = dfs(v, u, root);
        ans = min(ans, temp);
        dp[u][v] = dp[v][u] = min(dp[u][v], temp);
    }
    return ans;
}

void dfs1(int u, int fa) {
    for (int i = head[u]; ~i; i = edge[i].next) {
        int v = edge[i].to;
        if (v == fa) continue;
        if (fa) ans = max(ans, sum - d[u][v] + dp[u][v]);
        dfs1(v, u);
    }
}

int main() {
    int T;
    scanf("%d", &T);
    while (T--) {
        scanf("%d%d", &n, &m);
        for (int i = 0; i < n; i++) scanf("%lf%lf", &a[i].x, &a[i].y);
        for (int i = 0; i < n; i++) {
            for (int j = i + 1; j < n; j++) {
                d[i][j] = d[j][i] = distan(a[i], a[j]);
            }
        }
        cnt = 0;
        memset(head, -1, sizeof(head));
        Prim();
        for (int i = 0; i < n; ++i) 
            for (int j = 0; j < n; ++j) dp[i][j] = inf;
        for (int i = 0; i < n; i++) dfs(i, -1, i);
        ans = sum;
        dfs1(0, 0);
        ans *= 1.0 * m;
        printf("%.2f\n", ans);
    }
    return 0;
}