洛谷 P3128 [USACO15DEC] Max Flow P

maniubi發表於2024-08-28

洛谷 P3128 [USACO15DEC] Max Flow P

題意

給定一棵 \(n\) 個點的樹,給定 \(k\) 個點對 \((u,v)\),把 \(u\)\(v\) 路徑上所有點的點權加一,最後求最大點權。

思路

樹上差分模版。

維護 \(a_i\) 表示每個點到根的加法標記。

對於每個點對 \((u,v)\),把 \(a_u\)\(a_v\) 加一,\(a_{LCA(u,v)}\)\(a_{fa(LCA(u,v))}\) 減一。

因為 \(u\) 到根全部加一,\(v\) 到根全部加一,會讓 \(LCA(u,v)\) 多算一次(因為它本身也要加),\(fa(LCA(u,v))\) 多算一次(\(LCA(u,v)\) 減過一次)。

最後遍歷一遍數計算答案即可。

程式碼

#include <bits/stdc++.h>
#define ll long long
using namespace std;
const int N = 5e4 + 5;
int tot, ver[N << 1], nxt[N << 1], head[N];
int n, k, fa[N][20], de[N];
ll a[N], ans;
void add(int x, int y) {
	ver[++ tot] = y;
	nxt[tot] = head[x];
	head[x] = tot;
}
void DFS(int x) {
    de[x] = de[fa[x][0]] + 1;
    for (int i = 1; i <= 19; i ++)
        fa[x][i] = fa[fa[x][i - 1]][i - 1];
    for (int i = head[x], y; i; i = nxt[i]) {
        y = ver[i];
        if (y == fa[x][0]) continue;
        fa[y][0] = x;
        DFS(y);
    }
}
inline int LCA(int x, int y) {
    if (de[x] < de[y]) swap(x, y);
    for (int i = 19; i >= 0; i --) 
        if (de[fa[x][i]] >= de[y]) x = fa[x][i];
    if (x == y) return x;
    for (int i = 19; i >= 0; i --)
        if (fa[x][i] != fa[y][i]) x = fa[x][i], y = fa[y][i];
    return fa[x][0];
}
void dfs(int x) {
	for (int i = head[x], y; i; i = nxt[i]) {
		y = ver[i];
		if (y == fa[x][0]) continue;
		dfs(y); a[x] += a[y]; // 字首和
	}
}
int main() {
	cin >> n >> k;
	for (int i = 1, u, v; i < n; i ++) {
		cin >> u >> v;
		add(u, v); add(v, u);
	}
	DFS(1);
	for (int i = 1, u, v; i <= k; i ++) {
		cin >> u >> v;
		a[u] ++, a[v] ++, a[fa[LCA(u, v)][0]] --, a[LCA(u, v)] --; // 差分
	}
	dfs(1); // 統計答案
	for (int i = 1; i <= n; i ++) ans = max(ans, a[i]);
	cout << ans << "\n";
	return 0;
}