hdu 5184 Brackets && BestCoder Round #32

题解：1.卡特兰数的模型，维基百科上有证明

令1表示进栈，0表示出栈，则可转化为求一个2n位、含n个1、n个0的二进制数，满足从左往右扫描到任意一位时，经过的0数不多于1数。显然含n个1、n个0的2n位二进制数共有个，下面考虑不满足要求的数目。

考虑一个含n个1、n个0的2n位二进制数，扫描到第2m+1位上时有m+1个0和m个1（容易证明一定存在这样的情况），则后面的0-1排列中必有n-m个1和n-m-1个0。将2m+2及其以后的部分0变成1、1变成0，则对应一个n+1个0和n-1个1的二进制数。反之亦然（相似的思路证明两者一一对应）。

从而。证毕。

红字部分之后 进行一些小小的补充，帮助理解：

每一个变换之前的序列一一对应于变换之后的序列

变换之前的序列个数为n，之后的序列个数为m

也就是说，对于每一个变换前的序列可以唯一的找到变换后的一个序列相对应，n  <= m ;反之亦然m <= n

2.需要用到乘法逆元来进行取模操作：

 (b / a) % m = c,为了求c

a * x = 1(MOD m)(这里本来是应该三条线的等号打不出来。。。。)

第一个式子 (b / a) % m = c  变形后为  (b / a) = c(MOD m)

两式相乘得b * x = c(MOD m)

那么求出x就可以求出c啦

其实x就是a 的逆元

第二个式子a * x = 1(MOD m)     变形后为   a * x - 1 = qm    变形后为 a * x - qm = 1

用欧几里得扩展定理求x便可

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
typedef long long LL;
#define MAXN 1000005
#define MOD 1000000007
LL dp[MAXN];
char s[MAXN];
int n;
LL extgcd(LL a,LL b,LL & x,LL & y)
{
    LL d = a;
    if(b)
    {
        d = extgcd(b,a % b,y,x);
        y -= (a / b) * x;
    }
    else
    {
        x = 1;
        y = 0;
    }
    return d;
}
LL mod_inverse(LL a,LL m)
{
    LL x,y;
    extgcd(a,m,x,y);
    return (m + x % m) % m;
}
LL C(int a,int b)
{
    if(b < 0)return 0;
    LL inverse = (mod_inverse(dp[b],MOD) * mod_inverse(dp[a - b],MOD)) % MOD;
    return (dp[a] * inverse) % MOD;
}
LL solve()
{
    int len = strlen(s),p = 0,q = 0;
    if(n & 1)return 0;
    for(int i = 0;i < len;i++)
    {
        if(s[i] == '(')p++;
        else q++;
        if(p < q)return 0;
    }
    return (C(n - p - q,n / 2 - p) - C(n - p - q,n / 2 - p - 1) + MOD) % MOD;
}
int main()
{
    dp[0] = 1;
    for(int i = 1;i < MAXN;i++)
        dp[i] = (dp[i - 1] * i) % MOD;
    while(scanf("%d",&n) != EOF)
    {
        scanf("%s",s);
        LL ans = solve();
        printf("%I64d\n",ans);
    }
}