1. 程式人生 > >Qin Shi Huang's National Road System(列舉+最小瓶頸路)

Qin Shi Huang's National Road System(列舉+最小瓶頸路)

【題意】
       秦朝有n個城市,需要修建一些道路使得任意兩個城市之間連通,道士徐福可以用法術修路,不用花錢,但是隻能修一條路,因此需要慎重選擇道士徐福要修哪條路。秦始皇不僅希望其他道路的總長度B儘量短,同時還希望徐福修的路連線的兩個城市的人口之和A儘量大,因此要找到一個A/B最大的修路方案,並求出這個值。

【輸入格式】
       第一行是資料組數t(t<=10)每組資料第一行為城市數目n(3<=n<=1000)隨後n行每行有三個整數,x,y,p(0<=x,y<=1000 1<=p<=100000)分別是該城市的笛卡爾座標和人口數量,保證城市的位置不會互相重疊。

【輸出格式】
       A/B的最大值,保留兩位小數

【思路】
       因為n是1000的規模,所以我們可以考慮列舉道士徐福所修路的端點u,v 在此之前如果我們已經把最小生成樹和任意兩點的最小瓶頸路求出來之後,那麼每次的答案就是最小生成樹的權值減去u,v的最小瓶頸路(法術把u,v連通所以可以刪掉u,v之間最長的哪一條邊),取最小值就是答案。

#include<bits/stdc++.h>
using namespace std;

const int maxn = 1050;

struct Edge {
    int from, to;
    double dist;
    Edge(int f = 0, int t = 0, double
d = 0) :from(f), to(t), dist(d) {} bool operator<(const Edge& e)const { return dist < e.dist; } }; int n, m; int par[maxn]; bool used[maxn]; int x[maxn], y[maxn], p[maxn]; double f[maxn][maxn]; vector<Edge> edges, g[maxn]; double dis(int x1, int y1, int x2, int y2) { return
sqrt((x1 - x2)*(x1 - x2) + (y1 - y2)*(y1 - y2)); } int find(int x) { return x == par[x] ? x : par[x] = find(par[x]); } double kruscal() { for (int i = 0; i < n; ++i) { par[i] = i; g[i].clear(); } int cnt = 0; double mst = 0; for (int i = 0; i < m; ++i) { int u = edges[i].from; int v = edges[i].to; int x = find(u); int y = find(v); if (x != y) { par[x] = y; mst += edges[i].dist; ++cnt; g[u].push_back(Edge(u, v, edges[i].dist)); g[v].push_back(Edge(v, u, edges[i].dist)); if (cnt == n - 1) break; } } return mst; } void dfs(int u) { used[u] = 1; for (int i = 0; i < g[u].size(); ++i) { int v = g[u][i].to; if (!used[v]) { for (int x = 0; x < n; ++x) { if (used[x]) f[x][v] = f[v][x] = max(f[u][x], g[u][i].dist); } dfs(v); } } } int main() { int t; scanf("%d", &t); while (t--) { scanf("%d", &n); for (int i = 0; i < n; ++i) scanf("%d%d%d", &x[i], &y[i], &p[i]); edges.clear(); for (int i = 0; i < n; ++i) { for (int j = i + 1; j < n; ++j) { double d = dis(x[i], y[i], x[j], y[j]); edges.push_back(Edge(i, j, d)); } } m = edges.size(); sort(edges.begin(), edges.end()); double mst = kruscal(); memset(used, 0, sizeof(used)); for (int i = 0; i < n; ++i) for (int j = 0; j < n; ++j) f[i][j] = 0; dfs(0); double ans = 0; for (int u = 0; u < n; ++u) { for (int v = u + 1; v < n; ++v) { double A = p[u] + p[v]; double B = mst - f[u][v]; ans = max(ans, A / B); } } printf("%.2lf\n", ans); } return 0; }