2020CCPC威海站 L-Clock Master （分组背包）

题意：给定你一个数b，然后让你找到一组数x{}使得x1 + x2 + x3 + … xn在不超过b的范围的情况下，使得这组xi的最大公倍数的值最大，然后输出这个最大值的倍数。

思路：保证这一组xi的LCM最大，那么就要保证这组xi之间应该是互素的，保证这些x是素数或者是素数的幂次倍。而且每种素数和它的各幂次之间只能从中取一个，因为这样才是最优的，那么我们就可以看到到这一步之后就是一个分组背包的问题了。
状态转移方程：
dp[ i ] = max(dp[ i ],dp[ i - k ] * k),k即代表着素数或者是这个素数的幂次倍。
然后去枚举每一组即可更新。

代码：

#include <bits/stdc++.h>

using namespace std;
typedef long long ll;
const double eps = 1e-6;
const int MAXN = 3e4 + 7;
int prime[MAXN],cnt,vis[MAXN];
double dp[MAXN],_log[MAXN];

void init(){
	for(int i = 2;i < MAXN;i ++){
		if(!vis[i]) prime[++cnt] = i;
		for(int j = 1;j <= cnt && prime[j] * i < MAXN;j ++){
			vis[prime[j]*i] = 1;
			if(i % prime[j] == 0)
				break;
		}
	}
}
void work(){
	for(int i = 0;i < MAXN;i ++) dp[i] = 1;//_log[i] = log(i);或者可以用一个log来更新这个dp值
	for(int i = 1;i <= cnt;i ++){
		for(int j = MAXN-1;j >= prime[i];j --){//分组背包 每组之间的元素互斥的 所以容量从大到小枚举 保证了互斥的正确性即小的不会去更新大的
			for(int k = prime[i];k <= j;k *= prime[i]){
					dp[j] = max(dp[j],dp[j-k] * k); //dp[j-k] + _log[k];
			}
		}
	}
}

int main(){
	init();
	work();
	int t;
	scanf("%d",&t);
	while(t--){
		int n;
		scanf("%d",&n);
		printf("%.9f\n",log(dp[n]));
	}
	return 0;
}