bzoj1025[SCOI2009]游戏 DP 置换 筛法

题意：给你一个长度为n的排列，指定每一个数的对应关系，问能有多少种排法。

很明显不能直接做，转换成图论模型很显然，直接把对应关系看成边，然后求的就是所有环长度的lcm总数，然后。。。然后我就不是很懂怎么转化了= =
其实这个东西就等价于求一个和为n的数列的lcm种数。。因为环的总长度一定为n。。(其实也可以转化为置换的模型)
那么对于一个排数x，验证他的方法就是对x进行质因数分解，即x=$p1^{a1}*p2^{a2}.....pk^{ak}$.
那么有一个结论：假如$p1^{a1}+p2^{a2}+...+pk^{ak}<=n$，那么x就是可行的答案，具体证明省略，需要看请移步
然后背包搞搞就可以了。

#include<cstdio>
#include<cstring>
#include<algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
const int N=1e5+5;
int n,m;
int p[N],vis[N],tot;
long long f[1005][1005],ans;
int main()
{
    scanf("%d",&n);
    fo(i,2,n)
    {
        if (!vis[i])p[++tot]=i;
        fo(j,1,tot)
        {
            if (p[j]*i>1000)break;
            vis[p[j]*i]=1;
            if (i%p[j]==0)break;
        }   
    }
    f[0][0]=1;
    fo(i,1,tot)
    {
        fo(j,0,n)
        f[i][j]+=f[i-1][j];
        for(int j=p[i];j<=n;j*=p[i])
        fo(k,j,n)
        f[i][k]+=f[i-1][k-j];
    }   
    fo(i,0,n)ans+=f[tot][i];
    printf("%lld\n",ans);
}