图的绝对中心&最小直径生成树 GYM102391 I -Minimum Diameter Spanning Tree

题面
参考：
窝们要枚举绝对中心在哪一条边上
假设绝对中心S在边<u,v>上 离u的距离是x
根据定义 他到其他点的最短路的最大值最小 且这个最大值会出现两次
把他写成函数的形式：

$$d(i,x)=min(d[u][i]+x,d[v][i]+w-x)$$,其中x是S到u的距离，$w$是边$<u,v>$的边权
而窝们要对所有的$i$ 求出$$f(x)=\max d(i,x)$$
做出d的图像 可以发现它是-|x|的图像平移后得到的 在区间[0,w]上的部分
而它的纵截距是d[u][i], 在直线x=w上的截距是d[v][i]

而f的图像就是这样若干个折线构成的图像的最上边

观察图像可以发现所求的f的最小值在交点处或区间端点处取得
考虑怎么求交点
把点按到u的距离降序排序的话 对答案有贡献的点到v的距离一定是升序排序的
否则那个点的d图像会被其他点的d图像压在下面
那么按到u的距离降序排个序 然后扫一遍 用lst表示上一个有贡献的i
对于当前点cur 若d[u][cur]<=d[u][lst] 且d[v][cur]<=d[v][lst] 那么cur就无用了

否则若d[v][cur]>d[v][lst] 就有d[u][cur]+x==d[v][lst]+w - x 得x == (d[v][lst]-d[u][cur]+w)/2
但窝们要求的是f的最小值 甚至可以不求出x 直接在符合条件时用d[u][cur]+d[v][lst]+w和直径比大小就行了
取得最小值 就更新一下答案
端点处用d[u][rk[u][n]] + d[v][rk[v][n]]+w更新答案
在更新答案的时候维护一下是枚举哪一条边的时候取到了这个答案
然后以这条边为起点 建出最短路径树
可以得到一颗最小直径生成树

#include<bits/stdc++.h>
using namespace std;
const int maxn = 507;
#define ll long long
const ll inf = 1ll << 61;
int n, m;
struct edge {
    int u, v;
    ll w;
}e[maxn * maxn];
ll d[maxn][maxn], g[maxn][maxn];
vector<int> adj[maxn];
int pre[maxn];
int rk[maxn][maxn]; //rk[i][j]表示距离点i第j近的点
ll tmp[maxn], ans, dis[maxn], X;
bool cmp(int a, int b) { return tmp[a] < tmp[b];}
int rd() {
    int s = 0, f = 1; char c = getchar();
    while (c < '0' || c > '9') {if (c == '-') f = -1; c = getchar();}
    while (c >= '0' && c <= '9') {s = s * 10 + c - '0'; c = getchar();}
    return s * f;
}
typedef pair<int, int> pii;
set<pii, less<pii> > min_heap, res;
#define mk make_pair
void dij(int s1, int s2) {
    for (int i = 1; i <= n; i++) dis[i] = inf;
    dis[s1] = X; dis[s2] = g[s1][s2] - X;
    min_heap.insert(mk(dis[s1], s1));
    min_heap.insert(mk(dis[s2], s2));
    if (s1 != s2) {
        pre[s2] = s1;
    }
    while(!min_heap.empty()) {
        auto it = min_heap.begin();
        int u = it->second;
        min_heap.erase(*it);
        for (int i = 0; i < adj[u].size(); i++) {
            int v = adj[u][i];
            if (dis[v] > dis[u] + g[u][v]) {
                min_heap.erase(mk(dis[v], v));
                dis[v] = dis[u] + g[u][v];
                min_heap.insert(mk(dis[v], v));
                res.erase(mk(pre[v], v));
                pre[v] = u;
                res.insert(mk(pre[v], v));
            }
        }
    }
}
int main() {
    n = rd(); m = rd();
    for (int i = 1; i <= n; i++) for (int j = 1; j <= n; j++) d[i][j] = inf;
    for (int i = 1; i <= n; i++) d[i][i] = 0;
    for (int i = 1; i <= m; i++) {
        int u, v, w;
        u = e[i].u = rd(); v = e[i].v = rd(); w = e[i].w = g[u][v] = g[v][u] = rd()*2ll;
        d[u][v] = d[v][u] = w;
        adj[u].push_back(v); adj[v].push_back(u);
    }
    for (int k = 1; k <= n; k++)
        for (int i = 1; i <= n; i++) 
            for (int j = 1; j <= n; j++) 
                d[i][j] = min(d[i][j], d[i][k] + d[k][j]);
    for (int i = 1; i <= n; i++) {
        for (int j = 1; j <= n; j++) {
            rk[i][j] = j; tmp[j] = d[i][j];
        }
        sort(rk[i] + 1, rk[i] + n + 1, cmp);
    }
    ans = inf;//直径
    int s1, s2;
    for (int k = 1; k <= m; k++) {
        int u = e[k].u, v = e[k].v;
        if (d[u][rk[u][n]] == d[u][rk[u][n-1]] && d[u][rk[u][n]] * 2ll < ans) {
            ans = d[u][rk[u][n]] * 2ll;//u点作为中心
            s1 = s2 = u;
        }
        if (d[v][rk[v][n]] == d[v][rk[v][n-1]] && d[v][rk[v][n]] * 2ll < ans) {
            ans = d[v][rk[v][n]] * 2ll; 
            s1 = s2 = v;
        }
        int lst, cur;
        for (cur = n-1, lst = n; cur >= 1; cur--) {//d[u][i]减小 合法的d[v][i]应对应增加
            if (d[v][rk[u][lst]] < d[v][rk[u][cur]]) {
                if (ans > d[u][rk[u][cur]] + d[v][rk[u][lst]] + e[k].w) {
                    ans = d[u][rk[u][cur]] + d[v][rk[u][lst]] + e[k].w;
                    X = ans / 2 - d[u][rk[u][cur]];
                    s1 = u, s2 = v;                    
                }
                lst = cur;
            }
        }
    }
    cout << ans / 2ll << endl;
    dij(s1, s2);
    for (auto it = res.begin(); it != res.end(); ++it) {
        printf("%d %d\n", it->first, it->second);
    }
    if (s1 != s2) printf("%d %d\n",s1, s2);
    return 0;
}