洛谷 P1967 货车运输(lca+最小生成树) 题解

94 篇文章 0 订阅

题目来源:

https://www.luogu.org/problemnew/show/P1967

题目描述:

题目描述

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

输入输出格式

输入格式:

 

第一行有两个用一个空格隔开的整数n,mn,m,表示 AA 国有nn 座城市和 mm 条道路。

接下来 mm行每行33个整数 x, y, zx,y,z,每两个整数之间用一个空格隔开,表示从 xx号城市到yy号城市有一条限重为 zz 的道路。注意: xx 不等于 yy,两座城市之间可能有多条道路 

接下来一行有一个整数 q,表示有 q 辆货车需要运货。

接下来 q 行,每行两个整数 x、y,之间用一个空格隔开,表示一辆货车需要从 x 城市运输货物到 y 城市,注意:x 不等于 y 

 

输出格式:

 

共有 qq 行,每行一个整数,表示对于每一辆货车,它的最大载重是多少。如果货车不能到达目的地,输出-1−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\%30%的数据,0 < n < 1,000,0 < m < 10,000,0 < q< 1,0000<n<1,000,0<m<10,000,0<q<1,000;

对于 60\%60%的数据,0 < n < 1,000,0 < m < 50,000,0 < q< 1,0000<n<1,000,0<m<50,000,0<q<1,000;

对于 100\%100%的数据,0 < n < 10,000,0 < m < 50,000,0 < q< 30,000,0 ≤ z ≤ 100,0000<n<10,000,0<m<50,000,0<q<30,000,0≤z≤100,000。

解题思路:

       我本来想,这题应该还好吧,可以解出来这题的,结果竟然和我想象的完全不一样,需要这样解,我想需要好好消化一下。

首先我们可以先分析出,要先求出图的最大生成树,因为如果边权值不在最大生成树里,无论怎样,车都不会往那条路走,然后对于每个询问a,b,就是求a到b的最小边权值,因为在树里,a到b的路径只有一条,所以我们可以用倍增来记录一些数据,然后通过求lca的方式求出a到b路径中的最小值。(倍增也是现学的%%%),dfs里就是处理父子节点信息的。具体细节看代码,就行了,。。。这题没完全懂,特别是lca那段,以后还得在研究研究。。。

代码:

#include <iostream>
#include <iomanip>
#include <cstdio>
#include <cstring>
#include <string>
#include <vector>
#include <map>
#include <algorithm>
#include <queue>
#include <stack>
#include <cmath>
#define inf 0x3f3f3f3f
#define ll long long
using namespace std;
const int  maxn=1e6+10;
int n,m,fa[maxn],cnt=0,deep[maxn],head[maxn],f[maxn][21],w[maxn][21],vis[maxn];
struct newt
{
    int from,to,cost;
}edge1[maxn];
struct newtt
{
    int to,next,cost;
}e[maxn];
void addedge(int u,int v,int w)
{
    e[cnt].to=v;
    e[cnt].next=head[u];
    e[cnt].cost=w;
    head[u]=cnt++;
}
bool cmp(newt a,newt b)
{
    return a.cost>b.cost;
}
int fi(int x)
{
    if(x==fa[x])return x;
    return fa[x]=fi(fa[x]);
}
void kruscal()
{
    for(int i=1;i<=m;i++)
    {
        int f1=fi(edge1[i].from),f2=fi(edge1[i].to);
        if(f1!=f2)
        {
            addedge(edge1[i].from,edge1[i].to,edge1[i].cost);
            addedge(edge1[i].to,edge1[i].from,edge1[i].cost);
            fa[f1]=f2;
        }
    }
}
void dfs(int u)
{
    vis[u]=1;
    for(int i=head[u];i!=-1;i=e[i].next)
    {
        int v=e[i].to;
        if(vis[v])continue;
        deep[v]=deep[u]+1;
        f[v][0]=u;
        w[v][0]=e[i].cost;
        dfs(v);
    }
}
int lca(int a,int b)
{
    if(fi(a)!=fi(b))return -1;
    if(deep[a]>deep[b])swap(a,b);
    int ans=inf;
    for(int i=20;i>=0;i--)
    if(deep[f[b][i]]>=deep[a])
    {
        ans=min(ans,w[b][i]);
        b=f[b][i];
    }
    if(a==b)return ans;
    for(int i=20;i>=0;i--)
    {
        if(f[a][i]!=f[b][i])
        {
            ans=min(ans,min(w[a][i],w[b][i]));
            a=f[a][i];
            b=f[b][i];
        }
    }
    ans=min(ans,min(w[a][0],w[b][0]));
    return ans;
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)fa[i]=i,head[i]=-1;
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d%d",&edge1[i].from,&edge1[i].to,&edge1[i].cost);
    }
    sort(edge1+1,edge1+m+1,cmp);
    kruscal();
    for(int i=1;i<=n;i++)if(!vis[i])
    {
        deep[i]=1;
        dfs(i);
        f[i][0]=i;
        w[i][0]=inf;
    }
    for(int i=1;i<=20;i++)
    for(int j=1;j<=n;j++)
    {
        f[j][i]=f[f[j][i-1]][i-1];
        w[j][i]=min(w[f[j][i-1]][i-1],w[j][i-1]);
    }
    int q;
    scanf("%d",&q);
    for(int i=1,a,b;i<=q;i++)
    {
        scanf("%d%d",&a,&b);
        printf("%d\n",lca(a,b));
    }
    return 0;
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值