HNOI2009有趣的数列

首先next_permutation打表，发现Cat规律。

其实考试的时候这么作没什么问题，并且能够节省异常多的时间，那么如今咱们来想一下why。

首先我拿模型法解释一下，咱们把2n个数当作2n我的，既然分红奇数和偶数两种比较方式，那么我让他们站成两排，每一排有n我的，这n我的的身高递增，且，第二排的人必须高于第一排，那么这个问题就变成了：

有2n个身高互不相同人站成两排，每排n人，要求右边的人比左边的人高，后面的人比前面的人高，问我有几种排队方案。

这是一个Cat的模型，既然先站哪一排无所谓，我就让某个位置必须先站上第一排的人再站上第二排的人，若是我将站在第一排看作是0，站在第二排看作是1，那么既然每一个1前面必定有一个比他矮的人，则必定有一个0，那么就又转化成了求0，1序列，这是一个更加经典的Cat模型。（若是这里理解不了能够上网搜搜）

而后再拿折线法解释一下，咱们把偶数项看作x轴上的数，由于他们是单增的，把奇数项看作y轴上的数，因为奇数项小于与之对应偶数项，也就是不能越过y=x，函数的变化就好像只能向右走和向上走。这个问题在上一篇博客中有详细的解法。

因此咱们明白它是让咱们求Cat，但是P不必定是质数，逆元的问题很恶心。

因此咱们采用惟一分解来作。首先线性筛筛出2n之内的全部素数，而后咱们枚举每一个素数，对n执行如下操做：将n不断的除以这个素数，并将商加入s变量，最终s的值就是n！在算术基本定理拆分后，这个素数的指数。举个例子：

8!=2⁷*3²*5*7，8/2=4，4/2=2，2/2=1，1/2=0。4+2+1+0=7。

20!=2¹⁸……，20/2=10，10/2=5，5/2=2，2/2=1，1/2=0。10+5+2+1+0=18。

你们能够本身随便试两个。

这是为何呢？（下述i为质数）首先1~n中含有i这个因子的数有n/i个（1），含有i²这个因子的数有n/i²个（2），……含有i^m这个因子的数有n/i^m个（m）。那么咱们分层计算贡献，首先（1）中有n/i个i，加上，（2）中有2*n/i²个i，但不要忘了，咱们在（1）算过每一个数中的一个i，那么它们的贡献只有n/i²个i，同理，向后类推，最后n！中i的个数为∑n/pⁱ，与上述模拟过程一致。

那么咱们来证实一下复杂度，首先根据小于N的质数约有N/lnN个，咱们第一层枚举的代价就是O（N/lnN），而后观察上述过程，咱们的问题规模不断缩小，如上述二例，都是1/二、1/2的速度在缩小，对于其余素数相似，咱们取最坏O（log₂N）,那么总复杂度

O（N/lnN*log₂N），这玩意换换底就是O（N/ln2），1/ln2≈1.44，撇掉，大约O（N），（这是我本身证的，网上目测没有，若是有异议请指出，应该没什么问题吧……）

而后分子加分母减拆完了拿快速幂一乘就完事了。（快速幂并不影响上述复杂度，由于qpow也是O（logk）的，就当常数大了吧）。

（底下代码有表机，勾掉的调试略多，能够用来本身见证一下上面那个算法的正确性）

#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<cstdio>
#include<vector>
#include<queue>
#include<stack>
#include<set>
#include<map>
using namespace std;
int n,P,/*a[20],*/ans=1;
/*bool check(){
    for(int i=1;i<=n;i++)
        if(a[i*2-1]>a[i*2]) return 0;
    for(int i=3;i<=n*2;i++)
        if(a[i]<a[i-2]) return 0;
    return 1;
}*/
int prime[6000000],prime_num;
bool v[20050000];
void doprime(){
    for(int i=2;i<=2*n+5;i++){
        if(!v[i]) prime[++prime_num]=i;
        for(int j=1;j<=prime_num&&i*prime[j]<=2*n+5;j++){
            v[prime[j]*i]=1;
            if(i%prime[j]==0) break;
        }
    }
}
int qpow(int x,int k){
    int val=1;
    for(;k;k>>=1,x=1ll*x*x%P)
        if(k&1) val=1ll*val*x%P;
    return val%P;
}
int main(){
    //打表找规律系列。。。
/*    while(1){
        ans=0;
    scanf("%d%d",&n,&P);
    for(int i=1;i<=(n<<1);i++)
        a[i]=i;
    do{
        if(check()) {ans++;
            for(int i=1;i<=2*n;i++)
                cout<<a[i]<<" ";
            cout<<endl;
        }
    }while(next_permutation(a+1,a+1+2*n));
    printf("ANS=%d\n",ans);
    }*/
    scanf("%d%d",&n,&P);
    doprime();
    for(int i=1;i<=prime_num;i++){
        long long s=0;
        for(int j=2*n;j/=prime[i];) s+=j;
    //    cout<<"s1="<<s<<endl;
        for(int j=n;j/=prime[i];) s-=j;
        //cout<<"s2="<<s<<endl;
        for(int j=n+1;j/=prime[i];) s-=j;
    //    cout<<"s3="<<s<<endl;
        ans=1ll*ans*qpow(prime[i],s)%P;
    }
//    cout<<"Okprime"<<endl;
/*    for(int i=1;i<=prime_num;i++)
        cout<<prime[i]<<" ";cout<<endl;*/
/*    for(int i=1;i<=2*n;i++)
        mulfz(i);
    for(int i=1;i<=n;i++)
        mulfm(i);
    for(int i=1;i<=n+1;i++)
        mulfm(i);*/
/*    cout<<"OKfenjie"<<endl;
    for(int i=1;i<=prime_num;i++)
            ans=1ll*ans*qpow(prime[i],fz[i]-fm[i])%P;
    cout<<"Okqpow"<<endl;*/
    printf("%d",ans);
    return 0;
}

View Code

这道题取模，下道题高精。