dp学习笔记 6 没有上司的舞会2 (树形dp)

题目
题意: n个点的有根树，每个点有权值ai。选了父节点就不能选择子节点，并且最多选择m个节点。
思路: f[i][j][0/1]:以i为根，至多选择j个节点，0：选，1：不选。
参考上一题的什么背包，枚举树根i选择的节点个数j，之后枚举子树选择的节点0,1…j。如果不选树根i，那么子树都可以选；如果选树根i，只能不选子树的根节点。
 这里有点费解的是不知道为啥j要从大到小枚举，wls的解释是因为该题的模型和上一题背包一样，所以照搬从大到小。
 以后最后赋值的循环来看，如果从小到大会导致一个点会选择多次，每个i都能被i-1更新，i=1时正是选择它自己，i=2又能再选一次自己，寄。同样，递推方程时也有点这个意思，会重复计算的，学不明白了。
时间复杂度: O(nmm)
代码:

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<complex>
#include<cstring>
#include<cmath>
#include<vector>
#include<map>
#include<unordered_map>
#include<list>
#include<set>
#include<queue>
#include<stack>
#define OldTomato ios::sync_with_stdio(false),cin.tie(nullptr),cout.tie(nullptr)
#define fir(i,a,b) for(int i=a;i<=b;++i) 
#define mem(a,x) memset(a,x,sizeof(a))
#define p_ priority_queue
// round() 四舍五入 ceil() 向上取整 floor() 向下取整
// lower_bound(a.begin(),a.end(),tmp,greater<ll>()) 第一个小于等于的
// #define int long long //QAQ
using namespace std;
typedef complex<double> CP;
typedef pair<int,int> PII;
typedef long long ll;
// typedef __int128 it;
const double pi = acos(-1.0);
const int INF = 0x3f3f3f3f;
const ll inf = 1e18;
const int N = 502;
const int M = 1e6+10;
const int mod = 1e9+7;
const double eps = 1e-6;
inline int lowbit(int x){ return x&(-x);}
template<typename T>void write(T x)
{
    if(x<0)
    {
        putchar('-');
        x=-x;
    }
    if(x>9)
    {
        write(x/10);
    }
    putchar(x%10+'0');
}
template<typename T> void read(T &x)
{
    x = 0;char ch = getchar();ll f = 1;
    while(!isdigit(ch)){if(ch == '-')f*=-1;ch=getchar();}
    while(isdigit(ch)){x = x*10+ch-48;ch=getchar();}x*=f;
}
int n,m,k,T;
int a[N];
vector<int> va[N];
int f[N][N][2]; //以i为根、子树中选择的节点不超过j个。
void dfs(int cur)
{
	for(int son:va[cur])
	{
		dfs(son);
		for(int j=m;j;--j)
		for(int k=1;k<=j;++k)
		{
			f[cur][j][0] = max(f[cur][j][0],f[cur][j-k][0] + f[son][k][1]);
			f[cur][j][1] = max(f[cur][j][1],f[cur][j-k][1] + max(f[son][k][1],f[son][k][0]));
		}
	}
	for(int i=m;i;--i) f[cur][i][0] = f[cur][i-1][0] + a[cur];
}
void solve()
{
    cin>>n>>m;
    for(int i=2;i<=n;++i){int x;cin>>x;va[x].push_back(i);}
    for(int i=1;i<=n;++i) cin>>a[i];
    dfs(1);
    cout<<max(f[1][m][0],f[1][m][1]); 
}
signed main(void)
{  
   T = 1;
   // OldTomato; cin>>T;
   // read(T);
   while(T--)
   {
   	 solve();
   }
   return 0;
}