1. 程式人生 > >貨車運輸

貨車運輸

gis 避免 show blog getch main 描述 names tro

貨車運輸https://www.luogu.org/problemnew/show/P1967

題目描述

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 。

輸出格式:

輸出共有 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。


正解:先按邊權從大到小排序,跑Kruscal(即最大生成樹),樹上倍增維護lim[i][j]表示節點i到它的2^j祖先的路徑上的限重(邊權最小值),然後每個詢問跑LCA即可.
但可能會有多棵樹,比如這組數據:

數據輸入
5 2
1 2 3
3 4 5
1
3 4

數據輸出
5

具體處理見代碼:

#define RG register
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
const int N=1e5+5,M=5e5+5;
inline int read()
{
    RG int x=0,w=1;RG char ch=getchar();
    while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘)w=-1;ch=getchar();}
    while(ch>=‘0‘&&ch<=‘9‘)x=x*10+ch-‘0‘,ch=getchar();
    return x*w;
}
int n,m,Q,cnt,ct;
int par[N],last[N],dep[N],f[N][22],lim[N][22],s[N];//dep深度,f[i][j]表示i節點2^j祖先的序號,s為根節點數組
struct node{
    int u,v,w;
}e1[M];
struct edge{//鄰接表
    int to,next,w;
}e[N*2];
void insert(int u,int v,int w)
{
    s[++ct]=u;//將連入生成樹的節點中加入根節點數組(只放u是因為v必在u所在的樹中,若u為根節點,則v必不為根節點)
    e[++cnt]=(edge){v,last[u],w};last[u]=cnt;
    e[++cnt]=(edge){u,last[v],w};last[v]=cnt;
}
bool cmp(node a,node b)
{
    return a.w>b.w;
}
int find(int x){return x==par[x]?x:par[x]=find(par[x]);}
void Kruscal()//最大生成樹
{
    int num=0;
    for(int i=1;i<=m;i++)
    {
        if(num==n-1)return;
        int x=find(e1[i].u),y=find(e1[i].v);
        if(x==y)continue;
        par[x]=find(par[y]);//合並
        insert(e1[i].u,e1[i].v,e1[i].w);//連邊
        num++;
    }
}
void dfs(int now)
{
    for(int i=last[now];i;i=e[i].next)
    {
        int v=e[i].to;
        if(dep[v])continue;
        dep[v]=dep[now]+1;
        f[v][0]=now;
        lim[v][0]=e[i].w;//lim[i][0]即節點i到父親節點的邊權
        dfs(v);
    }
}
void init()
{
    for(int i=1;i<=ct;i++)
        if(dep[s[i]]==0)//保證每棵樹只dfs一遍
        {
            dep[s[i]]=1;//根節點深度為1
            dfs(s[i]);//為了避免有多棵樹,所以都dfs
        }
    for(int j=1;j<=20;j++)
        for(int i=1;i<=n;i++)
        {
            f[i][j]=f[f[i][j-1]][j-1];
            lim[i][j]=min(lim[f[i][j-1]][j-1],lim[i][j-1]);//因為是限重,所以取最小值
        }
}
int LCA(int x,int y)//倍增求LCA
{
    int LIM=1e9;
    if(dep[x]<dep[y])swap(x,y);
    for(int i=20;i>=0;i--)
        if(f[x][i]&&dep[f[x][i]]>=dep[y])
        {
            LIM=min(LIM,lim[x][i]);
            x=f[x][i];
        }
    if(x==y)return LIM;//特判y是x祖先
    for(int i=20;i>=0;i--)
        if(f[x][i]&&f[y][i]&&f[x][i]!=f[y][i])
        {
            LIM=min(LIM,min(lim[x][i],lim[y][i]));
            x=f[x][i];
            y=f[y][i];
        }
    LIM=min(LIM,min(lim[x][0],lim[y][0]));//註意還要取一次最小值
    return LIM;
}
int main()
{
    n=read();m=read();
    for(int i=1;i<=n;i++)par[i]=i;
    for(int i=1;i<=m;i++)e1[i].u=read(),e1[i].v=read(),e1[i].w=read();
    sort(e1+1,e1+m+1,cmp);
    Kruscal();
    memset(lim,127,sizeof(lim));//lim初始化為inf
    init();
    Q=read();
    while(Q--)
    {
        int x=read(),y=read();
        if(find(x)!=find(y))printf("-1\n");//不在同一棵樹上
        else printf("%d\n",LCA(x,y));
    }
    return 0;
}

貨車運輸