2020 China Collegiate Programming Contest, Weihai Site L. Clock Master（分组背包+预处理）

传送门

题意

给你一个正整数$n$，你需要将$n$拆成若干个正整数的和。设拆成了$t_1,t_2,...,t_s$这$s$个正整数，对于任意自然数$k$，存在向量$(k$ $mod$ $t_1,$$k$ $mod$ $t_2,...,$ $k$ $mod$ $t_s)$，这些向量中的不同向量个数记为$x$。
因此每种拆分对应一个不同的向量个数，现在让你求一个使$x$最大的拆分，输出$\ln x$。

数据范围：$1⩽n⩽3\times 10^4$

思路

① 一个结论： $x=lcm(t_1,t_2,...,t_s)$
由于本文主要侧重分组背包部分，对这个结论感兴趣的读者可以自己到网上寻找证明。
② 接下来就是考虑怎么样拆分。考虑最优解一定存在拆成的数两两互质的形式，多余的部分都可以拆成若干个$1$来凑数($1$与任何数互质)，这样计算答案就变成了计算$\ln t_1+\ln t_2+...+\ln t_s$。
证明： 假设最优解为$t_1,t_2,...,t_s$
1、若这些数两两互质，则成立。
2、若$gcd(t_x,t_y)=d>=2$，则将$t_x$从解中移除，将 $t_x/d$ 和 $t_x-t_x/d$ 个$1$加入到解中，原解的值不变。
3、重复步骤2直到所有数两两互质。
③考虑拆成若干个质数的幂次形式，即每个$t$都是$p^i$的形式，其中$p$为质数，$i$为自然数。
一个$p^i$的价值是$\ln p^i$，这不就是背包吗。
④对于每个$p$，只能选一次$p^i$，但不同的质数是可以任意选的，因此做一个分组背包，先筛一下质数，然后按$p$分组。预处理答案然后$O(1)$回答。
⑤总时间复杂度为$O(n$${\sum }_{i=2,i为质数}^n{\log }_{i}n$$)$ ，$\ln t$需要预处理不然会T。

#include<bits/stdc++.h>
#define ll long long
#define rep(i,a,b) for(int i=(a);i<=(b);i++)
#define dep(i,a,b) for(int i=(a);i>=(b);i--)
using namespace std;
const int maxn=4e4+5;
const double ln = log(exp(1));
template <typename T>
inline void read(T &X){
    X=0;int w=0; char ch=0;
    while(!isdigit(ch)) {w|=ch=='-';ch=getchar();}
    while(isdigit(ch)) X=(X<<3)+(X<<1)+(ch^48),ch=getchar();
    if(w) X=-X;
}
int n,m,k,cnt;
int a[maxn],c[maxn];
double dp[maxn];
double ans[maxn],lg[maxn];
int flag;
void init(){
    rep(i,2,30000){
        if(!c[i]){
            c[++cnt] = i;
            for (int j = i * 2; j <= 30000;j+=i)
                c[j] = 1;
        }
        lg[i] = log(i) / ln;
    }
}
void solve(){
    n = 30000;
    rep(i,1,cnt){
        dep(j,n,c[i]){
            ll tp = c[i];
            while(tp<=j){
                dp[j] = max(dp[j], dp[j - tp] + lg[tp]);
                tp *= c[i];
            }
        }
    }
    rep(i, 2, n) ans[i] = max(ans[i - 1], dp[i]);
}
int main(){
    init();
    solve();
    int T=1,cas=1;
    read(T);
    while (T--){
        read(n);
        cout << fixed << setprecision(12) << ans[n] << endl;
    }
    return 0;
}