[2020 China Collegiate Programming Contest Weihai Site - 重现赛] L Clock Master

分析：

这个题是一个分组背包问题。$v=a{1}^{p1}+a{2}^{p2}+...+a{n}^{pn}$
每一个质因数的幂次的集合就是一个分组。
另外，求log要先打表，不然会t

代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;//三年竞赛一场空，不开long long见祖宗 
//typedef __int128 lll;
#define print(i) cout << "debug: " << i << endl
#define close() ios::sync_with_stdio(0), cin.tie(0), cout.tie(0)
#define mem(a, b) memset(a, b, sizeof(a))
#define pb(a) push_back(a)
#define x first
#define y second
typedef pair<int, int> par;
const ll mod = 1e9 + 7;
const int maxn = 3e4 + 10;
const int inf = 0x3f3f3f3f;
int prime[maxn], cnt;
int vis[maxn];
double dp[maxn];
double Log[maxn];

void getprime()
{
    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 init()
{
    for(int i = 1; i < maxn; i++)
        Log[i] = log(i);
    for(int i = 1; i <= cnt; i++)
        for(int v = maxn - 1; v >= 1; v--)
            for(int j = prime[i]; j <= v; j *= prime[i])
                dp[v] = max(dp[v], dp[v - j] + Log[j]);
}

int main()
{
    getprime();
    init();
    int x;
    int n; cin >> n;
    for(int i = 1; i <= n; i++)
        cin >> x, printf("%.7lf\n", dp[x]);
}