卡特兰数 中南大学OJ1320

#include <iostream>
#include <cstdio>
#include <cstring>

using namespace std;

void gcd(int a,int b,int &d,int &x,int &y)
{
    if(!b)
    {
        d=a;
        x=1;
        y=0;
    }
    else
    {
        gcd(b,a%b,d,y,x);
        y-=x*(a/b);
    }
}

//计算模n下a的逆元，如果不存在逆元，返回-1
int inv(int a,int n)
{
    int d,x,y;
    gcd(a,n,d,x,y);
    return d==1? (x+n)%n: -1;
}

long long arr[10010];
long long mod=1000000007;

int main()
{
    int n;
    memset(arr,0,sizeof(arr));
    arr[0]=1;
    for(int i=1;i<=10000;i++)
    {
        arr[i]=((arr[i-1]*2*(2*i+1)%mod)*inv(i+2,mod))%mod;
    }
    while(scanf("%d",&n)!=EOF)
    {
        printf("%lld\n",arr[n-1]);
    }
    return 0;
}