2020 China Collegiate Programming Contest, Weihai Site.L Clock Master (数论，背包问题)

原题链接
题目大意:购买一个齿轮齿数为T的种需要T元，齿轮齿数为T的种可以表示T种时间状态。现给B元的预算，要求购买一些种，使得他们能表示的不同状态最多。
一开始，这些种都指向0。
一个3齿种与4齿种可以有12种状态。
但两个3齿种与一个4齿种也只有12种状态。
所以，想要组合数量更多，应该选择不同齿数的种。
如果选择的种的齿数为$a_1,a_2,a_3,...,a_n$，那么组合数就是最小公倍数$lcm(a_1,a_2,a_3,...,a_n)$
所以问题就转化成了，选择一些种，齿轮数不大于预算B，让lcm最大。
题目中要求，将结果以自然对数输出
也就是说每个中对答案的贡献是$log(齿轮数)$，为加法。
这样就变成了背包问题
我们只需要枚举质数的幂次，套用01背包即可。

#pragma GCC optmize(2)
#include <bits/stdc++.h>
#define int long long
#define endl '\n'
using namespace std;
const int N = 30010;
int primes[N], cnt;
bool st[N];

void get_primes(int n) {
	for (int i = 2; i <= n; i ++ ) {
		if (!st[i])
			primes[cnt ++ ] = i;
		for (int j = 0; primes[j] <= n / i; j ++ ) {
			st[primes[j] * i] = true;
			if (i % primes[j] == 0)
				break;
		}
	}
}
double lg[N];
double dp[N];

void init() {
	for (int i = 1; i < N; i++) {
		lg[i] = log(i);
	}
	for (int i = 0; i < cnt; i++) {//枚举各种质因子
		for (int j = N - 1; j >= primes[i]; j--) { //枚举容量
			for (int k = primes[i]; k < N; k *= primes[i]) {
				//每种质因子只能使用一种幂次
				if (j >= k) {
					dp[j] = max(dp[j], dp[j - k] + lg[k]);
				}
			}
		}
	}

}

void solve() {
	int b;
	scanf("%lld", &b);
	printf("%0.8lf\n", dp[b]);


}

signed main() {
	get_primes(N - 5);
	init();
	int t;
	cin >> t;
	while (t--)
		solve();
}