Skip to content

2846. Minimum Edge Weight Equilibrium Queries in a Tree 👍

  • Time: $O(n\log n + q(26 + \log n))$
  • Space: $O(n\log n)$
 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
class Solution {
 public:
  vector<int> minOperationsQueries(int n, vector<vector<int>>& edges,
                                   vector<vector<int>>& queries) {
    constexpr int kMax = 26;
    const int m = log2(n) + 1;
    vector<int> ans;
    vector<vector<pair<int, int>>> graph(n);
    // jump[i][j] := the node you reach after jumping 2^j from i
    vector<vector<int>> jump(n, vector<int>(m));
    // count[i][j] := the count of j from root to i, where 1 <= j <= 26
    vector<vector<int>> count(n);
    // depth[i] := the depth of i
    vector<int> depth(n);

    for (const vector<int>& edge : edges) {
      const int u = edge[0];
      const int v = edge[1];
      const int w = edge[2];
      graph[u].emplace_back(v, w);
      graph[v].emplace_back(u, w);
    }

    count[0] = vector<int>(kMax + 1);
    dfs(graph, 0, /*prev=*/-1, 0, jump, count, depth);

    // Calculate binary lifting.
    for (int j = 1; j < m; ++j)
      for (int i = 0; i < n; ++i)
        jump[i][j] = jump[jump[i][j - 1]][j - 1];

    for (const vector<int>& query : queries) {
      const int u = query[0];
      const int v = query[1];
      const int lca = getLCA(u, v, jump, depth);
      // the number of edges between (u, v).
      const int numEdges = depth[u] + depth[v] - 2 * depth[lca];
      // the maximum frequency of edges between (u, v)
      int maxFreq = 0;
      for (int j = 1; j <= kMax; ++j)
        maxFreq = max(maxFreq, count[u][j] + count[v][j] - 2 * count[lca][j]);
      ans.push_back(numEdges - maxFreq);
    }

    return ans;
  }

 private:
  void dfs(const vector<vector<pair<int, int>>>& graph, int u, int prev, int d,
           vector<vector<int>>& jump, vector<vector<int>>& count,
           vector<int>& depth) {
    if (prev != -1)
      jump[u][0] = prev;
    depth[u] = d;
    for (const auto& [v, w] : graph[u]) {
      if (v == prev)
        continue;
      // Inherit the count from the parent.
      count[v] = count[u];
      // Add one to this edge.
      ++count[v][w];
      dfs(graph, v, u, d + 1, jump, count, depth);
    }
  }

  // Returns the lca(u, v) via Calculate binary lifting.
  int getLCA(int u, int v, const vector<vector<int>>& jump,
             const vector<int>& depth) {
    // v is always deeper than u.
    if (depth[u] > depth[v])
      return getLCA(v, u, jump, depth);
    // Jump v to the same height of u.
    for (int j = 0; j < jump[0].size(); ++j)
      if (depth[v] - depth[u] >> j & 1)
        v = jump[v][j];
    if (u == v)
      return u;
    // Jump u and v to the node right below the lca.
    for (int j = jump[0].size() - 1; j >= 0; --j)
      if (jump[u][j] != jump[v][j]) {
        u = jump[u][j];
        v = jump[v][j];
      }
    return jump[v][0];
  }
};
 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
class Solution {
  public int[] minOperationsQueries(int n, int[][] edges, int[][] queries) {
    final int kMax = 26;
    final int m = (int) (Math.log(n) / Math.log(2)) + 1;
    int[] ans = new int[queries.length];
    List<Pair<Integer, Integer>>[] graph = new List[n];
    // jump[i][j] := the node you reach after jumping 2^j from i
    int[][] jump = new int[n][m];
    // count[i][j] := the count of j from root to i, where 1 <= j <= 26
    int[][] count = new int[n][];
    // depth[i] := the depth of i
    int[] depth = new int[n];

    for (int i = 0; i < n; ++i)
      graph[i] = new ArrayList<>();

    for (int[] edge : edges) {
      final int u = edge[0];
      final int v = edge[1];
      final int w = edge[2];
      graph[u].add(new Pair<>(v, w));
      graph[v].add(new Pair<>(u, w));
    }

    count[0] = new int[kMax + 1];
    dfs(graph, 0, /*prev=*/-1, 0, jump, count, depth);

    // Calculate binary lifting.
    for (int j = 1; j < m; ++j)
      for (int i = 0; i < n; ++i)
        jump[i][j] = jump[jump[i][j - 1]][j - 1];

    for (int i = 0; i < queries.length; ++i) {
      final int u = queries[i][0];
      final int v = queries[i][1];
      final int lca = getLCA(u, v, jump, depth);
      // the number of edges between (u, v).
      final int numEdges = depth[u] + depth[v] - 2 * depth[lca];
      // the maximum frequency of edges between (u, v)
      int maxFreq = 0;
      for (int j = 1; j <= kMax; ++j)
        maxFreq = Math.max(maxFreq, count[u][j] + count[v][j] - 2 * count[lca][j]);
      ans[i] = numEdges - maxFreq;
    }

    return ans;
  }

  private void dfs(List<Pair<Integer, Integer>>[] graph, int u, int prev, int d, int[][] jump,
                   int[][] count, int[] depth) {
    if (prev != -1)
      jump[u][0] = prev;
    depth[u] = d;
    for (Pair<Integer, Integer> pair : graph[u]) {
      final int v = pair.getKey();
      final int w = pair.getValue();
      if (v == prev)
        continue;
      // Inherit the count from the parent.
      count[v] = count[u].clone();
      // Add one to this edge.
      ++count[v][w];
      dfs(graph, v, u, d + 1, jump, count, depth);
    }
  }

  // Returns the lca(u, v) via Calculate binary lifting.
  private int getLCA(int u, int v, int[][] jump, int[] depth) {
    // v is always deeper than u.
    if (depth[u] > depth[v])
      return getLCA(v, u, jump, depth);
    // Jump v to the same height of u.
    for (int j = 0; j < jump[0].length; ++j)
      if ((depth[v] - depth[u] >> j & 1) == 1)
        v = jump[v][j];
    if (u == v)
      return u;
    // Jump u and v to the node right below the lca.
    for (int j = jump[0].length - 1; j >= 0; --j)
      if (jump[u][j] != jump[v][j]) {
        u = jump[u][j];
        v = jump[v][j];
      }
    return jump[v][0];
  }
}
 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
class Solution:
  def minOperationsQueries(
      self,
      n: int,
      edges: list[list[int]],
      queries: list[list[int]],
  ) -> list[int]:
    kMax = 26
    m = int(math.log2(n)) + 1
    ans = []
    graph = [[] for _ in range(n)]
    # jump[i][j] := the node you reach after jumping 2^j from i
    jump = [[0] * m for _ in range(n)]
    # count[i][j] := the count of j from root to i, where 1 <= j <= 26
    count = [[] for _ in range(n)]
    # depth[i] := the depth of i
    depth = [0] * n

    for u, v, w in edges:
      graph[u].append((v, w))
      graph[v].append((u, w))

    def dfs(u: int, prev: int, d: int):
      if prev != -1:
        jump[u][0] = prev
      depth[u] = d
      for v, w in graph[u]:
        if v == prev:
          continue
        # Inherit the count from the parent.
        count[v] = count[u][:]
        # Add one to this edge.
        count[v][w] += 1
        dfs(v, u, d + 1)

    count[0] = [0] * (kMax + 1)
    dfs(0, -1, 0)

    # Calculate binary lifting.
    for j in range(1, m):
      for i in range(n):
        jump[i][j] = jump[jump[i][j - 1]][j - 1]

    def getLCA(u: int, v: int) -> int:
      """Returns the lca(u, v) via Calculate binary lifting."""
      # v is always deeper than u.
      if depth[u] > depth[v]:
        return getLCA(v, u)
      # Jump v to the same height of u.
      for j in range(m):
        if depth[v] - depth[u] >> j & 1:
          v = jump[v][j]
      if u == v:
        return u
      # Jump u and v to the node right below the lca.
      for j in range(m - 1, -1, -1):
        if jump[u][j] != jump[v][j]:
          u = jump[u][j]
          v = jump[v][j]
      return jump[v][0]

    for u, v in queries:
      lca = getLCA(u, v)
      # the number of edges between (u, v).
      numEdges = depth[u] + depth[v] - 2 * depth[lca]
      # the maximum frequency of edges between (u, v)
      maxFreq = max(count[u][j] + count[v][j] - 2 * count[lca][j]
                    for j in range(1, kMax + 1))
      ans.append(numEdges - maxFreq)

    return ans