[ZJOI2010]排列计数
题目描述
称一个1,2,...,N的排列P1,P2...,Pn是Magic的,当且仅当2<=i<=N时,Pi>Pi/2. 计算1,2,...N的排列中有多少是Magic的,答案可能很大,只能输出模P以后的值
输入输出格式
输入格式:
输入文件的第一行包含两个整数 n和p,含义如上所述。
输出格式:
输出文件中仅包含一个整数,表示计算1,2,⋯, ���的排列中, Magic排列的个数模 p的值。
输入输出样例
输入样例#1:
20 23
输出样例#1:
16
说明
100%的数据中,1 ≤N ≤ 10^6, P≤ 10^9,p是一个质数。
题解
题目就是让你求出有多少种排列满足小根堆
由于每个点的限制条件构成了一颗二叉树
所以我们可以设f[i]表示以值i为根的小根堆的数目
那么显然ta的子树内的所有取值都要比他大
所以我们可以统计出它及它的子树有多少大Num[i]
这样\(f[i] =C(Num[i],Num[ls])*f[ls]*f[rs]\)就可以转移了
因为到i时有Num[i]个值待选
那么分配Num[ls]个值给左子树的方案数乘上f[ls]和f[rs]就是f[i]了
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
const int M = 1000005 ;
using namespace std ;
int n , mod ;
int f[M] , fac[M] , Num[M] ;
inline int fpw(int Base , int k) {
int temp = 1 ;
while(k) {
if(k & 1) temp = (1LL * temp * Base) % mod ;
Base = (1LL * Base * Base) % mod ; k >>= 1 ;
}
return temp ;
}
inline int C(int n , int m) {
return (1LL * fac[n] % mod * fpw((1LL * fac[m] * fac[n - m]) % mod , mod - 2) % mod) % mod ;
}
int main() {
cin >> n >> mod ;
fac[0] = 1LL ;
for(int i = 1 ; i <= n ; i ++) fac[i] = (1LL * fac[i - 1] * i) % mod ;
for(int i = n ; i >= 1 ; i --) {
Num[i] = 1 ;
if(i * 2 <= n) Num[i] += Num[i << 1] ;
if(i * 2 + 1 <= n) Num[i] += Num[i << 1 | 1] ;
if(i * 2 + 1 <= n) f[i] = (1LL * C(Num[i] - 1 , Num[i << 1]) * f[i << 1] % mod * f[i << 1 | 1]) % mod ;
else if(i * 2 <= n) f[i] = f[i << 1] ;
else f[i] = 1 ;
}
cout << f[1] << endl ;
return 0 ;
}