比赛 |
2022级DP专题练习赛4 |
评测结果 |
AAAAAAAAAA |
题目名称 |
排列计数 |
最终得分 |
100 |
用户昵称 |
op_组撒头屯 |
运行时间 |
0.249 s |
代码语言 |
C++ |
内存使用 |
12.02 MiB |
提交时间 |
2023-02-20 19:10:32 |
显示代码纯文本
#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define ull unsigned long long
#define ld long double
#define pii pair<int,int>
#define vi vector<int>
#define si set<int>
#define unsi unordered_set<int>
#define qi queue<int>
#define sti stack<int>
#define pqi priority_queue<int>
#define mii map<int,int>
#define unmii unordered_map<int,int>
#define fi first
#define se second
#define pb push_back
#define clr(f,n) memset(f,0,sizeof(int)*(n))
#define cpy(f,g,n) memcpy(f,g,sizeof(int)*(n))
const int N=1000000+5;
int head[N],ne=0;
ll n,mod;
ll f[N],sz[N];
ll fac[N],inv[N];
ll fst(ll x,ll y){
ll ans=1;
while(y){
if (y&1)ans=ans*x%mod;
x=x*x%mod;y>>=1;
}
return ans;
}
void init(){
fac[1]=1;int lim=min(n,mod-1);
for (int i=2;i<=lim;i++)fac[i]=fac[i-1]*i%mod;
inv[lim]=fst(fac[lim],mod-2);
for (int i=lim;i>=1;i--)inv[i-1]=inv[i]*i%mod;
}
ll lucas(ll x,ll y){
if (x<y)return 0;
if (!x)return 1;
ll ans=0,xx=x%mod,yy=y%mod;
if (xx>=yy) ans=fac[xx]*inv[yy]%mod*inv[xx-yy]%mod;
return lucas(x/mod,y/mod)*ans%mod;
}
int main(){
freopen ("permzj.in","r",stdin);
freopen ("permzj.out","w",stdout);
scanf("%lld%lld",&n,&mod);init();
f[0]=1;
for (int i=n;i>=1;i--){
int lc=i*2,rc=i*2+1;
if (lc>n)lc=0;
if (rc>n)rc=0;
sz[i]+=sz[lc]+sz[rc]+1;
f[i]=lucas(sz[i]-1,sz[lc])*f[lc]%mod*f[rc]%mod;
}
printf("%lld\n",f[1]);
return 0;
}