题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=3092
题意:给了一个数把它分成若干个数的和的形式,求这若干个数的最小公倍数的最大值是多少
分析:gcd(a,b) = 1时lcm最大,所以目标是尽可能将这个数处理成全为素数的形式。先用线性素数筛出3000里的所有质数,然后枚举每个质数,尽可能将这个质数填到n中,此时就类似于完全背包,每一个大于等于n的数都相当于一个背包,状态转移式子为dp[j] = max(dp[j - k] * k,dp[j])。填质数的时候因为要对m取模,会对dp的结果产生影响,因此对两边取对数,计算对数的值(取余数不会产生影响,3^10>3000),再开个数组保存真正的答案即可。
代码:
#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>
using namespace std;
int prime[1010],vis[3010];
int cnt;
double dp[3010];
int ans[3010];
void init()
{
for(int i = 2; i < 3010; i++)
{
if(!vis[i])
{
prime[cnt++] = i;
for(int j = i; j < 3010; j += i)
vis[j] = 1;
}
}
}
int main()
{
int s,m;
init();
//printf("%d\n",cnt);
while(~scanf("%d%d",&s,&m))
{
memset(dp,0,sizeof(dp));
for(int i = 0; i <= s; i++)
ans[i] = 1;
for(int i = 0; i < cnt && prime[i] < s; i++)
{
double t = log10(prime[i] * 1.0);
for(int j = s; j >= prime[i]; j--)
{
for(int k = prime[i],num = 1; k <= j; k *= prime[i],num++)
{
if(dp[j - k] + num * t > dp[j])
{
dp[j] = dp[j - k] + num * t;
ans[j] = (ans[j - k] * k) % m;
}
}
}
}
cout<<ans[s]<<endl;
}
return 0;
}