hdu3092(数论+完全背包)

题意：将一个数拆分成几个数的和，并求出这几个数最小公倍数模M后的最大值

代码

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;
typedef long long ll;
const int maxn=3e3+10;
void io() { ios::sync_with_stdio(0);cin.tie(0);cout.tie(0); }
double dp[maxn];
int ans[maxn],prime[maxn],cnt;
int not_p[maxn];
int s,m;
void init(){
    not_p[1]=1;
    for(int i=2;i<=3e3;i++){
        if(!not_p[i]){
            prime[++cnt]=i;
            for(int j=i+i;j<=3e3;j+=i)
                not_p[j]=1;
        }
    }
}
int main() {
    init();
    while(~scanf("%d%d",&s,&m)){
        for(int i=0;i<=s;i++){
            dp[i]=0.0,ans[i]=1;
        }
        for(int i=1;i<=cnt&&prime[i]<=s;i++){
            for(int j=s;j>=prime[i];j--){
                for(int k=prime[i],l=1;k<=j;k*=prime[i],l++){
                    double x=1.0*l*log2((double)(prime[i]));
                    if(dp[j]<dp[j-k]+x) dp[j]=dp[j-k]+x,ans[j]=(ans[j-k]*k)%m;
                }
            }
        }
        cout<<ans[s]<<endl;
    }

    return 0;
}