luogu1967 货车运输

似乎是什么最大瓶颈路径（顾名思义），它存在于最大生成树中。
证明：

如果最大瓶颈路径不存在与最大生成树中。
这些不在最大生成树中的边会和最大生成树形成环。
我们删掉环上最小的边，保留这一条边，会得到一棵新的更大的生成树。
这与原来那棵树是最大生成树矛盾了。
注意，最短路不一定在最小生成树上（如一个环的情况）
摘自洛谷题解。作者crf1596028912，摘自高天宇学长SD夏令营的PPT
我们用和倍增lca差不多的方法求之。
这代码贼长。

#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
int n, m, grand[10005][15], minq[10005][15], fa[10005], hea[10005], cnt;
int r1, r2, deep[10005], q, uu, vv;
struct Ed{
    int fro, too, val;
}ed[50005];
struct Edge{
    int too, nxt, val;
}edge[50005];
int rn(int &x){
    char ch=getchar();
    x = 0;
    while(ch<'0' || ch>'9') ch = getchar();
    while(ch>='0' && ch<='9'){
        x = x*10 + ch - '0';
        ch = getchar();
    }
    return x;
}
void init(){
    rn(n);
    rn(m);
    for(int i=1; i<=m; i++){
        rn(ed[i].fro);
        rn(ed[i].too);
        rn(ed[i].val);
    }
    for(int i=1; i<=n; i++)
        fa[i] = i;
}
bool cmp(Ed x, Ed y){
    return x.val>y.val;
}
int myfind(int x){
    return fa[x]==x?x:fa[x]=myfind(fa[x]);
}
void add_edge(int fro, int too, int val){
    edge[++cnt].too = too;
    edge[cnt].nxt = hea[fro];
    edge[cnt].val = val;
    hea[fro] = cnt;
}
void kruskal(){
    sort(ed+1, ed+1+m, cmp);
    int ecnt=0;
    for(int i=1; i<=m; i++){
        r1 = myfind(ed[i].fro);
        r2 = myfind(ed[i].too);
        if(r1!=r2){
            fa[r2] = r1;
            add_edge(ed[i].fro, ed[i].too, ed[i].val);
            add_edge(ed[i].too, ed[i].fro, ed[i].val);
        }
    }
}
void bd(int u){
    for(int i=hea[u]; i; i=edge[i].nxt){
        int t=edge[i].too;
        if(!deep[t]){
            grand[t][0] = u;
            minq[t][0] = edge[i].val;
            deep[t] = deep[u] + 1;
            bd(t);
        }
    }
}
void inits(){
    for(int i=1; i<=13; i++)
        for(int j=1; j<=n; j++){
            grand[j][i] = grand[grand[j][i-1]][i-1];
            minq[j][i] = min(minq[j][i-1], minq[grand[j][i-1]][i-1]);
        }
}
int lca(int x, int y){
    int ans = 0x3f3f3f3f;
    if(deep[x]<deep[y]){
        x ^= y;
        y ^= x;
        x ^= y;
    }
    for(int i=13; i>=0; i--)
        if(deep[grand[x][i]]>=deep[y]){
            ans = min(ans, minq[x][i]);
            x = grand[x][i];//请注意语句次序！！！
        }
    if(x==y)    return ans;
    for(int i=13; i>=0; i--)
        if(grand[x][i]!=grand[y][i]){
            ans = min(ans, min(minq[x][i], minq[y][i]));
            x = grand[x][i];
            y = grand[y][i];//请注意语句次序！！！
        }
    return min(ans, min(minq[x][0], minq[y][0]));
}
int main(){
    init();
    kruskal();
    for(int i=1; i<=n; i++)
        if(!deep[i]){
            deep[i] = 1;
            bd(i);
        }
    inits();
    rn(q);
    for(int i=1; i<=q; i++){
        rn(uu);
        rn(vv);
        if(myfind(uu)!=myfind(vv))  printf("-1\n");
        else                        printf("%d\n", lca(uu, vv));
    }
    return 0;
}