HDU 3092（完全背包）

省赛前疯狂补题中，博客潦草一点

HDU - 3092

这道题很难受，万万没想到是完全背包的变形

当是扩展思路记录一下，其实和一般的完全背包还是不一样的，他是把n个物品等效成一个物品，一次性加入背包中，这样做的好处是可以计算乘方，而如果一次性加入则会越界？感觉这个代码实现挺神奇的，当扩展思路了

题解链接

代码直接挂人家的，不想敲了。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
#define N 3500
#define LL long long
double dp[N];
LL ans[N];
LL prime[N],notprime[N];
LL gcd(LL a,LL b)
{
    return b?gcd(b,a%b):a;
}
void init() // 素数打表
{
    prime[0]=0;
    for(int i=2;i<=3000;i++)
    {
        if(!notprime[i])
        {
            prime[++prime[0]]=i;
            for(LL j=(LL)i*i;j<=3000;j+=i)
                notprime[j]=1;
        }
    }
}
int main()
{
    int n,m;
    init();
    while(scanf("%d %d",&n,&m)!=EOF)
    {
        for(int i=0;i<=n;i++)
        {
            dp[i]=0;
            ans[i]=1;
        }
        for(LL i=1;i<=prime[0]&&prime[i]<=n;i++)//遍历每一个物品
        {
            double tmp=log(prime[i]*1.0);//化乘法为除法
            for(LL j=n;j>=prime[i];j--) //实际是完全背包，为什么用01背包？
            {
                for(LL k=prime[i],q=1;k<=j;k*=prime[i],q++) //把k个物品等效成一个物品，这样写的好处是易于把乘法变成加法
                {
                    if(dp[j-k]+q*tmp>dp[j])
                    {
                        dp[j]=dp[j-k]+q*tmp;
                        ans[j]=ans[j-k]*k%m;
                    }
                }
            }
        }
        printf("%lld\n",ans[n]);
    }
    return 0;
}