ツリーの重心分解 (木の重心分解) の図解

# はじめに
大きなサイズの問題を扱っていると、元の問題の**サイズを半分以下**にしたい場面がしばしば登場します。**二分探索**もそうですし、**分割統治法**もそうです。

ツリーに関する問題に対しても、サイズを半分以下にしたい場面がしばしば登場します。そんなときに重要な概念になるのが**ツリーの重心**です。


# ツリーの重心とは
ツリーである頂点を取り除くと、元のツリーは何個かのツリーに分裂します。そのすべての部分ツリーのサイズが元のツリーのサイズの**半分以下**になるような頂点がツリーの重心です。

![ツリーの重心.jpg](https://qiita-image-store.s3.amazonaws.com/0/182963/5862e3e8-0f2f-02ee-dc3e-c75dcd884061.jpeg)


　  

なおこの図の場合、もう一つ重心があることがわかります。 

![ツリーの重心2個.jpg](https://qiita-image-store.s3.amazonaws.com/0/182963/fad64e0f-4083-3bc5-7e3d-9cd3ec647a25.jpeg)

　

一般に、**ツリーの重心は 1 個または 2 個**になります。多くのツリーでは重心は 1 個だけですが、この例については、2 つの重心を結んでいる枝の両端からなる部分木がともに元のサイズの半分ずつになっています。このような状況が発生する場合のみ、重心が 2 個になります。なお、重心が 2 個以下であることの詳しい証明は、以下のこうきさんの記事

- [木の重心列挙アルゴリズム](http://www.learning-algorithms.com/entry/2018/01/03/215559)

に記載があります。


<img src="https://qiita-image-store.s3.amazonaws.com/0/182963/7d630c40-a057-2981-cb40-ab7299a86dd8.jpeg" width=500>



# ツリーの重心が必ず存在すること
ツリーに対してサイズを半分以下にしながら問題を解きたい場面では、**ツリーの重心が存在するかどうか**はとても重要な問題です。ツリーの重心が必ず存在することは、パッと見ではそんなに自明ではないです。

しかしながら次の章に示すように、ツリーの重心の求め方を見れば、重心の存在が明らかになります。


# ツリーの重心を求める (ツリーの重心の存在証明)
ツリー DP の心得があれば簡単に求めることができます (N: ツリーの頂点数)。

-----
1. 根を適当に決めて根付き木にする
2. 各ノード v を根とした部分木のサイズ sizeSubtree[v] をツリー DP によって求めていく
3. 任意の v の子ノード ch に対して sizeSubtree[ch] <= N/2 が成り立ち、かつ N - sizeSubtree[v] <= N/2 を満たすような v を列挙する

-----

これだけです。以下が実装例になります。

```cpp
const int MAX_V = 2100;  // ツリーのサイズのありうる最大値

int N;  // ツリーのサイズ
vector<int> tree[MAX_V];  // ツリーを隣接リスト形式のグラフ構造で表したもの

int sizeSubtree[MAX_V];  // sizeSubtree[v] := v を根とする部分ツリーのサイズ
vector<int> centroids;  // 重心列挙の答えがここに入る

/* ツリーDP */
void SubFindCentroids(int v, int parent_of_v = -1) {
	sizeSubtree[v] = 1;
	bool isCentroid = true;
	for (auto ch : tree[v]) {
		if (ch == parent_of_v) continue;
		SubFindCentroids(ch, v);
		if (sizeSubtree[ch] > N / 2) isCentroid = false;
		sizeSubtree[v] += sizeSubtree[ch];
	}
	if (N - sizeSubtree[v] > N / 2) isCentroid = false;
	if (isCentroid) centroids.push_back(v);
}

void FindCentroids() {
	centroids.clear();
	SubFindCentroids(0, N); // これを呼び出すと centoroids に重心を列挙したものが入る
}
```


なお、[プログラミングコンテストチャレンジブック](https://www.amazon.co.jp/%E3%83%97%E3%83%AD%E3%82%B0%E3%83%A9%E3%83%9F%E3%83%B3%E3%82%B0%E3%82%B3%E3%83%B3%E3%83%86%E3%82%B9%E3%83%88%E3%83%81%E3%83%A3%E3%83%AC%E3%83%B3%E3%82%B8%E3%83%96%E3%83%83%E3%82%AF-%E7%AC%AC2%E7%89%88-%EF%BD%9E%E5%95%8F%E9%A1%8C%E8%A7%A3%E6%B1%BA%E3%81%AE%E3%82%A2%E3%83%AB%E3%82%B4%E3%83%AA%E3%82%BA%E3%83%A0%E6%B4%BB%E7%94%A8%E5%8A%9B%E3%81%A8%E3%82%B3%E3%83%BC%E3%83%87%E3%82%A3%E3%83%B3%E3%82%B0%E3%83%86%E3%82%AF%E3%83%8B%E3%83%83%E3%82%AF%E3%82%92%E9%8D%9B%E3%81%88%E3%82%8B%EF%BD%9E-%E7%A7%8B%E8%91%89%E6%8B%93%E5%93%89/dp/4839941068)の P.320 には、より**重心の存在証明**を意識した構成方法が以下のように示されています:

-----
1. 根を適当に決めて根付き木とし、v ← 根とする
2. While True:
3. 　　If v の子ノードたちからなる部分木たちのサイズがすべて N/2 以下: 
4. 　　　　v を重心として break
5. 　　v の子ノードのうち、最も部分木のサイズが大きいものを ch として、v ← ch

-----

つまり root から出発して、「子ノードを根とする部分木のうち一番大きいやつに侵入」を繰り返して、初めて「子供が全員 N/2 以下」になる部分で break します。このアルゴリズムが停止することは明らかです。root から出発して、最悪 leaf まで行けば部分木たちのサイズの最大値が 0 になるので、それまでにどこかで停止します。

そしてよくよく考えるとこのアルゴリズムが停止したときの v は重心になっています！
v の子供たちはみんなサイズ N/2 以下ですが、図のように v の親のサイズも N/2 以下になっていることがわかります。以上から、重心が常に存在することもわかりました。
　


![ツリーの重心の存在.jpg](https://qiita-image-store.s3.amazonaws.com/0/182963/9d1e5186-cda3-a189-f8b5-dd834ddfbdd9.jpeg)

　





# ツリー上の分割統治法: ツリーの重心分解を再帰的に行う方法
ツリーの重心分解は単発で行うケースもありますが、多くの場合は**再帰的に**重心分解していきます。
再帰的に重心分解していく上で難しいのは

- 再帰的に作られていく部分ツリーをどのように表現するか
- 重心分解によって作られた部分ツリーにさらに重心分解を適用するときに、その部分ツリー全体の頂点数を予め知ること

1 つ目については、既に取り除かれた頂点をフラグ管理しておけば実装することができます。
2 つ目については、ここでは重心を求めるために使用したツリーDPの結果をフル活用する方針で書いてみます。

再帰的な重心分解を意識した実装例を以下に示します。
より実践的な分割統治の書き方については、[こうきさんの記事](http://www.learning-algorithms.com/entry/2018/01/20/031005)にまとまっています。

```cpp
const int MAX_V = 2100;  // ツリーのサイズのありうる最大値

int N;  // ツリーのサイズ
vector<int> tree[MAX_V];  // ツリーを隣接リスト形式のグラフ構造で表したもの

int sizeSubtree[MAX_V];  // sizeSubtree[v] := v を根とする部分ツリーのサイズ (分割統治の毎ステップごとに再利用)
bool isRemoved[MAX_V];  // isRemoved[v] := v が既に取り除かれたかどうか
int whoIsParent[MAX_V];  // whoIsParent[v] := ツリーDP時に v の親が誰だったか

// メイン処理
vector<int> centroids;
void FindCentroidRecursive(int v, int size, int p = -1) {
	sizeSubtree[v] = 1;
	whoIsParent[v] = p;
	bool isCentroid = true;
	for (auto ch : tree[v]) {
		if (ch == p) continue;
		if (isRemoved[ch]) continue;
		FindCentroidRecursive(ch, size, v);
		if (sizeSubtree[ch] > size / 2) isCentroid = false;
		sizeSubtree[v] += sizeSubtree[ch];
	}
	if (size - sizeSubtree[v] > size / 2) isCentroid = false;
	if (isCentroid) centroids.push_back(v);
}

// 初期化
void Init() {
	for (int i = 0; i < MAX_V; ++i) {
		isRemoved[i] = false;
	}
}

// first: 重心, second: (重心の子ノード, 子部分木のサイズ) からなるベクトル
pair<int, vector<pair<int,int> > > FindCentroid(int root, int size) {
	vector<pair<int, int> > subtrees;
	centroids.clear();
	FindCentroidRecursive(root, size);
	int center = centroids[0];
	for (auto ch : tree[center]) {
		if (isRemoved[ch]) continue;
		if (ch == whoIsParent[center]) {
			subtrees.push_back(make_pair(ch, size - sizeSubtree[center]));
		}
		else {
			subtrees.push_back(make_pair(ch, sizeSubtree[ch]));
		}
	}
	return make_pair(center, subtrees);
}
```


# ツリーの重心分解を用いる具体的な問題たち
とてもよくまとまっている記事たちがあるのでそれらへのポインタを示します:

- [競技プログラミングにおける重心分解問題まとめ (はまやんさん)](http://hamayanhamayan.hatenablog.jp/entry/2017/12/19/152143)
- [木の重心列挙アルゴリズム (こうきさん)](http://www.learning-algorithms.com/entry/2018/01/03/215559)
- [重心分解による分割統治法の一般形について (こうきさん)](http://www.learning-algorithms.com/entry/2018/01/20/031005)


# おわりに
本記事ではツリーの重心という概念の図解を試みました。