换根DP

换根(DP)属于树形(DP),是二次扫描,因为这类问题通常不会指定根结点,并且根节点的变化会对一些值,例如子结点深度和、点权和等产生影响。
通常需要二次(DFS),第一次(DFS)来处理深度,以及点权和之类的问题,第二次(DFS)开始动态规划。

(u)为当前结点,(v)为当前结点的叶子结点。首先需要用(s_i)来表示以(i)为根的子树中的结点个数,并且有(s_u = 1 + sum s_v),那么这个通过一遍(DFS)就可以办得到。

考虑状态转移,令(f_u)表示以(u)为根时,所有结点深度之和,那么进行换根,让(v)为根,那么转移就是(f_v <- f_u),那么这时候各结点的深度可以分为一下两种情况:

  • 所有在(v)的子树上的结点深度都减少(1),那么总深度之和就减少(s_v)
  • 所有不在(v)的子树上的结点深度都增加(1),那么总深度和就增加(n - s_v)
    那么状态转移方程就可以推出来了,(f_v = f_u - s_v + n - s_v = f_u + n - 2 * s_v)
    于是第二遍(DFS)处理这个东西就好,这样就能求出以所有结点为根节点的深度和了。

例题:[POI2008]STA-Station
代码:

// Problem: P3478 [POI2008]STA-Station
// Contest: Luogu
// URL: https://www.luogu.com.cn/problem/P3478
// Memory Limit: 125 MB
// Time Limit: 2000 ms
// 
// Powered by CP Editor (https://cpeditor.org)

#include <bits/stdc++.h>

using namespace std;

const int N = 1E6 + 10, M = N * 2;

int Size[N];
int h[N], e[M], ne[M], idx;
int depth[N];
long long f[N];
int n;

void add(int a, int b) {
	e[idx] = b, ne[idx] = h[a], h[a] = idx++;
}

int dfs1(int u, int fa) {
	Size[u] = 1;
	if (fa == -1) depth[u] = 0;
	else depth[u] = depth[fa] + 1;
	
	for (int i = h[u]; ~i; i = ne[i]) {
		int j = e[i];
		if (j == fa) continue;
		int S = dfs1(j, u);
		Size[u] += S;
	}
	
	return Size[u];
}

void dfs2(int u, int fa) {
	for (int i = h[u]; i != -1; i = ne[i]) {
		int v = e[i];
		if (v == fa) continue;
		f[v] = f[u] - 2 * Size[v] + n;
		dfs2(v, u);	
	}
}

int main() {
	scanf("%d", &n);
	memset(h, -1, sizeof h);
	for (int i = 1; i <= n - 1; i++) {
		int a, b;
		scanf("%d%d", &a, &b);
		add(a, b), add(b, a);	
	}
	
	dfs1(1, -1);
	
	for (int i = 1; i <= n; i++) f[1] += depth[i];
	
	//换根
	//从u->v,在v子树上的深度减一,一共减去Size[v]
	//不在v子树上的,深度都增加1,总共增加n - Size[v]
	//那么从u到v答案变化,f[v] = f[u] - Size[v] + n - Size[v] = f[u] - 2 * Size[v] + n
	dfs2(1, -1);
	
	int res = 1;
	long long cnt = 0;
	for (int i = 1; i <= n; i++) {
		if (f[i] > cnt) {
			cnt = f[i];
			res = i;
		}
	}
	
	printf("%d
", res);
	
    return 0;
}
原文地址:https://www.cnblogs.com/ZhengLijie/p/15357762.html