【BZOJ3732】Network

Description

给你N个点的无向图 (1 <= N <= 15,000)，记为：1…N。
图中有M条边 (1 <= M <= 30,000) ，第j条边的长度为： d_j ( 1 < = d_j < = 1,000,000,000).

现在有 K个询问 (1 < = K < = 15,000)。
每个询问的格式是：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
HINT

1 <= N <= 15,000

1 <= M <= 30,000

1 <= d_j <= 1,000,000,000

1 <= K <= 15,000

Source

看见QDEZ Oxer大爷写这个题我就去膜一下
嗯…一遍MST然后链剖
剖毛线(╯‵□′)╯︵┻━┻
这TM就是NOIP货车运输啊QwQ

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define MAXN 30100
#define MAXINT 0x7fffffff
using namespace std;
int n,m,k,top,Top;
int f[MAXN];
int deep[MAXN],fa[MAXN][20],dis[MAXN][20],size[MAXN],chain[MAXN],num[MAXN];
bool vis[MAXN];
int find(int x)
{
    return f[x]==x?f[x]:f[x]=find(f[x]);
}
void Union(int a,int b)
{
    int x=find(a),y=find(b);
    if (x!=y)   f[y]=x;
}
struct Edge
{
    int u,v,w;
    bool operator <(const Edge &a)const
    {
        return w<a.w;
    }
}s[MAXN];
struct edge
{
    int to,w;
    edge *next;
}e[MAXN<<1],*prev[MAXN];
void insert(int u,int v,int w)
{
    e[++top].to=v;e[top].next=prev[u];prev[u]=&e[top];e[top].w=w;
}
void in(int &x)
{
    char ch=getchar();x=0;
    while (!(ch>='0'&&ch<='9')) ch=getchar();
    while (ch>='0'&&ch<='9')    x=x*10+ch-'0',ch=getchar();
}
void dfs(int x)
{
    vis[x]=1;
    for (int i=1;i<=19;i++)
    {
        if (deep[x]<(1<<i)) break;
        fa[x][i]=fa[fa[x][i-1]][i-1];
        dis[x][i]=max(dis[x][i-1],dis[fa[x][i-1]][i-1]);
    }
    for (edge *i=prev[x];i;i=i->next)
        if (!vis[i->to])
        {
            fa[i->to][0]=x;deep[i->to]=deep[x]+1;dis[i->to][0]=i->w;
            dfs(i->to);
        }
}
int lca(int x,int y)
{
    if (deep[x]<deep[y])    swap(x,y);
    int t=deep[x]-deep[y];
    for (int i=0;i<=19;i++)
        if ((1<<i)&t)   x=fa[x][i];
    for (int i=16;i>=0;i--)
        if (fa[x][i]!=fa[y][i]) x=fa[x][i],y=fa[y][i];
    if (x==y)   return x;
    return fa[x][0];
}
int query(int x,int y)
{
    int ret=-MAXINT,t=deep[x]-deep[y];
    for (int i=0;i<=19;i++)
        if (t&(1<<i))   ret=max(ret,dis[x][i]),x=fa[x][i];
    return ret;
}
int main()
{
    in(n);in(m);in(k);
    for (int i=1;i<=n;i++)
        for (int j=0;j<=19;j++) dis[i][j]=-MAXINT;
    for (int i=1;i<=m;i++)  in(s[i].u),in(s[i].v),in(s[i].w);
    sort(s+1,s+m+1);
    for (int i=1;i<=n;i++)  f[i]=i;
    for (int i=1;i<=m;i++)
        if (find(s[i].u)!=find(s[i].v))
        {
            Union(s[i].u,s[i].v);
            insert(s[i].u,s[i].v,s[i].w);insert(s[i].v,s[i].u,s[i].w);
        }
    dfs(1);
    for (int i=1;i<=k;i++)
    {
        int u,v;in(u);in(v);
        int t=lca(u,v);
        printf("%d\n",max(query(u,t),query(v,t)));
    }
}