BZOJ1017: [JSOI2008]魔兽地图DotR

Description

DotR (Defense of the Robots) Allstars是一个风靡全球的魔兽地图，他的规则简单与同样流行的地图DotA (Defense of the Ancients) Allstars。DotR里面的英雄只有一个属性——力量。他们需要购买装备来提升自己的力量值，每件装备都可以使佩戴它的英雄的力量值提高固定的点数，所以英雄的力量值等于它购买的所有装备的力量值之和。装备分为基本装备和高级装备两种。基本装备可以直接从商店里面用金币购买，而高级装备需要用基本装备或者较低级的高级装备来合成，合成不需要附加的金币。装备的合成路线可以用一棵树来表示。比如，Sange and Yasha的合成需要Sange, Yasha和Sange and Yasha Recipe Scroll三样物品。其中Sange又要用Ogre Axe, Belt of Giant Strength 和 Sange Recipe Scroll合成。每件基本装备都有数量限制，这限制了你不能无限制地合成某些性价比很高的装备。现在，英雄Spectre有M个金币，他想用这些钱购买装备使自己的力量值尽量高。你能帮帮他吗？他会教你魔法Haunt（幽灵附体）作为回报的。

Input

输入文件第一行包含两个整数，N (1 <= n <= 51) 和 m (0 <= m <= 2,000)。分别表示装备的种类数和金币数。装备用1到N的整数编号。接下来的N行，按照装备1到装备n的顺序，每行描述一种装备。每一行的第一个正整数表示这个装备贡献的力量值。接下来的非空字符表示这种装备是基本装备还是高级装备，A表示高级装备，B表示基本装备。如果是基本装备，紧接着的两个正整数分别表示它的单价（单位为金币）和数量限制（不超过100）。如果是高级装备，后面紧跟着一个正整数C，表示这个高级装备需要C种低级装备。后面的2C个数，依次描述某个低级装备的种类和需要的个数。

Output

第一行包含一个整数S，表示最多可以提升多少点力量值。

Sample Input

10 59
5 A 3 6 1 9 2 10 1
1 B 5 3
1 B 4 3
1 B 2 3
8 A 3 2 1 3 1 7 1
1 B 5 3
5 B 3 3
15 A 3 1 1 5 1 4 1
1 B 3 5
1 B 4 3

Sample Output

用传统的办法。。。

设f[i][j][k]表示以i为根的子树，要提前合成j个道具i提供给上面的节点，花费代价为k时最多可以提升多少点力量值。

那么我们转移时枚举一共要合成几个道具i。然后再做个背包，大概是前i个节点花费代价为k最多可以提升多少点力量值。

然后因为打错了一个字母竟然调了2h不能忍。

#include<cstdio>
#include<cctype>
#include<queue>
#include<cmath>
#include<cstring>
#include<algorithm>
#define rep(i,s,t) for(int i=s;i<=t;i++)
#define dwn(i,s,t) for(int i=s;i>=t;i--)
#define ren for(int i=first[x];i;i=next[i])
using namespace std;
inline int read() {
    int x=0,f=1;char c=getchar();
    for(;!isdigit(c);c=getchar()) if(c=='-') f=-1;
    for(;isdigit(c);c=getchar()) x=x*10+c-'0';
    return x*f;
}
const int maxn=55;
const int maxm=2005;
const int maxp=105;
const int inf=1e9;
char ch[12];
int n,m,first[maxn],next[maxm*10],to[maxm*10],dis[maxm*10],in[maxn],e;
void AddEdge(int u,int v,int w) {
    dis[++e]=w;to[e]=v;next[e]=first[u];first[u]=e;in[v]++;
}
int f[maxn][maxp][maxm],g[maxn][maxm],p[maxn],c[maxn],s[maxn];
void dp(int x) {
    if(!first[x]) {
        s[x]=min(s[x],m/c[x]);//QAQ,打成m/p[x]了2333
        rep(i,0,s[x]) rep(j,i,s[x]) f[x][i][j*c[x]]=(j-i)*p[x];
        return;
    }
    s[x]=inf;
    for(int i=first[x];i;i=next[i]) {
        dp(to[i]);
        s[x]=min(s[x],s[to[i]]/dis[i]);
        c[x]+=c[to[i]]*dis[i];
    }
    s[x]=min(s[x],m/c[x]);
    memset(g,-0x3f3f3f3f,sizeof(g));g[0][0]=0;
    dwn(y,s[x],0) {
        int cnt=0;
        for(int i=first[x];i;i=next[i]) {
            cnt++;
            dwn(j,m,0) rep(k,0,j) g[cnt][j]=max(g[cnt][j],g[cnt-1][j-k]+f[to[i]][y*dis[i]][k]);
        }    
        rep(j,0,y) rep(k,0,m) f[x][j][k]=max(f[x][j][k],g[cnt][k]+p[x]*(y-j));
    }
}
int main() {
    memset(f,-0x3f3f3f3f,sizeof(f));
    n=read();m=read();
    rep(i,1,n) {
        p[i]=read();scanf("%s",ch);
        if(ch[0]=='A') {
            int k=read(),v,w;
            while(k--) v=read(),w=read(),AddEdge(i,v,w);
        }
        else c[i]=read(),s[i]=read();
    }
    int ans=0;
    rep(i,1,n) if(!in[i]) {
        dp(i);rep(j,0,m) ans=max(ans,f[i][0][j]);
    }
    printf("%d
",ans);
    return 0;
}

View Code