标签：ch int Luogu ll register li MAXN SCOI2009 P4161

题意

定义一个长度为 \(n\) 的置换的步数为将 \(P=(1,2,\cdots,n)\) 在该置换操作下变回原样的最小次数。

求有多少个 \(K\) 使得存在一个长度为 \(n\) 置换使得其步数为 \(K\)。

\(\texttt{Data Range:}1\leq n\leq 10^3\)

题解

类比一下这题我们可以得到如下转移方程：

\[f_{i,j}=f_{i,j-1}+\sum\limits_{p_j^k\leq i}f_{i-p_j^k,j-1} \]

这里不需要乘 \(p_j^k\) 的原因是我只要求出 \(K\) 有多少而不是 \(K\) 的和，然后就没了。

代码

#include<bits/stdc++.h>
using namespace std;
typedef int ll;
typedef long long int li;
const ll MAXN=1e4+51;
ll n,ptot;
li res=1;
ll prime[MAXN],np[MAXN]; 
li f[MAXN];
inline ll read()
{
    register ll num=0,neg=1;
    register char ch=getchar();
    while(!isdigit(ch)&&ch!='-')
    {
        ch=getchar();
    }
    if(ch=='-')
    {
        neg=-1;
        ch=getchar();
    }
    while(isdigit(ch))
    {
        num=(num<<3)+(num<<1)+(ch-'0');
        ch=getchar();
    }
    return num*neg;
}
int main()
{
    n=read(),f[0]=1;
    for(register int i=2;i<=n;i++)
    {
        if(!np[i])
        {
            prime[++ptot]=i;
        }
        for(register int j=1;i*prime[j]<=n;j++)
        {
            np[i*prime[j]]=1;
            if(!(i%prime[j]))
            {
                break;
            }
        }
    }
    for(register int i=1;i<=ptot;i++)
    {
        for(register int j=n;j>=1;j--)
        {
            for(register int k=prime[i];k<=j;k*=prime[i])
            {
                f[j]+=f[j-k];
            }
        }
    }
    for(register int i=1;i<=n;i++)
    {
        res+=f[i];
    }
    printf("%lld\n",res);
}

标签：ch,int,Luogu,ll,register,li,MAXN,SCOI2009,P4161
来源： https://www.cnblogs.com/Karry5307/p/13455443.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。