hdu4828（卡特兰数+逆元）

这题的前几个数据分别为1,2,5,14,32......................然后确定这是个卡特兰数列

下面来介绍下卡特兰数,它的递推式为f[i+1] = f[i]*(4*n - 6)/n，其中f[2] = f[3] =1;f[4] = 2;f[5] = 14;f[6] = 32..................................

但是这题的n太大了，所以要用到逆元，这里除以一个数，就相当于乘以一个数的逆元。

代码如下（附注释）：

#include<iostream>
#include<algorithm>
#include<cstring>
#include<stack>
#include<queue>
#include<set>
#include<map>
#include<stdio.h>
#include<stdlib.h>
#include<ctype.h>
#include<time.h>
#include<math.h>

#define ll long long
#define N 1000010
#define MOD 1000000007
#define inf 0x7fffffff
#define eps 1e-9
#define pi acos(-1.0)
#define P system("pause")
using namespace std;
ll f[N];

void exgcd(ll a,ll b,ll& d, ll &x, ll &y)//扩展欧几里得算法，模板
{
    if(!b) {d = a; x = 1; y = 0;}
    else
    {
        exgcd(b,a%b,d,y,x);
        y -= x*(a/b);
    }
}

ll inv(ll a, ll n)//求a模n的逆元，也是模板
{
    ll d, x, y;
    exgcd(a,n,d, x, y);
    return d == 1 ? (x+n)%n : -1;
}

int main()
{
//freopen("input.txt","r",stdin);
//freopen("output.txt","w",stdout);
    ios::sync_with_stdio(false);
    ll i;
    f[3] = 1;//f[4] = 2;
    for(i = 3; i < N-5; i++){
         f[i+1] = (4*i-6)%MOD;
         f[i+1] = f[i+1]*f[i]%MOD;
         ll k = inv(i,MOD);
         f[i+1] = (f[i+1]*k)%MOD;
    }
    ll t, z = 1;
    cin>>t;
    while(t--)
    {
        ll n;
        cin>>n;
        printf("Case #%d:\n",z++);
        cout<<f[n+2]<<endl;

    }
    return 0;
}