(1)题意:
将一个数分为若干个数字X1……Xn,得到这若干个数字的最小公倍数LCM,
求LCM的最大值%m。
(2)思路:
让LCM最大,就是尽可能的让X1……Xn这些数字互质,所以先筛出素数(范围1~3000),然后将这些素数
视为物品,有一个背包的容积为n,就是完全背包了。
但是考虑n=7的情况,7 = 3+4,LCM(3,4) = 12.
4不是素数,但4时2的平方,所以,物品并不仅是所有的素数,还可能是素数的平方,所以还要遍历每个素数的k次方的值
(这里由于素数之间是互质的,所以素数的k次方之间也一定是互质的。)
(3)注意:
这里的结果对m取模,所以不能直接比较ans的大小,可以用log函数将乘法变为加法,
再设置一个dp数组,存储可以比较的值,ans数组记录结果。
(4)代码实现:
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
const int maxn = 1e5+10;
int tot,ans[maxn],prim[maxn],vis[maxn];
double dp[maxn];
void Init(){
tot = 0;memset(vis,0,sizeof(vis));
for(int i=2;i<=3000;i++){
if(!vis[i]) prim[tot++] = i;
for(int j=0;j<tot&&i*prim[j]<=3000;j++){
vis[i*prim[j]] = 1;
if(i%prim[j]==0) break;
}
}
}
int main(void){
Init();
int n,m;
while(~scanf("%d%d",&n,&m)){
for(int i=0;i<=n;i++){
dp[i] = 0;ans[i] = 1;
}
for(int i=0;i<tot&&prim[i]<=n;i++){
for(int j=n;j>=prim[i];j--){
double tp = log2(1.0*prim[i]);
for(int p1=1,p2=prim[i];p2<=j;p1++,p2*=prim[i])
if(dp[j-p2]+p1*tp>dp[j]){
dp[j] = dp[j-p2]+p1*tp;
ans[j] = ans[j-p2]*p2%m;
}
}
}
printf("%d\n",ans[n]);
}
return 0;
}