poj1742coins【单调队列优化多重背包】楼教主八题

早上好不容易想明白,就是不知道自己哪里错了==

我的错的

/***************
poj1742
2016.2.24
不知道哪里错了啊啊啊啊啊

***************/
#include <iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int n,m,c[102],a[102],q[102];
bool dp[100006];
int main()
{
   // freopen("cin.txt","r",stdin);
    while(~scanf("%d%d",&n,&m))
    {
        if(n==0&&m==0) break;
        for(int i=1;i<=n;i++) scanf("%d",&a[i]);
        for(int i=1;i<=n;i++) scanf("%d",&c[i]);
        int   sum=0;
        memset(dp,0,sizeof(dp));
        dp[0]=1;
        for(int i=1;i<=n;i++)
        {
            if(c[i]==1)
                for(int v=m;v>=a[i];v--)
                    if(dp[v-a[i]]) dp[v]=1;
            else if(a[i]*c[i]>=m)
                for(int v=a[i];v<=m;v++)
                    if(dp[v-a[i]]) dp[v]=1;
            else
            {
                for(int d=0;d<a[i];d++)
                {
                    int st=0,ed=-1;
                    sum=0;
                    for(int v=d;v<=m;v+=a[i])
                    {
                        if(ed-st==c[i]) sum-=q[st++];
                        q[++ed]=dp[v];
                        sum+=dp[v];
                        if(!dp[v]&&sum) dp[v]=1;
                    }
                }
            }
        }
        sum=0;
        for(int i=1;i<=m;i++) if(dp[i]) sum++;
        printf("%d\n",sum);
    }
    return 0;
}

ac的

#include <iostream>
#include<cstring>
using namespace std;
 
const int MAXN = 150;
int w[MAXN], c[MAXN];
int n, m;
bool dp[100000+10], queue[100000+10];
 
int solve_dp() {
    memset(dp, 0, sizeof(dp));
    memset(queue, 0, sizeof(queue));
    dp[0] = true;
 
    for(int i = 0; i < n; i ++) {
        if(c[i] == 1) {     // 仅允许一个包, 变成01背包问题
            for(int v = m; v >= w[i]; v--) {
                if(!dp[v] && dp[v-w[i]])
                    dp[v] = 1;
            }
        }else if(c[i]*w[i] >= m) {   // 完全背包问题, 即 w[i]*c[i] < m, 放入件数的限制是 c[i]
            for(int v = w[i]; v <= m; v++) {
                if(!dp[v] && dp[v-w[i]])
                    dp[v] = 1;
            }
        }else{  // 严格的分组背包问题
            for(int d = 0; d < w[i]; d++) {      // 对于所有余数 d [0, w[i])
                // 窗口大小为 c[i]
                int sum = 0, st = 0, ed = -1;   //st,ed 单调队列的开始和结尾, sum 队列中是否有一个 true
                for(int v = d; v <= m; v+= w[i]) {   // 完全背包 model, 但步长是 w[i]
                    if(ed - st == c[i]) {   // 窗口大小为0, 移除队首元素, 队首后移一位
                        sum -= queue[st++];
                    }
                    queue[++ed] = dp[v];
                    sum += dp[v];
                    if(!dp[v] && sum)
                        dp[v] = 1;
                }
            }
        }
    }
    int res = 0;
    for(int i = 1; i <= m; i ++)
        res += dp[i];
    return res;
 
}
int main() {
   // freopen("E:\\Copy\\ACM\\测试用例\\in.txt", "r", stdin);
    while(cin >> n >> m && n != 0) {
        for(int i = 0; i < n; i ++) {
            scanf("%d", &w[i]);
        }
        for(int i = 0; i < n; i ++) {
            scanf("%d", &c[i]);
        }
        // main function
        cout << solve_dp() << endl;
    }
    return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值