【LCA】Gym101669I Divide and Conquer

$Source:$
2017-2018 ACM-ICPC Southeastern European Regional Programming Contest (SEERC 2017)
$Problem:$
$n=1e5,m=2n-1$的图，其中$n-1$条边构成A树，另$n-1$条边构成B树。求全局最小割和方案数。
$Idea:$
容易发现最小割的大小不会超过3，否则假设每个点的度数为4，那么至少需要$2n$条边。
那么枚举A树的每一条割边，查找B树中有多少条树边连通着A树的两个连通块，反过来考虑，即B树的每一条树边$(u,v)$，对A树的哪些割边有贡献，显然是A树上$u$到$v$路径上的所有边。B对A做一遍，A再对B做一遍即可。
$Code:$

#include<bits/stdc++.h>
using namespace std;

#define I inline
#define pb push_back
typedef vector<int> VI;

const int N = 1e5+10;

VI G[N];
namespace LCA {
    int dep[N], Log[N], p[N][20];
    I void Log_init() { for(int i = 2; i < N; i++) Log[i] = Log[i>>1]+1; }
    void dfs(int u, int fa) {
        p[u][0] = fa, dep[u] = dep[fa]+1;
        for(int i = 1; i < 20; i++) p[u][i] = p[p[u][i-1]][i-1];
        for(int &v:G[u]) if(v != fa) dfs(v, u);
    }
    I int kfa(int x, int k) { while(k) x = p[x][Log[k]], k -= 1<<Log[k]; return x; }
    I int lca(int x, int y) {
        if(dep[x] < dep[y]) swap(x, y);
        x = kfa(x, dep[x]-dep[y]);
        if(x == y) return x;
        for(int k = Log[dep[x]-1]; k >= 0; k--)
            if(p[x][k] != p[y][k]) x = p[x][k], y = p[y][k];
        return p[x][0];
    }
}using LCA::lca;

int n, dp[N], mn = 4, ans;
VI A[N], B[N];

void dfs(int u, int fa) {
    for(int v:G[u]) if(v != fa) { dfs(v, u); dp[u] += dp[v]; }
    if(u == 1) return;
    if(dp[u]+1 < mn) { mn = dp[u]+1, ans = 1; }
    else if(dp[u]+1 == mn) ans++;
}

void solve(VI *A, VI *B) {
    for(int i = 1; i <= n; i++) G[i].clear(), dp[i] = 0;
    for(int u = 1; u <= n; u++)
    for(int &v:A[u]) G[u].pb(v), G[v].pb(u);
    LCA::dfs(1, 0);
    for(int u = 1; u <= n; u++)
    for(int &v:B[u]) { dp[u]++, dp[v]++; dp[lca(u, v)]-=2; }
    dfs(1, 0);
}

int main() {
    LCA::Log_init();
    scanf("%d", &n);
    for(int i = 1; i < n; i++) {
        int u, v; scanf("%d%d", &u, &v);
        A[u].pb(v);
    }
    for(int i = 1; i < n; i++) {
        int u, v; scanf("%d%d", &u, &v);
        B[u].pb(v);
    }
    solve(A, B); solve(B, A);
    if(mn == 3) printf("3 %d\n", ans);
    else printf("2 %d\n", ans/2);
    return 0;
}