hdu-3092（筛素数法+完全背包）

hdu-3092

 

（1）题意：

将一个数分为若干个数字X1……Xn，得到这若干个数字的最小公倍数LCM，

求LCM的最大值%m。

（2）思路：

让LCM最大，就是尽可能的让X1……Xn这些数字互质，所以先筛出素数（范围1~3000），然后将这些素数

视为物品，有一个背包的容积为n，就是完全背包了。

但是考虑n=7的情况，7 = 3+4，LCM（3,4） = 12.

4不是素数，但4时2的平方，所以，物品并不仅是所有的素数，还可能是素数的平方，所以还要遍历每个素数的k次方的值

（这里由于素数之间是互质的，所以素数的k次方之间也一定是互质的。）

 

（3）注意：

这里的结果对m取模，所以不能直接比较ans的大小，可以用log函数将乘法变为加法，

再设置一个dp数组，存储可以比较的值，ans数组记录结果。

 

（4）代码实现：

参考文章

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
const int maxn = 1e5+10;
int tot,ans[maxn],prim[maxn],vis[maxn];
double dp[maxn];
void Init(){
	tot = 0;memset(vis,0,sizeof(vis));
	for(int i=2;i<=3000;i++){
		if(!vis[i]) prim[tot++] = i;
		for(int j=0;j<tot&&i*prim[j]<=3000;j++){
			vis[i*prim[j]] = 1;
			if(i%prim[j]==0) break;
		}
	}
}
int main(void){
	Init();
	int n,m;
	while(~scanf("%d%d",&n,&m)){
		for(int i=0;i<=n;i++){
			dp[i] = 0;ans[i] = 1;
		}
		for(int i=0;i<tot&&prim[i]<=n;i++){
			for(int j=n;j>=prim[i];j--){
				double tp = log2(1.0*prim[i]);
				for(int p1=1,p2=prim[i];p2<=j;p1++,p2*=prim[i])
				if(dp[j-p2]+p1*tp>dp[j]){
					dp[j] = dp[j-p2]+p1*tp;
					ans[j] = ans[j-p2]*p2%m;
				}
			}
		}
		printf("%d\n",ans[n]);
	}
	
	return 0;
}