bzoj1025 [SCOI2009]游戏

http://www.elijahqi.win/archives/2912
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
考虑我的一个答案如果能够是符合题意的 那么显然我应该是我这个排列里所有置换的长度的lcm+1 那么问题转化有k个正整数他们的和为n求最小公倍数的数目 因为k不知道所以转化为验证一个答案是否是可行答案 于是有如下结论当 ans=p1^a1+p2^a2+…pn^an<=n的时候则该ans有解
证明详情：https://blog.csdn.net/wyfcyx_forever/article/details/40211739
于是分组背包做一下即可设dp[i][j]表示用前i个质数凑出和为j的方案数 这个j我从0开始只是每次带着走而已..
这里还有一篇看起来更为严谨的证明https://www.cnblogs.com/phile/p/4473192.html

#include<cstdio>
#include<cctype>
#include<algorithm>
#define N 1010
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S) {T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(!isdigit(ch)) {if (ch=='-') f=-1;ch=gc();}
    while(isdigit(ch)) x=x*10+ch-'0',ch=gc();
    return x*f;
}
long long dp[N][N];int prime[N],n,tot;bool not_prime[N]; 
int main(){
    freopen("bzoj1025.in","r",stdin);
    n=read();
    for (int i=2;i<=n;++i){
        if (!not_prime[i]) prime[++tot]=i;
        for (int j=1;prime[j]*i<=n;++j){
            not_prime[prime[j]*i]=1;
            if (i%prime[j]==0) break;
        }
    }long long ans=0;dp[0][0]=1;
    for (int i=1;i<=tot;++i)
        for (int j=0;j<=n;++j){
            dp[i][j]=dp[i-1][j];
            for (int k=prime[i];k<=j;k*=prime[i])
                dp[i][j]+=dp[i-1][j-k];
        }
    for (int i=0;i<=n;++i) ans+=dp[tot][i];printf("%lld\n",ans);
    return 0;
}