POJ1947 Rebuilding Roads（树形DP）

题目大概是给一棵树，问最少删几条边可以出现一个包含点数为p的连通块。

任何一个连通块都是某棵根属于连通块的子树的上面一部分，所以容易想到用树形DP解决：

dp[u][k]表示以u为根的子树中，包含根的大小k的连通块最少的删边数
要求答案就是min(dp[u][p],min(dp[v][p]+1))，u是整棵树的根，v是其他结点
转移从若干个子树各自选择要提供几个k转移，不过指数级时间复杂度，当然又是树上背包了。。

转移好烦，写得我好累好累。。还好1A了。。

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<algorithm>
 4 using namespace std;
 5 #define INF (1<<29)
 6 #define MAXN 155
 7 struct Edge{
 8     int u,v,next;
 9 }edge[MAXN];
10 int NE,head[MAXN];
11 void addEdge(int u,int v){
12     edge[NE].u=u; edge[NE].v=v; edge[NE].next=head[u];
13     head[u]=NE++;
14 }
15 int d[MAXN][MAXN],size[MAXN],son[MAXN];
16 void getSize(int u){
17     size[u]=1;
18     son[u]=0;
19     for(int i=head[u]; i!=-1; i=edge[i].next){
20         int v=edge[i].v;
21         getSize(v);
22         size[u]+=size[v];
23         ++son[u];
24     }
25 }
26 void dp(int u){
27     d[u][size[u]]=0;
28     d[u][1]=0;
29     bool first=1;
30     for(int i=head[u]; i!=-1; i=edge[i].next){
31         int v=edge[i].v;
32         dp(v);
33         ++d[u][1];
34         if(first){
35             for(int j=1; j<=size[v]; ++j) d[u][j+1]=d[v][j];
36             first=0;
37         }else{
38             for(int j=size[u]-2; j>=1; --j){
39                 ++d[u][j+1];
40                 for(int k=1; k<=min(j,size[v]); ++k){
41                     d[u][j+1]=min(d[u][j+1],d[v][k]+d[u][j+1-k]);
42                 }
43             }    
44         }
45     }
46 }
47 int main(){
48     for(int i=0; i<MAXN; ++i){
49         for(int j=0; j<MAXN; ++j) d[i][j]=INF;
50     }
51     memset(head,-1,sizeof(head));
52     int n,p,a,b;
53     scanf("%d%d",&n,&p);
54     int deg[MAXN]={0};
55     for(int i=1; i<n; ++i){
56         scanf("%d%d",&a,&b);
57         addEdge(a,b);
58         ++deg[b];
59     }
60     int root;
61     for(int i=1; i<=n; ++i){
62         if(deg[i]==0) root=i;
63     }
64     getSize(root);
65     dp(root);
66     int res=INF;
67     for(int i=1; i<=n; ++i){
68         if(root==i) res=min(res,d[i][p]);
69         else res=min(res,d[i][p]+1);
70     }
71     printf("%d",res);
72     return 0;
73 }