题目描述
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。
输入输出样例
4 3
1 2 4
2 3 3
3 1 1
3
1 3
1 4
1 3
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。
因为这题在bzoj里又是那啥啥，所以就又交luogu了
代码这么长又比别人慢真是太尴尬了
这么慢思路肯定很easy对吧
->先求用kruskal最大生成树，可忽略重边问题
->以重心为根造树
->对于每一组询问(ss,tt)，求其lca
->我这边模拟了它从ss和tt走到lca，程序秒速变慢 QwQ(懒癌晚期)
#include<cstdio>
#include<cstring>
#include<iostream>
#include<vector>
#include<algorithm>
#define inf 0x3f3f3f3f
using namespace std;
#define ss first
#define tt second inline int read(){
char ch;
int re=;
bool flag=;
while((ch=getchar())!='-'&&(ch<''||ch>''));
ch=='-'?flag=:re=ch-'';
while((ch=getchar())>=''&&ch<='') re=re*+ch-'';
return flag?-re:re;
} struct odge{
int from,to,w;
odge(int from=,int to=,int w=):
from(from),to(to),w(w){}
bool operator < (const odge &n1) const {
return w>n1.w;
}
};
struct edge{
int to,w,next;
edge(int to=,int w=,int next=):
to(to),w(w),next(next){}
};
const int maxn=,maxm=;
const int maxlogn=;
int n,m;
vector<odge> odges;
vector<edge> edges;
int head[maxn];
int par[maxn];
int cnt=;
int F[maxn],son[maxn];
bool vis[maxn];
int sum,root;
int depth[maxn],parent[maxlogn][maxn];
int q;
typedef pair<int,int> PII;
vector<PII> ques;
int length[maxn]; inline void add_edge(odge e){
edges.push_back(edge(e.to,e.w,head[e.from]));
head[e.from]=++cnt;
edges.push_back(edge(e.from,e.w,head[e.to]));
head[e.to]=++cnt;
} int find(int x){
return par[x]==x?x:par[x]=find(par[x]);
} void kruskal(){
memset(head,,sizeof head);
edges.push_back(edge(,,));
int cnt1=,cnt2=;
while(cnt2<m){
odge e=odges[cnt2++];
int dx=find(e.from),dy=find(e.to);
if(dx==dy) continue;
add_edge(e);
par[dy]=dx;
if(++cnt1==n) break;
}
} void getroot(int x,int fa){
son[x]=; F[x]=;
for(int ee=head[x];ee;ee=edges[ee].next)
if(edges[ee].to!=fa){
getroot(edges[ee].to,x);
son[x]+=son[edges[ee].to];
F[x]=max(F[x],son[edges[ee].to]);
}
F[x]=max(F[x],sum-son[x]);
if(F[x]<F[root]) root=x;
} void dfs(int x,int fa,int deep){
parent[][x]=fa;
depth[x]=deep;
for(int ee=head[x];ee;ee=edges[ee].next)
if(edges[ee].to!=fa){
length[edges[ee].to]=edges[ee].w;
dfs(edges[ee].to,x,deep+);
}
} void init(){
n=read(); m=read();
for(int i=;i<m;i++){
int from=read(),to=read(),w=read();
odges.push_back(odge(from,to,w));
}
for(int i=;i<=n;i++)
par[i]=i;
sort(odges.begin(),odges.end());
kruskal();
sum=F[root=]=n;
getroot(,);
dfs(root,-,);
for(int k=;k+<maxlogn;k++)
for(int v=;v<=n;v++)
if(parent[k][v]<) parent[k+][v]=-;
else parent[k+][v]=parent[k][parent[k][v]];
q=read();
for(int i=;i<q;i++){
int from=read(),to=read();
ques.push_back(make_pair(from,to));
}
} int getlca(int u,int v){
if(depth[u]>depth[v]) swap(u,v);
for(int k=;k<maxlogn;k++)
if(depth[v]-depth[u]>>k&) v=parent[k][v];
if(u==v) return u;
for(int k=maxlogn-;k>=;k--)
if(parent[k][u]!=parent[k][v]){
u=parent[k][u];
v=parent[k][v];
}
return parent[][u];
} void solve(){
for(int i=;i<q;i++){
int u=ques[i].ss,v=ques[i].tt;
if(find(u)!=find(v)){
puts("-1");
continue;
}
int ans=inf;
int lca=getlca(ques[i].ss,ques[i].tt); while(u!=lca){
ans=min(ans,length[u]);
u=parent[][u];
}
while(v!=lca){
ans=min(ans,length[v]);
v=parent[][v];
} printf("%d\n",ans);
}
} int main(){
//freopen("temp.in","r",stdin);
init();
solve();
return ;
}
天亮的时候 你还没回来
亲爱的 我们的孩子醒来了