1. 程式人生 > >luogu1967[NOIP2013D1T3] 貨車運輸 (最大生成樹+LCA)

luogu1967[NOIP2013D1T3] 貨車運輸 (最大生成樹+LCA)

etc 註意 bre 限制 文件名 con style ont 描述

題目描述

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。

kruskal和lca亂搞搞就好了qwq

 1 #include "bits/stdc++.h"
 2 using namespace std;
 3 typedef long long LL;
 4 const int MAX1=1e4+5,MAX2=1e5+5;
 5 int n,m,q,f[MAX1];
 6 int tot,head[MAX1],adj[MAX2<<1],wei[MAX2<<1],next[MAX2<<1
]; 7 int deep[MAX1],fa[MAX1][21],g[MAX1][21]; 8 struct Edge{int x,y,w;}edge[MAX2]; 9 bool cmp(Edge x,Edge y){return x.w>y.w;} 10 int getfather(int x){return f[x]==x?x:f[x]=getfather(f[x]);} 11 inline int read(){ 12 int an=0,x=1;char c=getchar(); 13 while (c<0 || c>9) {if (c==-) x=-1;c=getchar();} 14 while (c>=0 && c<=9) {an=(an<<3)+(an<<1)+c-0;c=getchar();} 15 return an*x; 16 } 17 void addedge(int u,int v,int w){tot++;adj[tot]=v,wei[tot]=w,next[tot]=head[u],head[u]=tot;} 18 void dfs(int x){ 19 int i,j; 20 for (i=1;i<=20;i++){ 21 if (deep[x]<(1<<i)) break; 22 fa[x][i]=fa[ fa[x][i-1] ][ i-1 ];g[x][i]=min(g[x][i-1],g[ fa[x][i-1] ][ i-1 ]); 23 } 24 for (i=head[x];i;i=next[i]) 25 if (deep[adj[i]]==-1){ 26 deep[adj[i]]=deep[x]+1,g[adj[i]][0]=wei[i],fa[adj[i]][0]=x;dfs(adj[i]); 27 } 28 } 29 int LCA(int x,int y){ 30 int an=1e9,i,j; 31 if (deep[x]==-1 || deep[y]==-1) return -1; 32 if (deep[x]<deep[y]) swap(x,y); 33 int dd=deep[x]-deep[y]; 34 for (i=20;i>=0;i--) 35 if (dd&(1<<i)) 36 an=min(an,g[x][i]),x=fa[x][i]; 37 for (i=20;i>=0;i--) 38 if (fa[x][i]!=fa[y][i]){ 39 an=min(min(an,g[x][i]),g[y][i]); 40 x=fa[x][i],y=fa[y][i]; 41 } 42 if (x!=y) an=min(g[x][0],min(g[y][0],an)); 43 return an; 44 } 45 int main(){ 46 freopen ("truck.in","r",stdin);freopen ("truck.out","w",stdout); 47 int i,j,x,y; 48 n=read(),m=read(); 49 for (i=1;i<=m;i++) edge[i].x=read(),edge[i].y=read(),edge[i].w=read(); 50 sort(edge+1,edge+m+1,cmp); 51 for (i=1;i<=n;i++) f[i]=i; 52 for (tot=i=1;i<=m;i++){ 53 int tx=getfather(edge[i].x),ty=getfather(edge[i].y); 54 if (tx!=ty){ 55 addedge(edge[i].x,edge[i].y,edge[i].w);addedge(edge[i].y,edge[i].x,edge[i].w); 56 f[tx]=ty; 57 } 58 } 59 memset(deep,-1,sizeof(deep));deep[1]=0; 60 q=read();dfs(1); 61 for (i=1;i<=q;i++){ 62 x=read(),y=read(); 63 printf("%d\n",LCA(x,y)); 64 } 65 return 0; 66 }

luogu1967[NOIP2013D1T3] 貨車運輸 (最大生成樹+LCA)