BZOJ 1025-游戏（数论+DP）

1025: [SCOI2009]游戏

Time Limit: 1 Sec   Memory Limit: 162 MB
Submit: 2355   Solved: 1541
[ Submit][ Status][ Discuss]

Description

　　windy学会了一种游戏。对于1到N这N个数字，都有唯一且不同的1到N的数字与之对应。最开始windy把数字按
顺序1，2，3，……，N写一排在纸上。然后再在这一排下面写上它们对应的数字。然后又在新的一排下面写上它们
对应的数字。如此反复，直到序列再次变为1，2，3，……，N。 
如： 1 2 3 4 5 6 对应的关系为 1->2 2->3 3->1 4->5 5->4 6->6 
windy的操作如下 
1 2 3 4 5 6 
2 3 1 5 4 6 
3 1 2 4 5 6 
1 2 3 5 4 6 
2 3 1 4 5 6 
3 1 2 5 4 6 
1 2 3 4 5 6 
这时，我们就有若干排1到N的排列，上例中有7排。现在windy想知道，对于所有可能的对应关系，有多少种可
能的排数。

Input

　　包含一个整数N,1 <= N <= 1000

Output

　　包含一个整数，可能的排数。

Sample Input

【输入样例一】
3
【输入样例二】
10

Sample Output

【输出样例一】
3
【输出样例二】
16

HINT

置换群的问题。对于每种对应关系，将其分解为循环节的形式，则变换次数为：各循环节长度的最小公倍数 
如 对应关系：1->2 2->3 3->1 4->5 5->4 6->6，转化后就是(1 2 3)(4 5)(6)，排数 = lcm(3,2,1)+1 = 7
循环节总长为N，所以问题转化为：
将N分解成若干正整数之和，求它们最小公倍数可能的取值个数 

#include<math.h>
#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
typedef long long ll;
#define  maxn 2005
ll dp[maxn][maxn],a[maxn]={1,1},b[maxn];
int  main(void)
{
	ll ans=0;
	ll n,i,j,k=0,p;
	scanf("%lld",&n);
	for(i=2;i<=n;i++)
	{
		if(a[i])
			continue;
		b[++k]=i;
		for(j=i*i;j<=n;j+=i)
			a[j]=1;
	}
	dp[0][0]=1;
	for(i=1;i<=k;i++)
	{
		for(j=0;j<=n;j++)
		{
			dp[i][j]=dp[i-1][j];
			for(p=b[i];j-p>=0;p*=b[i])
				dp[i][j]+=dp[i-1][j-p];
		}
	}
	for(i=0;i<=n;i++)
		ans+=dp[k][i];
	printf("%lld\n",ans);
	return 0;
}