SCOI2009[游戏] dp

bzoj1025:http://www.lydsy.com/JudgeOnline/problem.php?id=1025

看了题，手推了前4组，发现这道题貌似就是求所有构成环的最小公倍数有多少个。

一看就乐了，这道题拐个弯就简单了，于是想要dfs+一个数组储存lcm的出现标记，然后扫一遍就好了，然而理想很丰满，现实很残酷，按照类似string递推式估算，发现总次数爆了int，顿时就蒙了。

怎么做捏，百思不得其解了半个小时，搜了一下题解，发现都是dp做的，熬~~

对于任何一个lcm来说，都可以用唯一分解定理分解开来，按照题意如果这些因子的一些幂的值之和<=n(不足可以补1)，就找到了一个lcm.          不互素的两个环长的lcm如果和互素的两个环长的lcm相等，前者的和>后者，就不能组成更多种的lcm,所以互素的才是最优的

dp[i][j]表示前i个素数组成的和为j的lcm的个数dp[i][j]=∑dp[i-1][j-xi^k]    因为互素，所以前面的*一个质数方案数不变，把所有能转移的加起来(还要处理这个xi不选)

然后ans=∑dp[cnt][i]  i<=n

AC代码：

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<queue>
#include<cmath>
#define LL long long
using namespace std;
/*
其实是求组成环的最小公倍数有哪些，+1就是排数 
转化一下就是尝试把N成k个不小于1的整数，求它的最小公倍数的个数 
*/
int n;
LL f[1005][1005];//前i个因子搞出来和为j的方案数
//f[i][j]=∑f[i-1][j-x]  x为第i个因子的幂 
bool flag[1005];
int cnt;
int prime[1005];
void init()
{
	for(int i=2;i<=n;i++)
	{
		if(!flag[i])
		{
			prime[++cnt]=i;		
			//cout<<prime[cnt]<<",";
		}
		for(int j=1;j<=cnt&&i*prime[j]<=n;j++)
		{
			flag[i*prime[j]]=1;
			if(i%prime[j]==0)break;
		}
	}
}
int main()
{
	cnt=0;
	scanf("%d",&n);
	init();
	f[0][0]=1;
	for(int i=1;i<=cnt;i++)
	{
		for(int j=0;j<=n;j++)f[i][j]+=f[i-1][j];
		for(int j=prime[i];j<=n;j*=prime[i])
		{
			for(int k=0;k<=n-j;k++)
			{
				f[i][j+k]+=f[i-1][k];
			}
		}
	}
	LL ans=0;
	for(int i=0;i<=n;i++)ans+=f[cnt][i];
	printf("%lld\n",ans);
	return 0;
}