记录编号 |
589303 |
评测结果 |
AAAAAAAAAA |
题目名称 |
玩具 |
最终得分 |
100 |
用户昵称 |
小金 |
是否通过 |
通过 |
代码语言 |
C++ |
运行时间 |
0.522 s |
提交时间 |
2024-07-04 17:24:50 |
内存使用 |
2.02 MiB |
显示代码纯文本
#include<iostream>
#include<cstdio>
using namespace std;
int n,p;
long long inv[210];
long long dp[210][210],f[210][210],g[210][210];
long long ans;
long long qpow(long long x,long long y)
{
long long res=1;
while(y)
{
if(y&1) res=res*x%p;
x=x*x%p;
y>>=1;
}
return res;
}
int main()
{
freopen("toyy.in","r",stdin);
freopen("toyy.out","w",stdout);
scanf("%d%d",&n,&p);
dp[1][1]=1;
f[1][0]=1;
for(int i=1;i<=n;i++)
{
inv[i]=qpow(i,p-2);
}
for(int i=2;i<=n;i++)
{
for(int j=1;j<=i;j++)
{
dp[i][j]=(dp[i-1][j-1]*(j-1)%p+dp[i-1][j]*(i-j)%p)*inv[i]%p;
}
}
for(int i=0;i<=n;i++)
{
g[0][i]=1;
}
for(int i=1;i<=n;i++)
{
for(int j=0;j<=n;j++)
{
if(j) f[i][j]=g[i-1][j-1];
for(int k=1;k<=i;k++)
{
g[i][j]=(g[i][j]+f[k][j]*g[i-k][j]%p*dp[i][k]%p)%p;
}
}
}
for(int i=1;i<=n;i++)
{
ans=(ans+(f[n][i]-f[n][i-1]+p)%p*i%p)%p;
}
printf("%lld",ans);
return 0;
}