【NOIP2016提高A组集训第12场11.10】图的半径 解方程找最优解 枚举关键点

题目大意

给你一幅$n$个点，$m$条边的无向图。现在要求你找到一个点，可以在边上。是节点中到达这和点的最大值最小，并且输出这个最小值。

$n\leq200$
$m\leq19900$

解题思路

可以直接枚举每条边，找出每条边的最有点，然后再在所有边中找到最有的边。
对于一条连接$u,v$的边。设$a_i$表示从节点$i$到$u$的最短路$b_v$为节点$i$到$v$的最短路。这可可以用$floyd$预处理。那么这条边上必定存在一个点。设我们选的位置是$x$，那么这个点到$x$的最短路就是$\max(a_i+x,b_i+w-x)$，用$A_x$记录符合$a_i+x\leq b_i+w-x$的点，$B_x$表示$A_x$的补集。那么对于这个点的答案就是$\max(\max(A_x)+x,max(B_x)+w-x)$，那么我们可以枚举这$n$个分界点，一种集合划分的情况就对应着两个分界点之间的区域。那么现在的问题是在这段区域找到最优点。很显然当满足$\max(a_i)+x=\max(b_i)+w-x$时$x$取到最优优值。在判断$x$是否在这个区域。不在就取端点的最小值。然后动态维护这两个集合就行了。

程序

//YxuanwKeith
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <cmath>

using namespace std;

const int MAXN = 205, MAXM = 2e4 + 5;

struct Node {
    int u, v, w;
} p[MAXM];

struct Node2 {
    double num;
    int id;
} t[MAXN];

int n, m, dis[MAXN][MAXN], mx[MAXN];
double ans;

bool cmp(Node2 a, Node2 b) {
    return a.num < b.num;
}

void floyd() {
    for (int k = 1; k <= n; k ++)
        for (int i = 1; i <= n; i ++) {
            if (i == k) continue;
            for (int j = 1; j <= n; j ++) {
                if (i == j || k == j) continue;
                dis[i][j] = min(dis[i][j], dis[i][k] + dis[k][j]);
            }               
        }
}

double calc(double s, int m1, int m2, int w) {
    return max(m1 + s, m2 + w - s);
}

void solve() {
    for (int i = 1; i <= m; i ++) {
        int u = p[i].u, v = p[i].v, w = p[i].w;
        for (int j = 1; j <= n; j ++) {
            t[j].num = 1.0 * (dis[v][j] + w - dis[u][j]) / 2;
            t[j].id = j;
        }
        sort(t + 1, t + n + 1, cmp);
        mx[n + 1] = 0;
        for (int j = n; j; j --) mx[j] = max(mx[j + 1], dis[u][t[j].id]);
        int s = 0;
        t[n + 1].num = w;
        for (int j = 1; j <= n + 1; j ++) {
            double side = 1.0 * (s - mx[j] + w) / 2;
            if (side >= t[j - 1].num && side <= t[j].num) ans = min(ans, calc(side, mx[j], s, w)); else
                ans = min(ans, max(calc(t[j - 1].num, mx[j], s, w), calc(t[j].num, mx[j], s, w)));
            s = max(s, dis[v][t[j].id]);    
        }
    }
}

int main() {
    scanf("%d%d", &n, &m);
    double l = 0, r = 0;
    memset(dis, 60, sizeof dis);
    for (int i = 1; i <= m; i ++) {
        int u, v, w;
        scanf("%d%d%d", &u, &v, &w);
        dis[u][v] = min(dis[u][v], w);
        dis[v][u] = min(dis[v][u], w);
        r = r + w;
        p[i].u = u, p[i].v = v, p[i].w = w;
    }
    floyd();
    for (int i = 1; i <= n; i ++) dis[i][i] = 0;
    ans = 1e9 + 7;
    solve();
    printf("%.2lf", ans);
}