[LUOGU] P1967 货车运输

题目描述
A 国有 n 座城市，编号从 1 到 n，城市之间有 m 条双向道路。每一条道路对车辆都有重量限制，简称限重。现在有 q 辆货车在运输货物， 司机们想知道每辆车在不超过车辆限重的情况下，最多能运多重的货物。

输入输出格式
输入格式：
输入文件名为 truck.in。

输入文件第一行有两个用一个空格隔开的整数 n，m，表示 A 国有 n 座城市和 m 条道

路。 接下来 m 行每行 3 个整数 x、 y、 z，每两个整数之间用一个空格隔开，表示从 x 号城市到 y 号城市有一条限重为 z 的道路。注意： x 不等于 y，两座城市之间可能有多条道路 。

接下来一行有一个整数 q，表示有 q 辆货车需要运货。

接下来 q 行，每行两个整数 x、y，之间用一个空格隔开，表示一辆货车需要从 x 城市运输货物到 y 城市，注意： x 不等于 y 。

输出格式：
输出文件名为 truck.out。

输出共有 q 行，每行一个整数，表示对于每一辆货车，它的最大载重是多少。如果货

车不能到达目的地，输出-1。

输入输出样例
输入样例#1： 
4 3
1 2 4
2 3 3
3 1 1
3
1 3
1 4
1 3
输出样例#1： 
3
-1
3
说明
对于 30%的数据，0 < n < 1,000，0 < m < 10,000，0 < q< 1,000；

对于 60%的数据，0 < n < 1,000，0 < m < 50,000，0 < q< 1,000；

对于 100%的数据，0 < n < 10,000，0 < m < 50,000，0 < q< 30,000，0 ≤ z ≤ 100,000。

最大生成树+LCA+树上倍增
求最小值的最大值，和那个星际航行类似
求出最大生成树，在线处理LCA，同时用数组g倍增统计最小值

有这么几个脑残错误
1.最大生成树，建树的时候用的查找过的并查集代表建树
2.LCA忘记判断x和y在同一深度后是否为同一元素（即LCA），并且在调整x的深度时就要沿途更新
3.倍增范围 1-n（？）

//Stay foolish,stay hungry,stay young,stay simple
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cmath>
using namespace std;

const int MAXN=60005;

int n,m,q;

inline int read_d(){
    int s=0;
    char c;
    while(c=getchar(),c<'0'||c>'9');
    while(c<='9'&&c>='0'){
        s=s*10+c-'0';
        c=getchar();
    }
    return s;
}

struct Edge{
    int next,to,w;
}e[MAXN<<2];
int ecnt,head[MAXN];
inline void add(int x,int y,int w){
    e[++ecnt].next = head[x];
    e[ecnt].to = y;
    e[ecnt].w = w;
    head[x]=ecnt;
}

int uni[MAXN>>2];
int fnd_uni(int x){
    return x==uni[x]?x:uni[x]=fnd_uni(uni[x]);
}
void cat_uni(int x,int y){
    x=fnd_uni(x);
    y=fnd_uni(y);
    uni[y]=uni[x];
}

struct oEage{
    int from,next,to,w;
}oe[MAXN<<2];
bool cmp(const oEage &x,const oEage &y){
    return x.w > y.w ;
}
int oecnt,ohead[MAXN>>2];
inline void add_o(int x,int y,int w){
    oe[++oecnt].from = x;
    oe[oecnt].next = ohead[x];
    oe[oecnt].to = y;
    oe[oecnt].w = w;
    ohead[x]=oecnt;
}
int fa[MAXN>>2];
int fnd(int x){
    return x==fa[x]?x:fa[x]=fnd(fa[x]);
}
void cat(int x,int y){
    x=fnd(x);y=fnd(y);
    fa[y]=x;
}
void kls(){//longest
    int cnt=0;
    for(int i=1;i<=oecnt;i++){
        oEage v=oe[i];
        int x=fnd(v.from),y=fnd(v.to);
        if(x==y) continue;
        cnt++;
        add(v.from,v.to,v.w);//!! 
        add(v.to,v.from,v.w);
        cat(x,y);
        if(cnt==n-1) break;
    }
}

int dep[MAXN];
int f[MAXN>>2][15];
int g[MAXN>>2][15];
void dfs(int now,int pre){
    for(int i=head[now];i;i=e[i].next){
        int v=e[i].to;
        if(dep[v]) continue;
        if(v==pre) continue;
        dep[v]=dep[now]+1;
//      f[i][0]=now;
//      g[i][0]=e[i].w;
        f[v][0]=now;
        g[v][0]=e[i].w;
        dfs(v,now);
    }
}

void ReDouble(){
    for(int i=1;i<=n;i++){
        for(int j=1;j<=14;j++){
            g[i][j]=1<<28;
        }
    }
    for(int j=1;(1<<j)<=n;j++){
        for(int i=1;i<=n;i++){//?!
            f[i][j]=f[f[i][j-1]][j-1];
            g[i][j]=min(g[i][j-1],g[f[i][j-1]][j-1]);
        }
    }
}

int solve(int x,int y){
    if(dep[x]<dep[y]) swap(x,y);
    int ret=1<<28;
    for(int j=14;j>=0;j--){
        if(dep[f[x][j]]<dep[y]) continue;
        //!
        ret=min(ret,g[x][j]);
        x=f[x][j];
    } 
    if(x==y) return ret;//
    for(int j=log2(n);j>=0;j--){
        if(f[x][j]==f[y][j]) continue;
        ret=min(ret,min(g[x][j],g[y][j]));
        x=f[x][j];y=f[y][j];
    }
    ret=min(ret,g[x][0]);
    ret=min(ret,g[y][0]);
//  cout<<"LCA:"<<x<<" "<<y<<" "<<f[x][0]<<endl;
    return ret;
}

int main(){
    n=read_d();
    m=read_d();
    int x,y,w;
    for(int i=1;i<=n;i++){
        uni[i]=i;
        fa[i]=i;
    }
    for(int i=1;i<=m;i++){
        x=read_d();
        y=read_d();
        w=read_d();
        add_o(x,y,w);
        add_o(y,x,w);
        cat_uni(x,y);
    }
    sort(oe+1,oe+1+oecnt,cmp);
    kls();
    for(int i=1;i<=n;i++){
        if(!dep[i]) dep[i]=1,dfs(i,0);
    }
    ReDouble();
    q=read_d();
    for(int i=1;i<=q;i++){
        x=read_d();
        y=read_d();
        if(fnd_uni(x)!=fnd_uni(y)){
            puts("-1");
            continue;
        } 
        printf("%d
",solve(x,y));
    }
    return 0;
}