HDU4828 卡特兰数+逆元

Description
　　度度熊最近很喜欢玩游戏。这一天他在纸上画了一个2行N列的长方形格子。他想把1到2N这些数依次放进去，但是为了使格子看起来优美，他想找到使每行每列都递增的方案。不过画了很久，他发现方案数实在是太多了。度度熊想知道，有多少种放数字的方法能满足上面的条件？
Input
　　第一行为数据组数T(1<=T<=100000)。
　　然后T行，每行为一个数N(1<=N<=1000000)表示长方形的大小。
Output
　　对于每组数据，输出符合题意的方案数。由于数字可能非常大，你只需要把最后的结果对1000000007取模即可。
Sample Input
2
1
3
Sample Output
Case #1:
1
Case #2:
5

思路真的没什么思路
百度题解发现所有人都是打表发现这个是个卡特兰数….
好那就卡特兰数吧….

卡特兰数就是c(2n,n)/n+1
递推关系也很简单 katelan[a]=(katelan[a-1]*(4*a-2))/a+1
就这样了但是那个/a+1需要逆元
因为他模过了不知道能不能整除
所以再求逆元
两种求法都能用
但是似乎扩展GCD更快一点….大概500ms吧

#include<iostream>
#include<algorithm>
#include<cmath>
#include<string>
#include<cstdio>
using namespace std;
long long mod=1000000007;
long long katelan[1000001];
void exgcd(long long a,long long b,long long& d, long long &x, long long &y)
{  
    if(!b) {d = a; x = 1; y = 0;}  
    else  
    {  
        exgcd(b,a%b,d,y,x);  
        y -= x*(a/b);  
    }  
}  
long long inv(long long a, long long n)
{  
    long long d, x, y;  
    exgcd(a,n,d, x, y);  
    return d == 1 ? (x+n)%n : -1;  
}//讲道理这个板子我没有看懂....不过能用就好了管他能不能看懂的 这不重要 
long long ksm(long long a, long long n,long long mod)
{
    long long fanhui=1;
    while(n)
    {
        if(n&1)
        {
            fanhui*=a;
            fanhui%=mod;
        }
        a*=a;
        a%=mod;
        n>>=1;
    }
    return fanhui%mod;
}
long long feima(long long a,long long n)
{
    return ksm(a,n-2,n);
}
int main()
{
    katelan[0]=katelan[1]=1;
    for(int a=2;a<=1000000;a++)
    {
        katelan[a]=(4*a-2)*katelan[a-1];
        katelan[a]%=mod;
        katelan[a]*=feima(a+1,mod);
        katelan[a]%=mod; 
    }
    int T;
    cin>>T;
    int u=0;
    long long n;
    while(T--)
    {
        scanf("%lld",&n);
        printf("Case #%d:\n",++u);
        printf("%lld\n",katelan[n]);
    }
    return 0;
}