Least common multiple HDU - 3092

http://acm.hdu.edu.cn/showproblem.php?pid=3092

看的题解 首先 拆出的数必须是素数的次幂 假如拆出了14这样的数 他在lcm中的贡献与拆出2和7是一样的 反而浪费了5个数

然后就是分组背包了 即x x^1 x^2 x^3算一组 只能选其中一个

最后 不能在DP过程中取模 这样无法比较大小 可以想到取个对数 化为浮点加法 同时另开一数组记录答案

 

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
const int maxn=3e3+10;

double dp[maxn];
int book[maxn],prime[maxn],ans[maxn];
int tot;

void init()
{
    int i,j;
    book[1]=1;
    for(i=2;i*i<=3000;i++){
        if(!book[i]){
            for(j=i+i;j<=3000;j+=i){
                book[j]=1;
            }
        }
    }
    for(i=1;i<=3000;i++){
        if(!book[i]) prime[++tot]=i;
    }
}

int main()
{
    double det;
    int s,m,i,j,k,l;
    init();
    while(scanf("%d%d",&s,&m)!=EOF){
        for(i=0;i<=s;i++){
            dp[i]=0.0,ans[i]=1;
        }
        for(i=1;i<=tot&&prime[i]<=s;i++){
            for(j=s;j>=prime[i];j--){
                for(k=prime[i],l=1;k<=j;k*=prime[i],l++){
                    det=(double)(l)*log2((double)(prime[i]));
                    if(dp[j]<dp[j-k]+det){
                        dp[j]=dp[j-k]+det,ans[j]=(ans[j-k]*k)%m;
                    }
                }
            }
        }
        printf("%d\n",ans[s]);
    }
    return 0;
}