NOIP提高组 2013货车运输

时间:2023-07-23 22:29:44
NOIP提高组 2013货车运输

觉得题目水的离开
不屑的大佬请离开
不会图论的请离开
…….
感谢您贡献的访问量

————————————华丽的分割线————————————

题面:

题目描述
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做。
首先题目要求尽可能多的运输货物,运输货物只与两个点之间任意一条路径的最小边有关,贪心的想一下,让选择的路径最小边最大不就行了吗?只要对图做一遍最大生成树形成最大生成树的森林,然后用一个minn数组在倍增处理fa的时候倍增处理一下边的最小值,最后当询问的两个点不属于一棵树输出-1,否则做一遍LCA,在minn数组中取出最小边的值即可。
常数巨大的代码:

# include <math.h>
# include <stdio.h>
# include <stdlib.h>
# include <iostream>
# include <string.h>
# include <algorithm>
using namespace std; # define N 100001
# define IL inline
# define RG register
# define UN unsigned
# define ll long long
# define oo 2147483647
# define mem(a, b) memset(a, b, sizeof(a))
# define max(a, b) ((a) > (b)) ? (a) : (b)
# define min(a, b) ((a) < (b)) ? (a) : (b) IL int Get(){
RG char c = '!'; RG int num = 0, z = 1;
while(c != '-' && (c > '9' || c < '0')) c = getchar();
if(c == '-') z = -1, c = getchar();
while(c >= '0' && c <= '9') num = num * 10 + c - '0', c = getchar();
return num * z;
} int n, m, ft[N], nt[N], to[N], f[N], cnt;
int Fa[N], deep[N], fa[N][20], minn[N][20];
struct Edge{
int u, v, f;
IL bool operator < (Edge b) const{
return f > b.f;
}
} edge[N]; IL int Find(RG int x){
return (x == Fa[x]) ? x : Fa[x] = Find(Fa[x]);
} IL void Add(RG int u, RG int v, RG int ff){
nt[cnt] = ft[u]; to[cnt] = v; f[cnt] = ff; ft[u] = cnt++;
} IL void Dfs(RG int u){
for(RG int i = ft[u]; i != -1; i = nt[i]){
RG int v = to[i];
if(!deep[v]){
deep[v] = deep[u] + 1;
fa[v][0] = u; minn[v][0] = f[i];
Dfs(v);
}
}
} IL int Lca(RG int u, RG int v){
RG int ans = oo;
if(deep[v] < deep[u]) swap(u, v);
for(RG int j = 14; j >= 0; j--)
if(fa[v][j] && deep[fa[v][j]] >= deep[u]){
ans = min(ans, minn[v][j]);
v = fa[v][j];
}
if(u == v) return ans;
for(RG int j = 14; j >= 0; j--)
if(fa[u][j] != fa[v][j]){
ans = min(ans, minn[u][j]);
ans = min(ans, minn[v][j]);
u = fa[u][j];
v = fa[v][j];
}
ans = min(ans, minn[u][0]);
ans = min(ans, minn[v][0]);
return ans;
} int main(){
mem(ft, -1);
n = Get(); m = Get();
for(RG int i = 1; i <= n; i++) Fa[i] = i;
for(RG int i = 1; i <= m; i++)
edge[i] = (Edge){Get(), Get(), Get()};
sort(edge + 1, edge + m + 1);
RG int t = 0, i = 0;
while(i <= m){
i++;
RG int u = Find(edge[i].u), v = Find(edge[i].v);
if(u != v){
t++;
Fa[u] = Fa[v];
Add(edge[i].u, edge[i].v, edge[i].f);
Add(edge[i].v, edge[i].u, edge[i].f);
}
if(t == n - 1) break;
}
for(RG int i = 1; i <= n; i++)
if(!deep[i]) deep[i] = 1, Dfs(i);
for(RG int j = 1; (1 << j) <= n; j++)
for(RG int i = 1; i <= n; i++){
minn[i][j] = min(minn[i][j - 1], minn[fa[i][j - 1]][j - 1]);
fa[i][j] = fa[fa[i][j - 1]][j - 1];
}
RG int q = Get();
while(q--){
RG int u = Get(), v = Get();
if(Find(u) != Find(v)) printf("-1\n");
else printf("%d\n", Lca(u, v));
}
return 0;
}

龟速285ms过