hdu 3092 Least common multiple

题目链接：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;
}