NOIP2013货车运输

NOIP2013货车运输

---

题目描述

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

输入

第一行有两个用一个空格隔开的整数n，m，表示A国有n座城市和m条道路。接下来m行每行3个整数x、y、z，每两个整数之间用一个空格隔开，表示从x号城市到y号城市有一条限重为z的道路。注意：x不等于y，两座城市之间可能有多条道路。接下来一行有一个整数q，表示有q辆货车需要运货。接下来q行，每行两个整数x、y，之间用一个空格隔开，表示一辆货车需要从x城市运输货物到y城市，注意：x不等于y。
数据范围：n最大10000，m最大50000，q为30000，z为100000.

输出

共有q行，每行一个整数，表示对于每一辆货车，它的最大载重是多少。如果货车不能到达目的地，输出-1。

样例输入

4 3
1 2 4
2 3 3
3 1 1
3
1 3
1 4
1 3

样例输出

3
-1
3

做法

两点之间尽可能大的最小承重路径一定在最大生成树上，于是问题转化为查询最大生成树上，两点之间的最小边的查询，可以倍增，也可以树剖。注意还有可能不能互相到达。

#include <bits/stdc++.h>
#define rep(i,a,b) for(int i=a;i<=b;++i)
typedef long long ll;
const int N = 10010;
const int M = 50100;
const int inf = 0x3f3f3f3f;
inline int read() {char c=getchar(); int x=0,f=1;while(!isdigit(c)){if(c=='-')f=-1;c=getchar();}while(isdigit(c)){x=x*10+c-'0';c=getchar();}return x*f;
}
using namespace std;
struct edge{int u,v,nxt,w;}E[M<<1],E0[M<<1];
int h[N],cc;
void add(int u,int v,int z) {E[cc].v=v;E[cc].w=z;E[cc].nxt=h[u];h[u]=cc;++cc;
}
bool cmp(edge a,edge b) {return a.w > b.w;}
int n,m;
int f[N][22],mn[N][22],fa[N];
int find(int x) {if(x==fa[x])return x;return fa[x]=find(fa[x]);
}
void merge(int x,int y) {x=find(x), y =find(y);if(x!=y) fa[x]=y;
}
void krus() {sort(E0+1,E0+1+m,cmp);rep(i,1,m) {int tx = find(E0[i].u), ty = find(E0[i].v);if(tx!=ty) {merge(tx,ty);add(E0[i].u,E0[i].v,E0[i].w);add(E0[i].v,E0[i].u,E0[i].w);}}
}
bool vis[N];
int dep[N];
void dfs(int u,int pre) {vis[u] = 1;f[u][0] = pre;for(int i=h[u];~i;i=E[i].nxt) {int v = E[i].v;if(!vis[v]&&v!=pre) {dep[v]=dep[u]+1;mn[v][0] = E[i].w;dfs(v,u);for(int j=1;j<=20;++j) f[v][j] = f[f[v][j-1]][j-1];for(int j=1;j<=20;++j) mn[v][j] = min(mn[v][j-1],mn[f[v][j-1]][j-1]);}}
}
void init() {memset(mn,inf,sizeof(mn));krus();rep(i,1,n)if(!vis[i])dfs(i,0);
//    rep(i,1,n) {
//        printf("%d :\n",i);
//        printf("dep[%d] = %d\n",i,dep[i]);
//        rep(j,0,3) printf("fa[%d][%d] = %d\n",i,j,f[i][j]);
//        rep(j,0,3) printf("mn[%d][%d] = %d\n",i,j,mn[i][j]);
//    }
}
int lca(int x,int y) {if(dep[x] > dep[y]) swap(x,y);int tx=x,ty=y;for(int d=dep[y]-dep[x],i=0;d;d>>=1,++i) if(d&1) ty = f[ty][i];if(tx == ty) return tx;for(int i=20;i>=0;--i) {if(f[tx][i]==f[ty][i]) continue;tx = f[tx][i], ty = f[ty][i];}return f[tx][0];
}
int ask_mn(int x,int p) {int ans = inf, tx=x;for(int d=dep[x]-dep[p],i=0;d;d>>=1,++i) if(d&1) {ans = min(ans,mn[tx][i]);tx = f[tx][i];}return ans;
}
int ask(int x,int y) {if(find(x)!=find(y)) return -1;int LCA = lca(x,y);return min(ask_mn(x,LCA),ask_mn(y,LCA));
}
int main() {n = read(), m =read();rep(i,1,n) h[i]=-1,fa[i]=i;rep(i,1,m) E0[i].u=read(),E0[i].v=read(),E0[i].w=read();init();int q = read();while(q--) {int x=read(), y =read();printf("%d\n",ask(x,y));}return 0;
}
//9 8
//1 2 4
//1 5 3
//2 3 3
//2 4 2
//3 9 2
//3 8 7
//5 7 6
//5 6 7