树的重心

这是找树的重心的经典题目。

树的重心有下面几条常见性质：

定义1：找到一个点，其所有的子树中最大的子树节点数最少，那么这个点就是这棵树的重心。
定义2：以这个点为根，那么所有的子树（不算整个树自身）的大小都不超过整个树大小的一半。
性质1：树中所有点到某个点的距离和中，到重心的距离和是最小的；如果有两个重心，那么他们的距离和一样。
性质2：把两个树通过一条边相连得到一个新的树，那么新的树的重心在连接原来两个树的重心的路径上。
性质3：把一个树添加或删除一个叶子，那么它的重心最多只移动一条边的距离。

具体树的重心求法请看：算法竞赛入门经典281页（紫书）

题目:http://poj.org/problem?id=1655

题意:给定一棵树,求树的重心的编号以及重心删除后得到的最大子树的节点个数size,如果size相同就选取编号最小的.

分析:首先要知道什么是树的重心,树的重心定义为:找到一个点,其所有的子树中最大的子树节点数最少,那么这个点就是这棵树的重心,删去重

心后，生成的多棵树尽可能平衡.  实际上树的重心在树的点分治中有重要的作用, 可以避免N^2的极端复杂度（从退化链的一端出发），保证

NlogN的复杂度, 利用树型dp可以很好地求树的重心.

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
using namespace std;
const int MAXN = 20005;
const int INF = 0x3f3f3f3f;
int head[MAXN], tol;
int son[MAXN], ans, siz, n;//son[i]表示以i为根的子树中节点总数
bool vis[MAXN];
struct Edge
{
    int to, nxt;
}edge[MAXN*2];
void init()
{
    tol = 0;
    siz = INF;
    memset(vis, false, sizeof(vis));
    memset(head, -1, sizeof(head));
}
void addedge(int u, int v)
{
    edge[tol].to = v;
    edge[tol].nxt = head[u];
    head[u] = tol++;
}
void dfs(int cur)
{
    vis[cur] = true;
    son[cur] = 0;
    int tmp = 0;
    for(int i = head[cur]; i != -1; i = edge[i].nxt)
    {
        int v = edge[i].to;
        if(!vis[v])
        {
            dfs(v);
            son[cur] += son[v] + 1;
            tmp = max(tmp, son[v]+1);
        }
    }
    tmp = max(tmp, n-son[cur]-1);
    if(tmp < siz || tmp == siz && cur < ans)
    {
        ans = cur;
        siz = tmp;
    }
}
int main()
{
    int T;
    scanf("%d", &T);
    while(T--)
    {
        init();
        scanf("%d", &n);
        for(int i = 1; i <= n-1; i++)
        {
            int u, v;
            scanf("%d%d", &u, &v);
            addedge(u, v);
            addedge(v, u);
        }
        dfs(1);
        printf("%d %d\n", ans, siz);
    }
    return 0;
}

题目:http://poj.org/problem?id=3107

题意:给定一棵树,求树的所有重心,按照编号从小到大的顺序输出.

分析:本题与上题基本上一样,只是求的量不同,既然我们在找树的重心的时候用的树型dp,而且是求的子树中节点数的最大值,然后求所有最

大值的最小值,那么就有可能存在多个重心,我们每更新到一个最小值的时候就记录其它的最小值也为这个最小值的重心,这样下去就会找到所

有的重心.

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
using namespace std;
const int MAXN = 50005;
const int INF = 0x3f3f3f3f;
int n;
int head[MAXN], tol;
int son[MAXN], ans[MAXN], siz, num;
bool vis[MAXN];
struct Edge
{
    int to, nxt;
}edge[MAXN*2];
void addedge(int u, int v)
{
    edge[tol].to = v;
    edge[tol].nxt = head[u];
    head[u] = tol++;
}
void dfs(int cur)
{
    son[cur] = 0;
    vis[cur] = true;
    int tmp = 0;
    for(int i = head[cur]; i != -1; i = edge[i].nxt)
    {
        int v = edge[i].to;
        if(!vis[v])
        {
            dfs(v);
            son[cur] += son[v] + 1;
            tmp = max(tmp, son[v]+1);
        }
    }
    tmp = max(tmp, n-son[cur]-1);
    if(tmp < siz)
    {
        num = 1;
        ans[0] = cur;
        siz = tmp;
    }
    else if(tmp == siz)
        ans[num++] = cur;
}
void init()
{
    siz = INF;
    tol = num = 0;
    memset(vis, false, sizeof(vis));
    memset(head, -1, sizeof(head));
}
int main()
{
    while(~scanf("%d", &n))
    {
        init();
        for(int i = 1; i <= n-1; i++)
        {
            int u, v;
            scanf("%d%d", &u, &v);
            addedge(u, v);
            addedge(v, u);
        }
        dfs(1);
        sort(ans, ans+num);
        for(int i = 0; i < num; i++)
            printf("%d%c", ans[i], i == num-1 ? '\n':' ');

    }
    return 0;
}