部分树上问题及图的联通性（图论学习总结部分内容）

smile是对你的礼貌~@济南大学2024-05-15 8:05

文章目录

    • 前言
    • 三、部分树上问题及图的联通性
      • 最小生成树
        • 知识点
        • 例题
          • [e g 1 : eg1: eg1: 走廊泼水节（克鲁斯卡尔思想的灵活运用）](#e g 1 : eg1: eg1: 走廊泼水节（克鲁斯卡尔思想的灵活运用）)
          • [e g 2 ： eg2： eg2： B-Picnic Planning](#e g 2 ： eg2： eg2： B-Picnic Planning)
          • [e g 3 eg3 eg3：L - Classic Problem（缩点+ B o r u v k a Boruvka Boruvka 算法， 2023 C C P C 2023CCPC 2023CCPC广东防AK题）](#e g 3 eg3 eg3：L - Classic Problem（缩点+ B o r u v k a Boruvka Boruvka 算法， 2023 C C P C 2023CCPC 2023CCPC广东防AK题）)
      • 生成树计数
      • LCA、树上差分应用与树上问题
        • 知识点
        • 例题
          • [e g 1 : eg1: eg1: E. Boomerang（2024ICPC全国邀请赛(武汉)，树的直径，银牌题）](#e g 1 : eg1: eg1: E. Boomerang（2024ICPC全国邀请赛(武汉)，树的直径，银牌题）)
          • [e g 2 : eg2: eg2: 闇の連鎖 (树上差分)](#e g 2 : eg2: eg2: 闇の連鎖 (树上差分))
          • [e g 3 : eg3: eg3: 天天爱跑步(蓝书例题)](#e g 3 : eg3: eg3: 天天爱跑步(蓝书例题))
          • [e g 4 : eg4: eg4:​ 异象石(蓝书例题)](#e g 4 : eg4: eg4: 异象石(蓝书例题))
          • [e g 5 : eg5: eg5:​ 次小生成树(蓝书例题)](#e g 5 : eg5: eg5: 次小生成树(蓝书例题))
          • [e g 6 : eg6: eg6: 疫情控制(蓝书例题)](#e g 6 : eg6: eg6: 疫情控制(蓝书例题))
      • 基环树
          • [e g 1 : eg1: eg1: Island （洛谷P4381 [IOI2008]）](#e g 1 : eg1: eg1: Island （洛谷P4381 [IOI2008]）)
          • [e g 2 : eg2: eg2: 创世纪(蓝书例题)](#e g 2 : eg2: eg2: 创世纪(蓝书例题))
          • [e g 3 : eg3: eg3: 牛镇公务员考试（2024牛客寒假营，基环树一般的处理方式）](#e g 3 : eg3: eg3: 牛镇公务员考试（2024牛客寒假营，基环树一般的处理方式）)
          • [e g 4 : S u g a r S w e e t I I eg4:Sugar\ Sweet II eg4:Sugar SweetII​ （基环树 + 概率，2023ICPC杭州站1/2个银牌题）](#e g 4 : S u g a r S w e e t I I eg4:Sugar\ Sweet II eg4:Sugar SweetII （基环树 + 概率，2023ICPC杭州站1/2个银牌题）)
      • Tarjan与无向图的连通性
        • 例题
          • [e g 1 : eg1: eg1: 嗅探器](#e g 1 : eg1: eg1: 嗅探器)
          • [e g 3 eg3 eg3 ： [ H N O I 2012 ] [HNOI2012] [HNOI2012]​​ 矿场搭建](#e g 3 eg3 eg3 ： [ H N O I 2012 ] [HNOI2012] [HNOI2012] 矿场搭建)
      • Tarjan与有向图的连通性
        • 例题
          • [e g 1 : eg1: eg1: [ZJOI2007] 最大半连通子图](#e g 1 : eg1: eg1: [ZJOI2007] 最大半连通子图)
      • 练习题

前言

由于图论学习总结内容过多，全放在一篇博客过于冗长现进行拆分，本文是部分树上问题及图的联通性 部分，其他部分地址见：图论学习总结(For XCPC)

三、部分树上问题及图的联通性

最小生成树

知识点
  • K r u s k a l Kruskal Kruskal 算法
  • P r i m Prim Prim 算法
  • B o r u v k a Boruvka Boruvka 算法
  • 次小生成树
  • 瓶颈生成树
  • 最小瓶颈路
  • K r u s k a l Kruskal Kruskal 重构树
例题
e g 1 : eg1: eg1: 走廊泼水节（克鲁斯卡尔思想的灵活运用）

A-走廊泼水节(nowcoder.com)

题目大意

给定一棵N个节点的树，要求增加若干条边，把这棵树扩充为完全图，并满足图的唯一最小生成树仍然是这棵树。 求增加的边的权值总和最小是多少。

e g 2 ： eg2： eg2： B-Picnic Planning

B-Picnic Planning (nowcoder.com)（码量较大）

题目大意

e g 3 eg3 eg3：L - Classic Problem（缩点+ B o r u v k a Boruvka Boruvka 算法， 2023 C C P C 2023CCPC 2023CCPC广东防AK题）

Problem - L - Classic Problem

题目大意

Given an undirected complete graph with n n n vertices and m m m triples P 1 , P 2 , ⋯   , P m P_1, P_2, \cdots, P_m P1,P2,⋯,Pm where P i = ( u i , v i , w i ) P_i = (u_i, v_i, w_i) Pi=(ui,vi,wi), it's guaranteed that 1 ≤ u i < v i ≤ n 1 \leq u_i < v_i \leq n 1≤ui<vi≤n, and for any two triples P i P_i Pi and P j P_j Pj with different indices we have ( u i , v i ) ≠ ( u j , v j ) (u_i, v_i) \ne (u_j, v_j) (ui,vi)=(uj,vj).

For any two vertices x x x and y y y in the graph ( 1 ≤ x < y ≤ n 1 \leq x < y \leq n 1≤x<y≤n), define the weight of the edge connecting them as follows:

  • If there exists a triple P i P_i Pi satisfying u i = x u_i = x ui=x and v i = y v_i = y vi=y, the weight of edge will be w i w_i wi.
  • Otherwise, the weight of edge will be ∣ x − y ∣ |x - y| ∣x−y∣.

Calculate the total weight of edges in the minimum spanning tree of the graph.

注：这题还未ac 先占坑

生成树计数

LCA、树上差分应用与树上问题

知识点

  • 树的直径

    • 树形DP
    • 两次 B F S BFS BFS （注意运用条件）

  • L C A LCA LCA

    • 树上倍增法
    • T a r j a n Tarjan Tarjan 算法 (离线算法)

  • 树链剖分

  • 树上hash

    • 在判断一些树是否同构时，常常把这些树转成哈希存储起来，降低复杂度
    • 这里给出 O I W i k i OI\ Wiki OI Wiki 的推荐实现:

  • 树分治

例题
e g 1 : eg1: eg1: E. Boomerang（2024ICPC全国邀请赛(武汉)，树的直径，银牌题）

E. Boomerang
题目大意

给定一个 n 个节点的树，根节点为 r, 设 V ( r , t ) = v ∣ d i s ( r , v ) ⩽ t V(r, t) = {v|dis(r, v) ⩽ t} V(r,t)=v∣dis(r,v)⩽t，其中 d i s ( u , v ) dis(u, v) dis(u,v) 表示树上两点 u 和 v 间的唯一简单路径的边数。再给定一个 t 0 t_0 t0，对每个 k k k，你需要选择一个点 r ′ r' r′，定义 V ′ ( r ′ , t ) = v ∣ d i s ( r ′ , v ) ⩽ k ( t − t 0 ) V'(r', t) = {v|dis(r', v) ⩽ k(t−t0)} V′(r′,t)=v∣dis(r′,v)⩽k(t−t0)。找到一个最小的 t t t，使得 V ( r , t ) ⊆ V ′ ( r ′ , t ) V(r, t) ⊆ V'(r', t) V(r,t)⊆V′(r′,t)。
题目分析

  • 对题目进行分析，我们很容易想到一点，就是随着 k k k 的增大， t t t 非严格递减。
  • 我们考虑在 k = i ∈ [ 1 , n ] k=i \in [1,n] k=i∈[1,n] 时 r ′ r' r′ 怎么选择会最优。
  • 假设 t t t 是一个可能的答案，那么显然，子树 V ( r , t ) V(r, t) V(r,t) 的直径的中点就是 r ′ r' r′ 的最佳选择。
  • 到这里，这道题的解法已经很明显了，我们可以维护出每个时刻 t t t，子树 V ( r , t ) V(r, t) V(r,t) 的直径是多少。
  • 从小到达枚举 t t t 判断在 k = i k=i k=i 的情况下是否满足条件。(这里在具体实现时用到了第一条提到的性质)
    参考了jls思路，ac代码(Python，FastIO已省略，其中手写栈 模板已给出)：
python 复制代码 
# --------------------
# 手写栈模板
# 克服py栈太浅的问题
def bootstrap(f, stack=[]):
    def wrappedfunc(*args, **kwargs):
        if stack:
            return f(*args, **kwargs)
        else:
            to = f(*args, **kwargs)
            while True:
                if type(to) is GeneratorType:
                    stack.append(to)
                    to = next(to)
                else:
                    stack.pop()
                    if not stack:
                        break
                    to = stack[-1].send(to)
            return to

    return wrappedfunc


def solve():
    n = I()
    g = [[] for _ in range(n+1)]
    for _ in range(n-1):
        x,y = MI()
        g[x].append(y)
        g[y].append(x)
    r,t0 = MI()

    par, dep, h, dia = [0]*(n+1), [0]*(n+1), [0]*(n+1), [0]*n
    final = 0
    @bootstrap
    def dfs(x) -> None:
        nonlocal final
        v = [0,0]
        for y in g[x]:
            if y == par[x]: continue
            par[y] = x
            dep[y] = dep[x] + 1  # 深度
            yield dfs(y)
            d = h[y] + 1  # 高度
            for i in range(2):
                if d > v[i]: d,v[i] = v[i],d

        h[x] = v[0]
        dia[dep[x] + v[1]] = max(dia[dep[x] + v[1]], 2*v[1])
        final = max(v[0]+v[1], final)
        yield None
    dfs(r)

    for i in range(1,n):
        dia[i] = max(dia[i], min(dia[i-1]+1, final))
    for i in range(n-2,-1,-1):
        dia[i] = max(dia[i], dia[i+1]-2)

    t = t0
    ans = [0]*(n+1)
    for k in range(n,0,-1):
        while (t - t0) * k * 2 < dia[min(t, n-1)]:
            t += 1
        ans[k] = t

    print_arr(ans[1:])


solve()
e g 2 : eg2: eg2: 闇の連鎖 (树上差分)

闇の連鎖

题目大意

​ 根据题意，"主要边"构成一棵树，"附加边"是非树边。把一条附加边 ( x , y ) (x,y) (x,y) 添加到树上会构成一个环，如过一开始切断 x x x~ y y y 路径上的一条边，那么第二条必须切断 ( x , y ) (x,y) (x,y) 才能将 D a r k Dark Dark 分为不连通的两部分。我们可以考虑 ( x , y ) (x,y) (x,y) 把 x x x~ y y y 路径上的树边覆盖了一次，那么如果我们第一次切掉覆盖数为 0 0 0 的树边， 那么之后切任意一条非树边即可，如果我们切掉覆盖数为 1 1 1 的树边，第二步的切分唯一，其余情况无法将 D a r k Dark Dark 分为不连通的两部分。

​ 由上述分析可知，我们要解决问题的模型就是在给定的一张无向图和一课生成树，求非树边将树边覆盖了多少次 。解决这一问题的经典做法就是"树上差分 "。给定每个结点的权值初值为 0 0 0，对于每条非树边 ( x , y ) (x,y) (x,y) 我们让 v a l x + = 1 , v a l y + = 1 , v a l l c a ( x , y ) − = 2 val_x+=1,\ val_y+=1,\ val_{lca(x,y)}-=2 valx+=1, valy+=1, vallca(x,y)−=2，最后用树形DP做一次 "子树前缀和" 即可得到每条树边的覆盖次数，然后统计即可得到答案。

​ 关于时间复杂度：考虑用树上倍增法求 L C A LCA LCA，复杂度为 O ( n log ⁡ n ) O(n\log n) O(nlogn)。对于覆盖次数的预处理和答案求解需要遍历整张图，复杂度为 O ( n + m ) O(n+m) O(n+m)。总的时间复杂度为 O ( n + m + n log ⁡ n ) O(n+m+n\log n) O(n+m+nlogn)。

​ 启发： 如果我们需要对树上的一条路径的每条边加上一个 d d d，我们可以转化为 v a l x + = d , v a l y + = d , v a l l c a ( x , y ) − = 2 d val_x+=d,\ val_y+=d,\ val_{lca(x,y)}-=2d valx+=d, valy+=d, vallca(x,y)−=2d。

e g 3 : eg3: eg3: 天天爱跑步(蓝书例题)

E-天天爱跑步

e g 4 : eg4: eg4:​ 异象石(蓝书例题)

F-异象石_

e g 5 : eg5: eg5:​ 次小生成树(蓝书例题)

G-次小生成树_

e g 6 : eg6: eg6: 疫情控制(蓝书例题)

H-疫情控制_

基环树

e g 1 : eg1: eg1: Island （洛谷P4381 [IOI2008]）

P4381 [IOI2008] Island - 洛谷

题目大意

​ n n n 个点 n n n 条边，可以知道这题给你的是一个基环树森林，而根据题意中渡船的条件可知，当你从一棵基环树到另一棵的时候，你不能再回去了。所以我们可以对每个基环树单独计算。我们不能重复到达同一岛屿，所以我们要求的就是基环树上的一条最长简单路径的长度，那么所有基环树的和就是我们的答案。

​ 现在我们考虑每一棵基环树答案的情况1、不经过环（该情况可树形DP求解，即求树的直径）2、经过环。对于情况二我们要找到环上的两个结点 s i , s j s_i,s_j si,sj，使得 D [ s i ] + D [ s j ] + d i s t [ s i , s j ] D[s_i]+D[s_j]+dist[s_i,s_j] D[si]+D[sj]+dist[si,sj] 最大，其中 d i s t [ s i , s j ] dist[s_i,s_j] dist[si,sj] 是环上两点的距离，我们把环拆开复制一倍，再做线性 D P DP DP 即可求得答案，用单调队列优化，该DP时间复杂度为 O ( P ) O(P) O(P) P P P 为环的点数。总的时间复杂度为 O ( n ) O(n) O(n)。

​ 关于实现，这道题需要 d f s dfs dfs 或 b f s + 拓扑排序 bfs+拓扑排序 bfs+拓扑排序 找连通块和环，还要树形 D P DP DP 求直径，再对环做 D P DP DP 等等，常数比较大， 1 0 6 10^6 106 点 d f s dfs dfs 还可能爆栈c++实现如下：记录详情 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)（ac代码）

python实现如下：

记录详情 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

结果：

e g 2 : eg2: eg2: 创世纪(蓝书例题)

B-创世纪

题目大意

e g 3 : eg3: eg3: 牛镇公务员考试（2024牛客寒假营，基环树一般的处理方式）

K-牛镇公务员考试

题目大意

e g 4 : S u g a r S w e e t I I eg4:Sugar\ Sweet II eg4:Sugar SweetII​ （基环树 + 概率，2023ICPC杭州站1/2个银牌题）

2023_ICPC_杭州站 - H - Codeforces

题目大意

Sugar is sweet. There are n n n children asking for sugar. Prof. Chen gives out sugar to the children. The i i i-th child initially has a i a_{i} ai bags of sugar. There are n n n events happening in uniformly randomized order . The i i i-th event is:

  • If the i i i-th child has strictly less bags of sugar than the b i b_{i} bi-th child, then the i i i-th child will get extra w i w_{i} wi bags of sugar. Otherwise, nothing happens.

Now, since the events happen in random order, Randias, which is the assistant of Prof. Chen, wants to know the expected number of bags of sugar each child will have after all the events happen.

It can be shown that the answer can be expressed as an irreducible fraction x y \frac{x}{y} yx where x x x and y y y are integers and y ≢ 0 ( m o d 1 0 9 + 7 ) y \not \equiv 0 \pmod {10^9 + 7} y≡0(mod109+7). Output the integer equal to x ⋅ y − 1 ( m o d 1 0 9 + 7 ) x \cdot y^{-1} \pmod {10^9 + 7} x⋅y−1(mod109+7). In other words, output such an integer a a a that 0 ≤ a < 1 0 9 + 7 0\leq a < 10^9 + 7 0≤a<109+7 and a ⋅ y ≡ x ( m o d 1 0 9 + 7 ) a \cdot y \equiv x \pmod {10^9 + 7} a⋅y≡x(mod109+7)​.

题目解析

  • 我们分析一下题目，第 i i i 个事件的发生与否只与 a b i a_{b_i} abi 有关，事件的依赖关系构成基环树森林。
  • 我们分析一下依赖 关系，如果 a i < a b i a_i < a_{b_i} ai<abi 那么事件 i i i 必定发生，如果 a i ≥ a b i a_i \ge a_{b_i} ai≥abi 且 a i < a b i + w i a_i < a_{b_i} + w_i ai<abi+wi 那么事件 i i i 必定在事件 b i b_i bi 之后发生，否则一定不发生。
  • 我们要求的是每个小朋友最后的期望糖果数，其实只要求出每个人加糖的期望即可。
  • 对于必定会发生的事件我们让他的概率为 1 1 1，对于可能发生的事件我们对其先后关系进行连边。
  • 假设从一个必定发生的事件 x x x 走向 可能发生的事件 y y y 是基环树上的一条简单路径，那么 y y y 发生的概率就是这条路径长度的阶乘分之一： 1 L i ! \frac{1}{L_i!} Li!1，因为只有事件发生的顺序刚好满足这个排列才行。

ac代码参考(FastIO已省略)

python 复制代码 
mod = 10**9 + 7
N = 5*(10**5)+1
inv = lambda x:pow(x,mod-2,mod)
jc = [1]*N  # 阶乘
jc[0] = 0
for i in range(2,N):
    jc[i] = i*jc[i-1] % mod
 
def solve():
    n = I()
    a = [0]+LI()
    b = [0]+LI()
    w = [0]+LI()
    g = [[] for _ in range(n+1)]
    li = [0]*(n+1)
    deg = [0]*(n+1)
    for i in range(1,n+1):
        if a[i] >= a[b[i]] and a[i] < a[b[i]] + w[b[i]]:
            g[b[i]].append(i)
            deg[i] += 1
        if a[i] < a[b[i]]:
            li[i] = 1
    q = deque()
    for i in range(1,n+1):
        if deg[i] == 0: q.append(i)
    while q:
        x = q.popleft()
        if li[x] == 0: continue
        for y in g[x]:
            li[y] += li[x]+1
            deg[y] -= 1
            if deg[y] == 0: q.append(y)
 
    for i in range(1,n+1):
        a[i] = (a[i] + w[i]*inv(jc[li[i]])) % mod
    print_arr(a[1:])
 
for _ in range(int(input())):
    solve()

Tarjan与无向图的连通性

注：python 的Tarjan算法与连通性相关代码模板和知识点可以看我的这篇博客

例题
e g 1 : eg1: eg1: 嗅探器

嗅探器

题目大意（来源：牛客网）

某军搞信息对抗实战演习，红军成功地侵入了蓝军的内部网络，蓝军共有两个信息中心，红军计划在某台中间服务器上安装一个嗅探器，从而能够侦听到两个信息中心互相交换的所有信息，但是蓝军的网络相当的庞大，数据包从一个信息中心传到另一个信息中心可以不止有一条通路。现在需要你尽快地解决这个问题，应该把嗅探器安装在哪个中间服务器上才能保证所有的数据包都能被捕获？

input

输入文件的第一行一个整数 n ( n ≤ 100 ) n\ (n\le 100) n (n≤100)，表示蓝军网络中服务器的数目。

接下来若干行是对蓝军网络的拓扑结构描述，每行是两个整数 i , j i,j i,j，表示编号为i和编号为j的两台服务器间存在连接（显然连接是双向的），服务器的编号从 1 1 1 开始，一行两个 0 0 0 表示网络的拓扑结构描述结束，再接下来是两个整数 a , b a,b a,b，分别表示两个中心服务器的编号。

output

输出编号。如果有多个解输出编号最小的一个，如果找不到任何解，输出No solution。

题目解析

读完题目很容易想到这道题就是让我们求从 s s s 到 t t t 的必经点或判断无解。首先，一个点如果是必经点，那么其必定是一个割点，如果这个割点将 s s s 和 t t t​ 恰好分隔在两个双连通分量中，那么这就是我们答案的候选之一。

具体实现细节方面，我们可以从 s s s 出发对 T a r j a n Tarjan Tarjan 求割点的算法进行的时候做一个判断，看看搜索树中，如果当前点 x x x 是割点，且 x ! = s , x ! = t , d f n [ y ] ≤ d f n [ t ] x != s,\ x!=t,\ dfn[y]\le dfn[t] x!=s, x!=t, dfn[y]≤dfn[t] 。如果满足 d f n [ y ] ≤ d f n [ t ] dfn[y]\le dfn[t] dfn[y]≤dfn[t] 那么意味着从当前割点出发沿着 y y y 这一方向搜索树的子树中有 t t t，那么此时割点 x x x 必定是必经点。

e g 3 eg3 eg3 ： [ H N O I 2012 ] [HNOI2012] [HNOI2012]​​ 矿场搭建

[ H N O I 2012 ] [HNOI2012] [HNOI2012]​​ 矿场搭建

Tarjan与有向图的连通性

注：python 的Tarjan算法与连通性相关代码模板和知识点可以看我的这篇博客

例题
e g 1 : eg1: eg1: [ZJOI2007] 最大半连通子图

ZJOI2007 最大半连通子图

练习题

B-树网的核_0x63 图论-树的直径与最近公共祖先

C-最优比率生成树_0x62 图论-最小生成树(0/1分数规划)

D-雨天的尾巴 (树上差分的综合运用)

相关推荐
nuclear201126 分钟前
使用Python 在Excel中创建和取消数据分组 - 详解
python·excel数据分组·创建excel分组·excel分类汇总·excel嵌套分组·excel大纲级别·取消excel分组
Lucky小小吴41 分钟前
有关django、python版本、sqlite3版本冲突问题
python·django·sqlite
带多刺的玫瑰44 分钟前
Leecode刷题C语言之统计不是特殊数字的数字数量
java·c语言·算法
爱敲代码的憨仔1 小时前
《线性代数的本质》
线性代数·算法·决策树
yigan_Eins1 小时前
【数论】莫比乌斯函数及其反演
c++·经验分享·算法
GIS 数据栈1 小时前
每日一书 《基于ArcGIS的Python编程秘笈》
开发语言·python·arcgis
爱分享的码瑞哥1 小时前
Python爬虫中的IP封禁问题及其解决方案
爬虫·python·tcp/ip
阿史大杯茶1 小时前
AtCoder Beginner Contest 381（ABCDEF 题）视频讲解
数据结构·c++·算法
დ旧言~2 小时前
【高阶数据结构】图论
算法·深度优先·广度优先·宽度优先·推荐算法
张彦峰ZYF2 小时前
投资策略规划最优决策分析
分布式·算法·金融
热门推荐
01(欧拉)openEuler系统添加网卡文件配置流程、(欧拉)openEuler系统手动配置ipv6地址流程、(欧拉)openEuler系统网络管理说明02PyTorch机器学习实现液态神经网络03【HarmonyOS】HUAWEI DevEco Studio 下载地址汇总04玄机平台应急响应—webshell查杀05Coze扣子平台完整体验和实践（附国内和国际版对比）06Ubuntu 20.04使用Livox mid 360 测试 FAST_LIO07【目标跟踪】相机运动补偿08Unity中PICO实现 隔空取物 和 接触抓取物体09RAG 实践- Ollama+RagFlow 部署本地知识库10怎样让音频速度变慢？请跟随以下方法进行操作