2020ccpc 威海 Clock Master（分组背包）

题目大意：
$T$ 组数据，每组数据给出一个 $n$ ，构造若干个数的和相加不超过 $n$ ，即 $n\ge a_1+a_2+...+a_m$ ，求 $max(lcm(a_1,a_2,...,a_m))$ 对结果取以自然对数为底的对数
题目分析：
我们知道一个数$n$可通过质因数分解表示为$n=p_1^{cn{t}_1}\ast p_2^{cn{t}_2}\ast ...\ast p_x^{cn{t}_x}$
$m$ 个数的 $lcm$ 的本质是将这 $m$ 个数质因数分解后相同的质因子指数取最大值相乘，即 $$lcm=p_1^{max(cn{t}_1)}\ast p_2^{max(cn{t}_2)}\ast ...\ast p_y^{max(cn{t}_y)}$$
因此要使 $lcm$ 最大为了避免数的浪费就只需将 $n$ 拆成多个质数的幂次相加，即 $n=p_1^{cn{t}_1}+p_2^{cn{t}_2}+...+p_m^{cn{t}_m}$
此时这 $m$ 个数的 $lcm$ 即为 $p_1^{cn{t}_1}\ast p_2^{cn{t}_2}\ast ...\ast p_m^{cn{t}_m}$
此模型可以理解为每个质数的幂次作为一个组，体积为 $p^i$ 的数的价值是 $log(p^i)$ 然后就可以跑分组背包求解这个最优价值了
求解过程中由于答案很大所以对其取对数处理，利用 $ln(a\ast b)=lna+lnb$ 避免爆精度
要把3e4内的 $dp$ 数组直接递推好，询问时 $O(1)$ 查询，此时你会发现你还是 $tle$ ，此刻你离成功只差处理一个 $log$ 值的表，提前打好这个表就可以通过了
时间复杂度为 $O(\frac{n\ast n}{logn})$
具体细节见代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#define ll long long

#define inf 0x3f3f3f3f
using namespace std;
ll read()
{
	ll res = 0,flag = 1;
	char ch = getchar();
	while(ch<'0' || ch>'9')
	{
		if(ch == '-') flag = -1;
		ch = getchar();
	}
	while(ch>='0' && ch<='9')
	{
		res = (res<<3)+(res<<1)+(ch^48);//res*10+ch-'0';
		ch = getchar();
	}
	return res*flag;
}
const int maxn = 3e4+5;
const int mod = 1e9+7;
const double pi = acos(-1);
const double eps = 1e-8;
int pri[maxn],cnt;
bool vis[maxn];
void get_prime()//线筛 
{
	vis[1] = true;
	for(int i = 2;i < maxn;i++)
	{
		if(!vis[i]) pri[++cnt] = i;
		for(int j = 1;j <= cnt && i*pri[j] < maxn;j++)
		{
			vis[i*pri[j]] = true;
			if(i%pri[j] == 0) break;
		}
	}
}
double logg[maxn],dp[maxn];
void init()//预处理 
{
	for(int i = 1;i < maxn;i++)
		logg[i] = log(i);
	for(int i = 1;i <= cnt;i++)
			for(int v = maxn-1;v >= pri[i];v--)
				for(int j = pri[i];j <= v;j *= pri[i])
					dp[v] = max(dp[v],dp[v-j]+logg[j]);
}
int main()
{
	get_prime();
	init();
	int t = read();
	while(t--)
	{
		int n = read();
		printf("%.10f\n",dp[n]);
	}
	return 0;
}