[bzoj1017]:[JSOI2008]魔兽地图DotR

传送门
JSOI2008 All Clear!
这个题简直了,我整个代码都是抄的,抄完还是一脸懵逼
我就tm没见过这么难的树形dp
一开始还以为水题,然后发现根本不知道怎么dp
题解还是直接抄hzwer的吧。。

设P[x],L[x],M[x],表示物品x的能量,购买上限与价格
L[x]=min(L[x],m/M[x])
高级装备的M和L随便dp一下
再用f[i][j][k]表示第i个物品,有j件用于上层的合成,花费金钱是k所能获得的最大力量
对于以x为根的子树,枚举合成 l 个 x 物品,然后再用其余的钱买一些 x 子树内的装备不用于合成
枚举合成物品数量 l ,用g[i][j]表示x的前i个儿子的子树,花费j的钱,所能获得的最大力量
g[tot][j]=max{g[tot-1][j-k]+f[e[i].to][l*e[i].v][k]}
//e[i].to是儿子结点,e[i].v是所需数量,就是从j中拿出k的钱在e[i].to的子树内购买
最后再枚举合成的 l 个 x 物品中有 j 个是直接用于增加力量,剩余用于合成的
f[x][j][k]=max{g[tot][k]+P[x]*(l-j)}

代码:

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<cstdlib>
#define ll long long
using namespace std;
inline int read(){
    int x=0;char ch=' ';int f=1;
    while(ch!='-'&&(ch<'0'||ch>'9'))ch=getchar();
    if(ch=='-')f=-1,ch=getchar();
    while(ch>='0'&&ch<='9')x=(x<<3)+(x<<1)+(ch^48),ch=getchar();
    return x*f;
}
const int N=55,M=2005,K=105;
struct edge{int to,next,l;}e[M];
int n,m,tot;
int head[N],f[N][K][M],g[N][M],h[N][M];
inline void addedge(int x,int y,int l){e[++tot].to=y;e[tot].next=head[x];head[x]=tot;e[tot].l=l;}
int v[N],cost[N],mx[N],in[N],out[N];
inline void dp(int x){
    if(!out[x]){
        mx[x]=min(mx[x],m/cost[x]);
        for(int i=0;i<=mx[x];++i)
            for(int j=i;j<=mx[x];++j)
                f[x][i][j*cost[x]]=(j-i)*v[x];
        return;
    }
    mx[x]=0x3f3f3f3f;
    for(int i=head[x];i;i=e[i].next){
        int u=e[i].to;
        dp(u);
        mx[x]=min(mx[x],mx[u]/e[i].l);
        cost[x]+=cost[u]*e[i].l;
    }
    mx[x]=min(mx[x],m/cost[x]);
    memset(g,-0x3f,sizeof(g));
    g[0][0]=0;
    for(int l=mx[x];l>=0;--l){
        int tot=0;
        for(int i=head[x];i;i=e[i].next){
            tot++;
            for(int j=0;j<=m;++j)
                for(int k=0;k<=j;++k)
                    g[tot][j]=max(g[tot][j],g[tot-1][j-k]+f[e[i].to][l*e[i].l][k]);
        }
        for(int j=0;j<=l;++j)
            for(int k=0;k<=m;++k)
                f[x][j][k]=max(f[x][j][k],g[tot][k]+v[x]*(l-j));
    }
}
int main(){
    memset(f,-0x3f,sizeof(f));
    n=read();m=read();
    for(int i=1;i<=n;++i){
        char ch[3];
        v[i]=read();
        scanf("%s",ch);
        if(ch[0]=='A'){
            int C=read();
            out[i]=C;
            while(C--){
                int to=read(),num=read();
                addedge(i,to,num);
                in[to]++;
            }
        }
        else cost[i]=read(),mx[i]=read();
    }
    int tot=0;
    for(int x=1;x<=n;++x){
        if(!in[x]){
            tot++;
            dp(x);
            for(int i=0;i<=m;++i)
                for(int j=0;j<=i;++j)
                    for(int k=0;k<=mx[x];++k)
                        h[tot][i]=max(h[tot][i],h[tot-1][j]+f[x][k][i-j]);
        }
    }
    int ans=0;
    for(int i=0;i<=m;++i)ans=max(ans,h[tot][i]);
    printf("%d",ans);
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值