前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >BZOJ3732: Network(Kruskal重构树)

BZOJ3732: Network(Kruskal重构树)

作者头像
attack
发布2018-07-27 15:27:48
3560
发布2018-07-27 15:27:48
举报

题意

Link

给出一张$n$个点的无向图,每次询问两点之间边权最大值最小的路径

$n \leqslant 15000, m \leqslant 30000, k \leqslant 20000$

Sol

很显然答案一定在最小生成树上,但是此题还有一个更为玄学的做法—Kruskal重构树

它是在Kruskal算法上改进而来的。

算法流程:

  1. 对于此题来说,将边权从小到大排序
  2. 用并查集维护两点的联通性,若祖先不相同,那么新建一个节点,权值为边权。左右儿子分别为两个点

这样建出来的树,我们称之为Kruskal重构树

它有许多美妙的性质

  • 是一颗二叉树
  • 两点的LCA的点权为原图中最大值最小的路径上的最大值
  • 任意点的权值大于左右儿子的权值,是一个大根堆

对于此题的样例来说,建出来的图大概长这样

#include<cstdio>
#include<algorithm>
using namespace std;
const int MAXN = 1e5 + 10, INF = 1e9, B = 19;
inline int read() {
    char c = getchar(); int x = 0, f = 1;
    while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
    while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
    return x * f;
}
int fa[MAXN], f[MAXN][21], ch[MAXN][2], cnt, val[MAXN], deep[MAXN];
int N, M, K;
struct Edge {
    int u, v, w;
    bool operator < (const Edge &rhs) const {
        return w < rhs.w;
    }
}E[MAXN];
int head[MAXN], num = 0;
inline void AddEdge(int x, int y, int z) {
    E[++num] = (Edge){x, y, z};
}
int find(int x) {
    if(fa[x] == x) return fa[x];
    else return fa[x] = find(fa[x]);
}
void Kruskal() {
    sort(E + 1, E + num + 1);
    int tot = 0;
    for(int i = 1; i <= M; i++) {
        int x = E[i].u, y = E[i].v;
        int fx = find(x), fy = find(y);
        if(fx == fy) continue;
        ch[++cnt][0] = fx, ch[cnt][1] = fy;
        fa[fa[x]] = fa[fa[y]] = f[fa[x]][0] = f[fa[y]][0] = cnt;
        val[cnt] = E[i].w;
        
    }
}
void dfs(int x) {
    if(!ch[x][0] && !ch[x][1]) return ;
    deep[ch[x][0]] = deep[ch[x][1]] = deep[x] + 1;
    dfs(ch[x][0]); dfs(ch[x][1]);
}
int LCA(int x, int y) {
    if(deep[x] < deep[y]) swap(x, y);
    for(int i = B; i >= 0; i--) 
        if(deep[f[x][i]] >= deep[y])
            x = f[x][i];
    if(x == y) return x;
    for(int i = B; i >= 0; i--)
        if(f[x][i] != f[y][i])
            x = f[x][i], y = f[y][i];
    return f[x][0];
}
main() {
    //freopen("a.in", "r", stdin);
    cnt = N = read(); M = read(); K = read();
    for(int i = 1; i <= N << 1; i++) fa[i] = i;
    for(int i = 1; i <= M; i++) {
        int x = read(), y = read(), z = read();
        AddEdge(x, y, z);
    }
    Kruskal();
    deep[cnt] = 1; dfs(cnt);
    for(int i = 1; i <= B; i++)
        for(int j = 1; j <= 2 * N; j++)
            f[j][i] = f[f[j][i - 1]][i - 1];
    while(K--) {
        int  x = read(), y = read();
    //    printf("%d\", LCA(x, y));
        printf("%d\n", val[LCA(x, y)]);
    }
    return 0;
}
本文参与 腾讯云自媒体分享计划,分享自作者个人站点/博客。
原始发表:2018-07-22 ,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 作者个人站点/博客 前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体分享计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
目录
  • 题意
  • Sol
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档