[SCOI2009]游戏

题目：
https://ac.nowcoder.com/acm/problem/20271

给定$n$，在满足${\sum }_{i=1}^k{a}_i=n,a_i\ge 1$的情况下。求有多少种$lcm(a_1,a_2...a_k)$

思路：
假设$lcm(a_1,a_2...a_k)=m={\prod }_{i=1}^x{p}_i^{b_i}$，因为$p_i\le n$，所以$x<n$。可以构造令$a_1=p_1^{b_1},a_2=p_2^{b_2}...a_x=p_x^{b_x}$，剩下的用$1$补就能满足${\sum }_{i=1}^k{a}_i=n$（最优的情况）。此时${\sum }_{i=1}^x{p}_i^{b_i}\le n$
只要再证明${\sum }_{i=1}^x{p}_i^{b_i}>n$无解。

证：
理想情况$a_1=p_1^{b_1}{k}_1,a_2=p_2^{b_2}{k}_2...a_x=p_x^{b_x}{k}_x(gcd(p_i,k_i)=1)$，此时${\sum }_{i=1}^x{a}_i>n$
更坏的情况比如$a_1=p_1^{b_1}{k}_1,a_2=p_2^{b_2}{k}_2...a_{x-1}=p_{x-1}^{b_{x-1}}{p}_x^{b_x}{k}_{x-1}$
因为$a_1\ge p_1^{b_1}...a_{x-2}\ge p_{x-2}^{b_{x-2}},a_{x-1}\ge p_{x-1}^{b_{x-1}}+p_x^{b_x}(p_{x-1}^{b_{x-1}}{p}_x^{b_x}\ge p_{x-1}^{b_{x-1}}+p_x^{b_x})$，所以${\sum }_{i=1}^{x-1}{a}_i>n$。所以无解。

综上只需要判断有多少种$b$数组满足${\sum }_{i=1}^x{p}_i^{b_i}\le n$，$p_x$为小于等于$n$的最大质数。可以用背包$dp$。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1000+5;
int pri[N],tot;
ll f[N];
bool vis[N];
int n;
void sieve(){
	for(int i=2;i<=n;i++){
		if(!vis[i]){
			pri[++tot]=i;
		}
		for(int j=1;j<=tot;j++){
			if(pri[j]*i>n) break;
			vis[pri[j]*i]=1;
			if(i%pri[j]==0) break;
		}
	}
}
int main()
{
	scanf("%d",&n);
	sieve();
	f[0]=1;
	for(int i=1;i<=tot;i++)
	 for(int j=n;j>=0;j--){
	 	for(int k=pri[i];k<=j;k*=pri[i]){
	 		f[j]+=f[j-k];
		 }
	 }
	ll ans=0;
	for(int i=0;i<=n;i++)ans+=f[i];
	printf("%lld",ans);
	return 0;
}