[唯一分解+分组dp] 2020ccpc威海 L

题目

一个钟表有n个指针，每一个指针有t个刻度每个刻度。所有制真的刻度不能超过d,求所有指针合起来的不同方向最多有多少种。

思路

每个指针的不同方向种类即他们的$lcm$
题意转换为：
求$MAX(lcm(t_1,t_2,t_3,t_4,...,t_n))$
$t_1+t_2+t_3+t_4+...+t_n<=b$
由$lcm(a,b)=a\ast b/gcd(a,b)$
要让$lcm$的值最大，那就让$a,b$互质，$gcd(a,b)=1$
由唯一分解定理，一个数可由有限个质数组成，那么我们先于处理$30000$内的所有质数，用$p[i][k]$表示第$i$个质数的$k$次方，用$a[i][k]$来储存当前$ln(p[i][k])$，最后分组背包即可。

代码

#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<cctype>
#include<ctime>
#include<iostream>
#include<string>
#include<map>
#include<queue>
#include<stack>
#include<set>
#include<vector>
#include<iomanip>
#include<list>
#include<bitset>
#include<sstream>
#include<fstream>
#include<complex>
#include<algorithm>
#if __cplusplus >= 201103L
#include <unordered_map>
#include <unordered_set>
#endif
#define int long long
using namespace std;
const int INF = 0x3f3f3f3f;
int prime[30100],b[3510],p[3510][210];
double dp[30010],a[3500][210];
bool vis[30010];
int tot=0;
void getprime(int n){
	vis[1]=1;
	vis[0]=1;
	for(int i=2;i<=n;i++){
		if(!vis[i]){
			prime[++tot]=i;
		}
		for(int j=1;j<=tot&&prime[j]*i<=n;j++){
			vis[prime[j]*i]=1;
			if(i%prime[j]==0){
				break;
			}
		}
	}
}
signed main(){
//	ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
	int t;
	cin>>t;
	getprime(30010);
	for(int i=1;i<=tot;i++){
		int ans=1;
		for(int j=1;ans<=30005;j++){
			ans*=prime[i];
			p[i][j]=ans;
			a[i][j]=log(ans);
			b[i]=j;
		}
	}
	for(int i=1;i<=tot;i++){
		for(int j=30000;j>0;j--){
			for(int k=1;k<=b[i];k++){
				if(j-p[i][k]>=0) dp[j]=max(dp[j],dp[j-p[i][k]]+a[i][k]);
			}
			dp[j]=max(dp[j],dp[j-1]);
		}
	}
	while(t--){
		int n;
		cin>>n;
		printf("%.9lf\n",dp[n]);
	}
    return 0;
}