The Balance(母函数或01背包)

本文探讨了如何通过01背包问题和母函数方法解决特定的砝码称量问题,即找出在给定砝码数量及重量的情况下,无法称量出的所有重量值。文中提供了两种算法的具体实现代码。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题目链接:The Balance

题目大意:给你n个砝码,问在砝码和的范围内,有多少个数值是无法称量出来的。
注意:砝码可以放在称量物一边,例如:4和5,能称量出1。

思路:
首先我会想到01背包,其次是母函数。

母函数思路:
这里我们要在创造一些砝码,就是他们的负砝码。
给你 m1,m2,m3,m4….mn,创造出-m1,-m2,-m3,-m4….-,mn.
这样我们就可以写母函数方程了。
f(x)=(1+x^m1)(1+x^m2)(1+x^m3)…(1+x^mn)(1+x^(-m1))(1+x^(-m2))(1+x^(-m3))…(1+x^(-mn)).
计算过程中指数出现负数的就不用记录了。

代码:

#include<stdio.h>
#include<string.h>

const int maxn=10009;
int a[maxn],dp[maxn],dp1[maxn],cun[maxn];
int main()
{
    int n;
    while(~scanf("%d",&n))
    {
        int N=2*n,sum=0;
        for(int i=0;i<n;i++)
            scanf("%d",&a[i]),sum+=a[i],a[i+n]=-a[i];//记录砝码和,创造负砝码
        memset(dp,0,sizeof(dp));
        dp[0]=1,dp[a[0]]=1;
        for(int i=1;i<N;i++)
        {
            for(int j=0;j<=sum;j++)
            {
                for(int k=0;k<=1;k++)
                {
                    if(j+k*a[i]>=0)//只记录指数为非负的即可
                        dp1[j+k*a[i]]+=dp[j];
                }
            }
            for(int j=0;j<=sum;j++)
            {
                dp[j]=dp1[j];
                dp1[j]=0;
            }
        }
        int inq=0;
        for(int i=1;i<=sum;i++)//记录无法称量的克数
            if(!dp[i])
            cun[inq++]=i;
        printf("%d\n",inq);
        for(int i=0;i<inq-1;i++)
            printf("%d ",cun[i]);
        if(inq)
            printf("%d\n",cun[inq-1]);

    }
    return 0;
}

01背包的思路:
这里与其他01背包不同的就是可以将克数减少为负数,要特别记录一下。
或者你也可以将所有的数值都右移sum(砝码和),这样负数就成了整数。

代码:

#include<stdio.h>
#include<string.h>
const int maxn=1e5+9;
int dp[maxn],dp_fu[maxn];//dp数组记录能够组成的正的克数,dp_fu数组记录能够组成的负砝码克数(相当于放到了称量物一边),
int a[maxn],b[maxn];//a数组记录给出的砝码,b数组记录无法称量的克数,
int main()
{
    int n;
    while(~scanf("%d",&n))
    {
        int sum=0,len=0;
        for(int i=0; i<n; i++)
            scanf("%d",&a[i]),sum+=a[i];
        memset(dp,0,sizeof(dp));
        memset(dp_fu,0,sizeof(dp_fu));
        dp[0]=1,dp_fu[0]=1;
        for(int i=0; i<n; i++)
        {
            for(int j=sum-a[i]; j>=0; j--)//先减去给出的砝码
                if(dp_fu[j]) dp_fu[j+a[i]]=1;
            for(int j=0; j<=sum; j++)//负数的区域,
                if(dp[j]&&j>=a[i]) dp[j-a[i]]=1;
                else if(dp[j]&&j<a[i]) dp_fu[a[i]-j]=1;

            for(int j=sum; j>=a[i]; j--)//再加上给出的砝码
                if(dp[j-a[i]]) dp[j]=1;
            for(int j=0; j<=sum; j++)//负数的区域
                if(j<=a[i]&&dp_fu[j]) dp[a[i]-j]=1;
                else if(dp_fu[j]&&j>a[i]) dp_fu[j-a[i]]=1;
        }
        int inq=0;
        for(int i=1; i<=sum; i++)//查找无法称量的克数
            if(!dp[i])
                b[inq++]=i;
        printf("%d\n",inq);
        for(int i=0; i<inq; i++)
            if(i==0) printf("%d",b[i]);
            else printf(" %d",b[i]);
        if(inq)printf("\n");
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值