树形DP入门

给定一颗有N个节点的树（一般是无根树，就有N-1条无向边），可以任选一个节点作为根节点

一般以节点从深到浅（子树从小到大）的顺序作为dp阶段顺序

dp的状态表示中，第一维通常是节点编号（节点编号代表了以该节点为根的子树）

对于每个节点x，先递归在它的每个子节点上进行dp，回溯时，从子节点向x进行状态转移

A - Anniversary part

N个员工，编号为1~N

他们之间有从属关系，也就是说他们的关系就像一棵以校长为根的树，父结点就是子结点的直属上司。现在有个宴会，宴会每邀请来一个员工 i 都会增加一定的快乐指数 Ri，但如果某个员工的直属上司来了，那么这个员工就不会来。计算邀请哪些员工可以使快乐指数最大，输出最大的快乐指数

设 dp [x] [0] 表示在 x 为根的子树中邀请部分员工，并且 x 不参加时，快乐指数总和的最大值，此时 x 子节点（直接下属）可以参加也可以不参加 (s表示x子节点)

设 dp [x] [1] 表示在 x 为根的子树中邀请部分员工，并且 x 参加时，快乐指数总和的最大值，此时 x 子节点（直接下属）都不能参加，H[x] 表示当前节点（x）的快乐指数 (s表示x子节点)

这个题给的是有根树，就可以从根节点开始，dp目标就是 max(F[root,0] , F[root,1]) 时间复杂度ON

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#include <vector>
using namespace std;

vector<int> son[10010];
int dp[10010][2];//0不参加，1参加
int v[10010];//记录有没有父节点
int h[10010];//记录快乐指数
int n;

void DFS(int x){
    dp[x][0] = 0;
    dp[x][1] = h[x];
    for (int i = 0; i < son[x].size(); i++)
    {
        int y = son[x][i];
        DP(y);
        dp[x][0] += max(dp[y][0],dp[y][1]);
        dp[x][1] += dp[y][0];
    }
}

int main(){
    cin>>n;
    for (int i = 1; i <=n ; i++)
        scanf("%d",&h[i]);
    for (int i = 1; i <n ; i++)
    {
        int x,y;
        cin>>x>>y;
        v[x] = 1;//x有爸爸
        son[y].push_back(x);//x是y的子节点
    }
    int root;
    for (int i = 1; i <= n; i++)
        if(!v[i]){ //i没有爸爸
            root = i;
            break;
        }
    DFS(root);
    cout << max(dp[root][0],dp[root][1]) << endl;
    return 0;
}

B - Strategic game

有n个点，在某些点上放置哨兵，每个哨兵可以监控和它有边相连的点，问监视所有的点需要的最少的哨兵数

也就是说一颗n个结点的树，要求选出其中的一些顶点，使得对于树中的每条边(u, v)，u和v至少有一个被选中，要求给出选中顶点数最少的方案

设 dp [x] [0] 表示在不选择节点 x 的情况下，以 x 为根节点的子树，最少需要选择的节点数

​ 当i为叶子节点时

​ 当i不为叶子节点时 (s表示x子节点)

设 dp [x] [1] 表示在选择节点 x 的情况下，以 x 为根节点的子树，最少需要选择的节点数
​ 当i为叶子节点时

​ 当i不为叶子节点时 (s表示x子节点)

#include <iostream>
#include <cstring>
#include <cstdio>
#include <cmath>
#define maxn 1508
using namespace std;

int dp[maxn][2];
int soncnt[maxn];
int parent[maxn];
int n;

void DFS(int x) {
    int i, d1=0, d0=0;
    if (soncnt[x] == 0) {
        dp[x][0] = 0;
        dp[x][1] = 1;
        return;
    }
    for (i=0; i < n; i++) {
        if (parent[i] == x) {
            DFS(i);
            d1 += min(dp[i][0], dp[i][1]);
            d0 += dp[i][1];
        }
    }
    dp[x][1] = d1 + 1;
    dp[x][0] = d0;
}

int main() {
    int dad, son, m;
    while (cin >> n) {
        memset(soncnt, 0, sizeof(soncnt));
        memset(parent, -1, sizeof(parent));
        int root = -1;
        for (int i = 0; i < n; i++) {
            scanf("%d:(%d)", &dad, &m);
            soncnt[dad] = m;
            if (root == -1) {
                root = dad;
            }
            while (m--) {
                scanf("%d", &son);
                parent[son] = dad;
            }
        }
        DFS(root);
        cout << min(dp[root][0], dp[root][1]) << endl;
    }
    return 0;
}

C - Tree Cutting

给一颗n个结点的树，节点编号为1~n

问：删除哪些结点后，剩余各个子树的大小均小于原总结点数的一半

拆除一个节点后，剩余部分为其若干儿子的子树以及该节点上层所连其余部分（n-size[i]），只要这些连接块大小都不超过n/2，该节点就满足条件。因而我们可以先求出每个节点所管辖的那棵树的大小，自下而上地为每个节点求出其子树规模（该点规模=其儿子的规模和+1）。

在建树的时候可以直接用连接表（vector）储存无向边，这时由于无法区分与每个点相连的是其父节点还是子节点，会引发问题：在dfs的时候把父节点误认作儿子节点，解决方法就是在递归的时候传入父节点编号，然后在递归，计算规模，比较大小的时候避开它就可以了

#include<iostream>
#include<cstring>
#include<vector>
#include<cstdio>
#include<algorithm>
using namespace std;
 
int n;
int root;
vector<int>G[10000+400];
vector<int>ans;
int sz[10000+400];
 
void dfs(int par,int u){
    sz[u]=1;
    for(int i=0;i<(int)G[u].size();i++){
        if(G[u][i]!=par)
            dfs(u,G[u][i]);
    }
    int piece=0;
    for(int i=0;i<(int)G[u].size();i++)
        if(par!=G[u][i]){
            sz[u]+=sz[G[u][i]];
            piece=max(piece,sz[G[u][i]]);
        }
    piece=max(piece,n-sz[u]);
    if(piece<=n/2)
        ans.push_back(u);
}
 
int main(){
    while(cin>>n){
        memset(sz,0,sizeof(sz));
        int x,y;
        for(int i=0;i<n-1;i++){
            scanf("%d%d",&x,&y);
            G[x].push_back(y);
            G[y].push_back(x);
        }
        dfs(0,1);
        sort(ans.begin(),ans.end());
        for(int i=0;i<(int)ans.size();i++)
            cout<<ans[i]<<endl;
    }
    return 0;
}

LCA 最近公共祖先

倍增（基于二分的方法）

假如树上的两个点，处于同一深度时，对于往上跳mid步，显然有单调性：

如果它们往上跳mid步是同一个点，则这个点是它们的共同祖先，但不一定是最近公共祖先

那么如果我们可以快速得到每个点往上跳若干步是谁，显然我们可以利用二分来求LCA

这个可以用倍增来做。

倍增的思想其实非常简单，就是利用二进制的思想。一个显然的结论有：

节点u往上跳2^(k+1)的布的祖先 = （节点u往上跳2^k布的祖先）往上跳2^k布的祖先

用代码表示就是：

fa[u][k + 1] = fa[ fa[u][k] ][k];

这个的预处理也非常简单，我们只需要初始化每个fa[u][0]就可以了就可以了，因为接下来的部分都可以根据上面的公式递推。

然后预处理好这个之后，我们就可以根据这个进行二分了，代码如下：

int fa[maxn][20], dep[maxn];

void dfs(int u, int f, int d)
{
    dep[u] = d;
    fa[u][0] = f;
    for(int i = 1; i < 20; ++i)
    {
        fa[u][i] = fa[fa[u][i -  1]][i] - 1;
    }
    for(int i = head[u]; ~i; i = nxt[i])
    {
        int t = to[i];
        if(t != f)
        {
            dfs(t, u, d + 1);
        }
    }
}

int lca(int u, int v)
{
    if(dep[u] < dep[v]) swap(u, v);
    int k = dep[u] - dep[v];
    for(int i = 0; i < k; ++i)
    {
        if((1 << i) & k) u = fa[u][i];
    }
    if(u == v) return u;
    for(int i = 19; i >= 0; --i)
    {
        if(fa[u][i] != fa[v][i])
        {
            u = fa[u][i];
            v = fa[v][i];
        }
    }
    return fa[u][0];
}