BZOJ1017

傳送門：BZOJ1017

比較複雜的樹形Dp。
複雜的是這道題Dp方程的定義:fi,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;
}