[BZOJ1025]SCOI2009游戏|置换的循环|DP

一开始傻逼看错题。。

不难发现这是一个置换的循环问题，那么一个给定置换的行数就是这个置换各循环节长度的lcm，所以答案就是求n的加法分解的lcm数。。但是暴力搞肯定要T。。考虑把每个循环节长度质因数分解，如果有两个长度不互质，可以让其中一个除以gcd，那么这几个长度的和变小，且lcm不变。。由此可知我们只要统计所有长度互质的情况即可，再进一步可以直接统计所有长度都只有一个质因子的情况。。f[i][j]表示用前i个素数，目前长度和为j的种数，YY一下转移发现这是一个背包DP。。然后以为1对lcm是没有贡献的，所有当长度不足n时用1补足即可，所以答案就是Σf[cnt][i](1<=i<=n,cnt为n以内的素数个数)。。贴个链接http://www.cnblogs.com/iwtwiioi/p/4104385.html

#include<iostream>
#include<cstdio>
#include<memory.h>
#define ll long long
#define N 1005
using namespace std;
ll f[N][N],ans=0;
int n,cnt=0,p[N],u[N];
void getp()
{
	int i,j;
	memset(u,0,sizeof(u));
	for (i=2;i<=n;i++)
	{
		if (!u[i]) p[++cnt]=i;
		for (j=1;i*p[j]<=n&&j<=cnt;j++)
		{
			u[i*p[j]]=1;
			if (i%p[j]==0) break;
		}
	}
}
void solve()
{
	int i,j,t;
	f[0][0]=1;
	for (i=1;i<=cnt;i++)
	{
		for (j=0;j<=n;j++)
		{
			f[i][j]=f[i-1][j];
			for (t=p[i];t<=j;t*=p[i]) f[i][j]+=f[i-1][j-t];
		}
	}
	for (i=0;i<=n;i++) ans+=f[cnt][i];
}
int main()
{
	scanf("%d",&n);
	getp();
	solve();
	printf("%I64d",ans);
}