1. 程式人生 > >[POI2005]BAN-Bank Notes [多重揹包]

[POI2005]BAN-Bank Notes [多重揹包]

傳送門

二進位制拆分 , 然後f[j] 表示到j的最少錢幣 

f[j]=min(f[j],f[j-v[i]]+w[i])

用一個bool陣列記錄是否轉移


#include<bits/stdc++.h>
#define N 205
#define M 20050
using namespace std;
int n,f[M],b[N],c[N],k,tot;
struct Node{int w,v,val;}a[N*16];
bool from[N*16][M]; int ans[M],vis[M];
void dfs(int now,int val){
    if(now==0 || val==0) return;
    while(!from[now][val]) now--;
    val -= a[now].val;
    ans[a[now].w] += a[now].v;
    now--; dfs(now,val);
}
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++) scanf("%d",&b[i]),vis[b[i]]=1;
    for(int i=1;i<=n;i++){
        scanf("%d",&c[i]); int x=c[i];
        for(int k=1;k<=x;k<<=1)
            a[++tot]=(Node){b[i],k,b[i]*k},x-=k;
        if(x) a[++tot]=(Node){b[i],x,b[i]*x};
    }
    memset(f,127,sizeof(f)); f[0]=0;
    scanf("%d",&k);
    for(int i=1;i<=tot;i++){
        for(int j=k;j>=a[i].val;j--){
            if(f[j-a[i].val] + a[i].v < f[j]){
                f[j] = f[j-a[i].val] + a[i].v;
                from[i][j] = 1;
            }
        }
    }
    printf("%d\n",f[k]); dfs(tot,k);
    for(int i=1;i<=20000;i++) if(vis[i]) printf("%d ",ans[i]);
    return 0;
}