HDU - 3092 Least common multiple 解题报告（数论+完全背包）

题目描述

题目：http://acm.hdu.edu.cn/showproblem.php?pid=3092
大概意思：给定一个数s和模m，将s分解为几个数相加的形式，求这几个加数的最小公倍数最大值%m的值。

思路分析

这道题我一开始是一道数论题，就想着：既然要最小公倍数的最大值，最大公倍数lcm=a*b/k(k=gcd(a,b))，那么，我们有限考虑的应该是k最小，也就是1的情况，那么，我们就想到分解出来的这些数应该是质数。因此，我就按照素数分解的方法去求，结果，wa。。。
这里我忽略了一点：当分解为这些素数后，仍可能剩下一个数为合数，可能导致这个分解的结果不是最优解。而我们知道，合数还能再写成素数的和，这样，才能保证得到的积是最大值。那么，这个时候，求积的时候，一个素数就可能就几次方。那么，这个问题就变成了完全背包问题：每一个素数就是物品，每一个素数可以取无限多个，而背包容量就是被分解的数，要求的就是在一定背包容量下，背包内的数的积的最大值。
另外这道题还有个处理：因为积可能会很大，所以，我们可以取对数，将数化小进行比较。

完整代码

#include <iostream>
#include <cmath>
using namespace std;
const int maxn=3005;
int primes[maxn]={0};
int cnt=0;
void make_list()//素数筛打表
{

	bool tmp[maxn]={0};
	for(int i=2;i<maxn;i++)
	{
		if(!tmp[i])
		{
			primes[cnt++]=i;
			for(int j=i*i;j<maxn;j+=i)
			{
				tmp[j]=1;
			}
		}
	}
	return ;
}

int main()
{
	make_list();
	int s,mod;
	while(cin>>s>>mod)
	{
		double dp[maxn]={0};
		int ans[maxn]={0};
		for(int i=0;i<=s;i++) ans[i]=1;
		for(int i=0;i<cnt;i++)
		{
			double tmp=log10(primes[i]*1.0);//取对数
			for(int j=s;j>=primes[i];j--)
			{
				for(int counter=1,k=primes[i];j>=k;k=k*primes[i],counter++)
				{
					if(dp[j]<dp[j-k]+counter*tmp)//比较
					{
						dp[j]=dp[j-k]+counter*tmp;
						ans[j]=(ans[j-k]*k)%mod;
					}
				} 
			}
		}
		cout<<ans[s]<<endl;
	}
	
	return 0;
}