bzoj1025 [SCOI2009]游戏（数论+背包dp）

考虑每一个置换的排数其实就是所有循环节长度的lcm+1。

考虑一个数x是否可能是排数，即是否可能表示为若干数的lcm，且和<=n。

设$x=p_1^{q_1}\cdots p_k^{q_k}$，则最小的和为$p_1^{q_1}+p_2^{q_2}+\cdots +p_k^{q_k}$，如果<=n则合法，否则不合法。

因此我们可以枚举每一个质因子，当做分组背包，f[i][j],表示考虑前i个质数,在最小规则下，和为j的排数有多少种。每组的物品为$p^1,p^2,\cdots,p^k$。
复杂度$O(n/lgn*logn*n)$

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
#define inf 0x3f3f3f3f
#define ll long long
#define N 1010
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,tot=0,prime[200];ll f[200][N],ans=0;//f[i][j],考虑前i个质数,在最小规则下，和为j的排数有多少种
bool notprime[N];
inline void getprime(){
    notprime[1]=1;
    for(int i=2;i<=n;++i){
        if(!notprime[i]) prime[++tot]=i;
        for(int j=1;prime[j]*i<=n;++j){
            notprime[prime[j]*i]=1;
            if(i%prime[j]==0) break;
        }
    }
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();getprime();f[0][0]=1;
    for(int i=1;i<=tot;++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=j;k<=n;++k)
                f[i][k]+=f[i-1][k-j];
    }
    for(int i=0;i<=n;++i) ans+=f[tot][i];
    printf("%lld\n",ans);
    return 0;
}