省赛前疯狂补题中,博客潦草一点
这道题很难受,万万没想到是完全背包的变形
当是扩展思路记录一下,其实和一般的完全背包还是不一样的,他是把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;
}