poj3659 ( 树形dp )

参考博客：https://blog.csdn.net/jiang199235jiangjj/article/details/7878473

动态规划的思路：
         根据题意知道每个节点有三种状态：
              1、点i建塔，i的所有孩子都覆盖，用dp[i][0]表示；
              2、点i不建塔，i和i的所有孩子都覆盖，用dp[i][1]表示；
              3、点i不建塔，i不覆盖，i的所有孩子都覆盖，用dp[i][2]表示；
         如果这样不好理解那么这样理解可能容易一点（参考别人的）：覆盖i，要么是父节点覆盖，要么是自己，要么是孩子，所以三种状态（和上面的对应）：             
              1、点i自己覆盖自己，i的所有孩子都覆盖，用dp[i][0]表示；
              2、点i被自己的孩子覆盖，i和i的所有孩子都覆盖，用dp[i][1]表示；
              3、点i被父节点覆盖，i的所有孩子都覆盖，用dp[i][2]表示；
         那么动态转移方程就是（v是i的孩子）：
              dp[i][0]+=min(dp[v][0],dp[v][1],dp[v][2]);
              dp[i][2]+=min(dp[v][0],dp[v][1]);
              对于dp[i][1],要考虑全面，也就是说：必须要有一个孩子建塔，才能保证i被覆盖（Min=sum(min(dp[v][0]-dp[i][1]))，也就是当所有孩子的dp[v][0]>dp[v][i]时，Min表示他们差值最小的那个差值）。
              所以方程是dp[i][1]+=min(dp[v][0],dp[1])（至少存在一个孩子的dp[v][0]<=dp[v][1]，否则要dp[i][1]+=Min）；

AC代码：

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<string>
#include<queue>
#include<vector>
#include<cmath>
#include<stack>
#include<map>
#include<set>
using namespace std;
#define LL long long
const int MOD=100000007;
const int inf=0x3f3f3f3f;
const LL inff=0x3f3f3f3f3f3f3f3f;
const LL N=300005;
const LL M=5005;
#define MEF(x) memset(x,-1,sizeof(x))
#define ME0(x) memset(x,0,sizeof(x))
#define MEI(x) memset(x,inf,sizeof(x))
int dp[10005][3];
bool vis[10005];
void Dp(int x)
{
    bool flag=true;
    vis[x]=false;
    int Min=(1<<30);
    dp[x][0]=1;
    dp[x][1]=dp[x][2]=0;
    for(int i=first[x];i!=-1;i=edge[i].next)
    {
        int lv=edge[i].v;
        if(vis[lv])
        {
            Dp(lv);
            dp[x][0]+=min(dp[lv][0],min(dp[lv][1],dp[lv][2]));
            dp[x][2]+=min(dp[lv][0],dp[lv][1]);
            if(dp[lv][0]<=dp[lv][1])
            {
                flag=false;
                dp[x][1]+=dp[lv][0];
            }
            else
            {
                dp[x][1]+=dp[lv][1];
                Min=min(Min,dp[lv][0]-dp[lv][1]);
            }
        }
    }
    if(flag)
    {
        dp[x][1]+=Min;
    }
}
int main()
{
    ios::sync_with_stdio(false);
    int n;
    cin>>n;
    init();
    for(int n1=1,a,b;n1<=n-1;n1++)
    {
        cin>>a>>b;
        add(a,b);
        add(b,a);
    }
    ME0(dp);
    memset(vis,true,sizeof(vis));
    Dp(1);
    cout<<min(dp[1][0],dp[1][1])<<endl;
    return 0;
}