[Codevs] 3287 货车运输

3287 货车运输

2013年NOIP全国联赛提高组

 时间限制: 1 s
 空间限制: 128000 KB
 题目等级 : 钻石 Diamond
 
题目描述 Description

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

输入描述 Input Description

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

输出描述 Output Description

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

样例输入 Sample Input

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

样例输出 Sample Output

3
-1
3

数据范围及提示 Data Size & Hint

对于 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。

分析 Analysis

= =从夏令营A结束的时候说要补

= =结果到夏令营B快结束时才补上

理论基础:树上操作 >> LCA / 树上倍增 + 最小生成树Kruskal

其实是道水题= =

既然边权是限重,我们可以构建出一棵最大生成树。

这样虽然绕了路,不过最低限重就可以保证不会太小了。

首先一波Kruskal选出最大生成树的边并存进新边集(原来的图是基本要拆掉的,所以也就无所谓用不用邻接表了,直接散存)

根据新边集DFS建一棵树,同样,因为某些原因,树根节点我们随便定,拿D6 Roll一次或者干脆选1都是不错的选择。

之后维护两个倍增组,一组是最近公共祖先钦定的父节点组,另一组储存最低限重。

然后在求LCA的过程中顺便求路径上的最低限重即可。

mincost即最低限重。

代码 Code

  1 #include<cstdio>
  2 #include<iostream>
  3 #include<algorithm>
  4 #define maxn_m 100000
  5 #define maxn 50000
  6 using namespace std;
  7 
  8 struct edge1{
  9     int u,v,len;
 10 }e1[maxn_m];
 11 
 12 struct edge{
 13     int from,v,len;
 14 }e[maxn_m];
 15 
 16 int tot,first[maxn],n,m,q,a,b;
 17 void insert(int u,int v,int len){
 18     tot++;
 19     e[tot].from = first[u];
 20     e[tot].v = v;
 21     e[tot].len = len;
 22     first[u] = tot;
 23 }
 24 
 25 bool cmp(const edge1 &a,const edge1 &b){return a.len>b.len;}
 26 int pre[maxn];
 27 int find(int x){if(!pre[x]){return x;}else{pre[x] = find(pre[x]);return pre[x];}}
 28 void unite(int x,int y){if(find(x)!=find(y))pre[find(x)]=y;}
 29 void kru(){
 30     sort(e1,e1+m,cmp);
 31     
 32     for(int i = 0;i < m;i++){
 33         if(find(e1[i].u) != find(e1[i].v)){
 34             insert(e1[i].u,e1[i].v,e1[i].len);
 35             insert(e1[i].v,e1[i].u,e1[i].len);
 36             unite(e1[i].u,e1[i].v);
 37         }
 38     }
 39 }
 40 
 41 int depth[maxn],fa[maxn][20],mincost[maxn][20];
 42 bool book[maxn];
 43 void dfs(int root){
 44     for(int i = first[root];i;i = e[i].from){
 45         int v = e[i].v;
 46         if(!book[v]){
 47             book[v] = true;
 48             depth[v] = depth[root]+1;
 49             fa[v][0] = root;
 50             mincost[v][0] = e[i].len;
 51             for(int i = 1;i <= 18;i++){
 52                 fa[v][i] = fa[fa[v][i-1]][i-1];
 53                 mincost[v][i] = min(mincost[fa[v][i-1]][i-1],mincost[v][i-1]);
 54             }
 55             dfs(v);
 56         }
 57     }
 58 }
 59 
 60 int lca(int u,int v){
 61     if(depth[u] < depth[v]) swap(u,v);
 62     
 63     int maxx = 999999999;
 64     for(int i = 18;i >= 0;i--){
 65         if(depth[fa[u][i]] >= depth[v]){
 66             if(mincost[u][i]) maxx = min(maxx,mincost[u][i]);
 67             u = fa[u][i];
 68         }
 69     }
 70     
 71     if(u == v) return maxx;
 72     
 73     for(int i = 18;i >= 0;i--){
 74         if(fa[u][i] != fa[v][i]){
 75             if(mincost[u][i]) maxx = min(maxx,mincost[u][i]);
 76             if(mincost[v][i]) maxx = min(maxx,mincost[v][i]);
 77             u = fa[u][i];
 78             v = fa[v][i];
 79         }
 80     }
 81     
 82     if(mincost[u][0]) maxx = min(mincost[u][0],maxx);
 83     if(mincost[v][0]) maxx = min(mincost[v][0],maxx);
 84     
 85     return maxx;
 86 }
 87 
 88 int main(){
 89     
 90     scanf("%d%d",&n,&m);
 91     
 92     for(int i = 0;i < m;i++){
 93         scanf("%d%d%d",&e1[i].u,&e1[i].v,&e1[i].len);
 94     }
 95     
 96     kru();
 97     
 98     book[1] = true;
 99     depth[1] = 1;
100     dfs(1);
101     
102     scanf("%d",&q);
103     
104     for(int i = 0;i < q;i++){
105         scanf("%d%d",&a,&b);
106         if(find(a) != find(b)) printf("-1
");
107         else printf("%d
",lca(a,b));
108     }
109     
110     return 0;
111 }
暴力下楼擦伤手腕qwq
转载请注明出处 -- 如有意见欢迎评论
原文地址:https://www.cnblogs.com/Chorolop/p/7347538.html