[luogu1485 HNOI2009] 有趣的数列 (组合数学 卡特兰数)

[luogu1485 HNOI2009] 有趣的数列 (组合数学 卡特兰数)

传送门

Solution

卡特兰数 排队问题的简单变化 答案为\(C_{2n}^n \pmod p\)
由于没有逆元，只好用分解质因数,易证可以整除

Code

//By Menteur_Hxy
#include <ctime>
#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
#define Re register
#define Ms(a,b) memset(a,(b),sizeof(a))
#define Fo(i,a,b) for(Re int i=(a),_=(b);i<=_;i++)
#define Ro(i,a,b) for(Re int i=(b),_=(a);i>=_;i--)
using namespace std;
typedef long long LL;

inline int read() {
    int x=0,f=1;char c=getchar();
    while(!isdigit(c)) {if(c=='-')f=-f;c=getchar();}
    while(isdigit(c)) x=(x<<1)+(x<<3)+c-48,c=getchar();
    return x*f;
}

const int N=2e6+10;
int tot,ans=1,n,MOD;
bool vis[N];
int pri[N],d[N];

void init(int Lim) {
    Fo(i,2,Lim) {
        if(!vis[i]) pri[++tot]=i;
        for(Re int j=1;i*pri[j]<=Lim&&j<=tot;j++) {
            vis[i*pri[j]]=1;
            if(i%pri[j]==0) break;
        }
    }
}

void Divid(int x,int w) {
    for(Re int i=1;i<=tot&&pri[i]*pri[i]<=x;i++) 
        while(x%pri[i]==0) x/=pri[i],d[pri[i]]+=w;
    if(x>1) d[x]+=w;//!!!
}

LL qpow(LL a,int b) {LL t=1;for(;b;b>>=1,a=a*a%MOD)if(b&1)t=t*a%MOD;return t;}

int main() {
    n=read(),MOD=read();init(n+n);Divid(n+1,-1);
    Fo(i,n+1,n+n) Divid(i,1); Fo(i,1,n) Divid(i,-1);
    Fo(i,1,n+n) ans=1ll*ans*qpow(i,d[i])%MOD;
    printf("%d",ans);
    return 0;
}

posted @ 2018-10-17 20:07 Menteur_Hxy 阅读( ...) 评论( ...) 编辑 收藏