2017.10.07【NOIP提高组】模拟赛B组 Heatwave 题解

时间:2021-08-04 06:21:04

传送门

Description

  给你N个点的无向连通图,图中有M条边,第j条边的长度为: d_j.
  现在有 K个询问。
  每个询问的格式是:A B,表示询问从A点走到B点的所有路径中,最长的边最小值是多少?

Input

  第一行: N, M, K。
  第2..M+1行: 三个正整数:X, Y, and D (1 <= X <=N; 1 <= Y <= N). 表示X与Y之间有一条长度为D的边。
  第M+2..M+K+1行: 每行两个整数A B,表示询问从A点走到B点的所有路径中,最长的边最小值是多少?

Output

  对每个询问,输出最长的边最小值是多少。

Sample Input

6 6 8
1 2 5
2 3 4
3 4 3
1 4 8
2 5 7
4 6 2
1 2
1 3
1 4
2 3
2 4
5 1
6 2
6 1

Sample Output

5
5
5
4
4
7
4
5

Data Constraint

50% 1<=N,M<=3000 其中30% K<=5000
​100% 1 <= N <= 15,000 1 <= M <= 30,000 1 <= d_j <= 1,000,000,000 1 <= K <= 20,000

Analysis

一眼看过去就觉得是SPFA,然而,我真不知道那些用spfa的人是怎么过的
最后想了想,发现既然是求最长边权的最小值,那就用最小生成树呗
先过一遍最小生成树,然后求A-B的最短路径,怎么求呢???
其实很简单,既然都是棵树了,那么求LCA就好啦,边求的时候顺便求出边权的最大值就是答案了
相通了就觉得好水……

Realization

#include<bits/stdc++.h>
using namespace std;
#define fo(i,a,b) for(i=a;i<=b;i++)
#define fd(i,a,b) for(i=a;i>=b;i--)
#define fin(x) freopen(""#x".in","r",stdin)
#define fout(x) freopen(""#x".out","w",stdout)
#define open(x) fin(x);fout(x)
#define ll long long
#define inf INT_MAX

struct M{
 int u,v,w;
}a[30001];
int n,m,k,gx,gy,x,y,len,i,j,fa[15001],d[15001],h[80001],ft[80001],next[80001],last[80001],g[15001][17],tw[15001][17];
bool cmp(M a,M h){
 return a.w<h.w;
}
int get(int k){
 if(fa[k]==0) return k;
 return (fa[k]=get(fa[k]));
}
void add(int x,int y,int z){
 h[++len]=y;
 next[len]=last[x];
 last[x]=len;
 ft[len]=z;
}     
void dg(int k,int y){
 int p=last[k];
 while(p!=0){
 if(h[p]!=y){
 g[h[p]][0]=k;
 d[h[p]]=d[k]+1;
 tw[h[p]][0]=ft[p];
 dg(h[p],k);
 }
 p=next[p];
 }
}
int lca(int x,int y){
 int k,ans=0;
 if(d[x]<d[y]) swap(x,y);
 k=trunc(log(d[x]-d[y]+1)/log(2));
 while(k>=0){
 if(d[g[x][k]]>d[y]) ans=max(ans,tw[x][k]),x=g[x][k];
 k--;
 }
 if(d[x]!=d[y]) ans=max(ans,tw[x][0]),x=g[x][0];
 k=trunc(log(d[x])/log(2));
 while(k>=0){
 if(g[x][k]!=g[y][k]) ans=max(ans,max(tw[x][k],tw[y][k])),x=g[x][k],y=g[y][k];
 k--;
 }
 if(x==y) return ans;
 else return max(ans,max(tw[x][0],tw[y][0]));
}
int main(){
 scanf("%d%d%d",&n,&m,&k);
 fo(i,1,m) scanf("%d%d%d",&a[i].u,&a[i].v,&a[i].w);
 sort(a+1,a+m+1,cmp);
 fo(i,1,m){
 gx=get(a[i].u);
 gy=get(a[i].v);
 if(gx!=gy){
 add(a[i].u,a[i].v,a[i].w);
 add(a[i].v,a[i].u,a[i].w);
 fa[gx]=gy;
 }
 }
 d[1]=1;
 dg(1,0);
 fo(j,1,trunc(log(n)/log(2)))
 fo(i,1,n){
 g[i][j]=g[g[i][j-1]][j-1];
 tw[i][j]=max(tw[i][j-1],tw[g[i][j-1]][j-1]);
 }
 fo(i,1,k){
 scanf("%d%d",&x,&y);
 printf("%d\n",lca(x,y));
 }
}