P2606 [ZJOI2010]排列计数

p[i]>p[i/2]可以看出是小根堆的性质，即儿子大于父亲，所以问题转化为了把1—n插入n个节点的小根堆，有多少种插入方法；

对于一个节点u,设siz[u]为包括u节点的子树的大小，为了维护小根堆的性质，u必为最小值,其左儿子有C(siz[u]-1,siz[u2])种取法，结合乘法原理，dp[u]=C(siz[u]-1,siz[u2])dp[siz[u2]]dp[siz[u2+1]];
求C时应求乘法逆元。
代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#define LL long long
#define MAXN 1000010
using namespace std;
LL n,p,siz[MAXN],dp[MAXN],jc[MAXN];
LL quick_pow(LL a,LL b,LL Mod){
	LL ans=1,sum=a;
	while(b){
		if(b&1) ans=(ans*sum)%Mod;
		sum=(sum*sum)%Mod;
		b>>=1;
	}
	return ans;
}
LL dfs(int u){
	if(siz[u]<=2) return 1;
	if(siz[u]==3) return 2;
	if(!dp[siz[u*2]]) dp[siz[u*2]]=dfs(u*2);
	if(!dp[siz[u*2+1]]) dp[siz[u*2+1]]=dfs(u*2+1);
	return jc[siz[u]-1]*quick_pow(jc[siz[u*2]],p-2,p)%p*quick_pow(jc[siz[u*2+1]],p-2,p)%p*dp[siz[u*2]]%p*dp[siz[u*2+1]]%p;
}
int main(){
	scanf("%lld%lld",&n,&p);
	for(int i=1;i<=n;i++){
		int t=i;
		while(t) siz[t]++,t>>=1;
	} 
	jc[0]=1;
	for(int i=1;i<=n;i++) jc[i]=(i*jc[i-1])%p;
	printf("%lld\n",dfs(1));
	return 0;
}