BZOJ1017

传送门：BZOJ1017

比较复杂的树形Dp。
复杂的是这道题Dp方程的定义:$f_{i,j,k}$表示以节点i为根的子树，向i的父节点提供j个物品，在子树上总花费k元能得到的最大力量值。

转移时我们先不考虑以i为根的子树合成了n个i却未完全上交的情况。i为树叶时，Dp是显然的，如果i不为树叶，则转换是显然的。但在这转移中，假设我们给子节点$j'$分配了$k'$元，则我们要考虑其它所有的情况。

避免的方法是，用g数组来记录，这一步不太好讲但是比较好懂，建议直接看代码。

代码上的小细节见下。

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

const int maxlongint=0x7fffffff,inf=maxlongint/3;

int n,m;
int power[105];
int num[105];
int cost[105];
int need[105],fa[105];
int f[105][105][2205];
int g[2205];
int ans;
vector<int> G[105];

void Dfs(int p)
{
    if(G[p].size()==0){
        num[p]=min(num[p],m/cost[p]);
        for(int i=0;i<=num[p];i++)
            for(int j=i;j<=num[p];j++)
                f[p][i][j*cost[p]]=(j-i)*power[p];
        return;
    }
    num[p]=inf;
    for(int i=0;i<G[p].size();i++){
        int v=G[p][i];
        Dfs(v);
        num[p]=min(num[p],num[v]/need[v]);
    }
    for(int i=0;i<=num[p];i++)
        f[p][i][0]=0;  
    for(int i=0;i<G[p].size();i++){
        int v=G[p][i];
        for(int j=0;j<=num[p];j++){
            memcpy(g,f[p][j],sizeof(f[p][j]));
            memset(f[p][j],-1,sizeof(f[p][j]));
            for(int k=0;k<=m;k++)
                for(int r=0;r<=k;r++)
                    if(g[k-r]!=-1&&f[v][j*need[v]][r]!=-1)
                        f[p][j][k]=max(f[p][j][k],g[k-r]+f[v][j*need[v]][r]);
        }
    }
    for(int i=0;i<=num[p];i++)
        for(int j=i;j<=num[p];j++)
            for(int k=0;k<=m;k++)
                if(f[p][j][k]!=-1)
                    f[p][i][k]=max(f[p][i][k],f[p][j][k]+(j-i)*power[p]),ans=max(ans,f[p][i][k]);
}

void Readdata()
{
    freopen("loli.in","r",stdin);
    scanf("%d%d",&n,&m);
    char crash;
    for(int i=1;i<=n;i++){
        scanf("%d%c",&power[i],&crash);
        scanf("%c",&crash);
        if(crash=='B')
            scanf("%d%d",&cost[i],&num[i]);
        else{
            int a;
            scanf("%d",&a);
            for(int j=1;j<=a;j++){
                int b,c;
                scanf("%d%d",&b,&c);
                G[i].push_back(b);
                need[b]=c;fa[b]=i;
            }
        }
    }
}

void Solve()
{
    memset(f,-1,sizeof(f));
    for(int i=1;i<=n;i++)
        if(!fa[i])
            Dfs(i);
    printf("%d\n",ans);
}

void Close()
{
    fclose(stdin);
    fclose(stdout);
}

int main()
{
    Readdata();
    Solve();
    Close();
    return 0;
}