• [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
    转载请注明出处 -- 如有意见欢迎评论
  • 相关阅读:
    全文检索(SOLR)前端应用浅析续 Rails前端分析
    全文搜索应用 企业搜索SearchBlox
    持续集成(CI) 几种测试的区别(摘录)
    图解全文检索SOLR的酷应用[Ajax Solr]
    Php如何实现下载功能超详细流程分析
    在MySQL字段中使用逗号分隔符
    session的垃圾回收机制
    【转】apache常用配置
    深入理解PHP之数组(遍历顺序)
    正则表达式的子模式详解
  • 原文地址:https://www.cnblogs.com/Chorolop/p/7347538.html
Copyright © 2020-2023  润新知