杨辉三角形

\quad \quad \quad \quad $C_n^m=\frac{n!}{m!(n-m)!}$

组合数和杨辉三角有着密切的关系。杨辉三角上的每一个数字都等于它的左上方和右上方的和（除了边界）

第$n$行，第$m$个就是，就是$C_n^m$ （从$0$开始）

所以以后求杨辉三角或者组合数都可以用到下面的递推公式：

Code:

const int N = 2000 + 5;
const int MOD = (int)1e9 + 7;
int comb[N][N];//comb[n][m]就是C(n,m)
void init()
{
    for(int i = 0; i < N; i ++)
    {
        comb[i][0] = comb[i][i] = 1;
        for(int j = 1; j < i; j ++)
        {
            comb[i][j] = comb[i-1][j] + comb[i-1][j-1];
            comb[i][j] %= MOD;
        }
    }
}

不过这个是$O（n^2）$，有没有更简单的呢？肯定有，因为我们学了逆序数，所以我们可以直接求解：

Code:

ll fac[N], invn[N], invfac[N];

int init()
{
    fac[0] = fac[1] = invfac[0] = invfac[1] = invn[0] = invn[1] = 1;
    for (int i = 2; i <= 5000005; ++i)
    {
        fac[i] = fac[i - 1] * i % mod;
        invn[i] = (mod - mod / i) * invn[mod % i] % mod; //逆元线性递推公式
        invfac[i] = invfac[i - 1] * invn[i] % mod;
    }
}

ll comb(ll n, ll m)
{
    if (m > n)
        return 0;
    if (m < 0 || n < 0)
        return 0;
    ll res = fac[n];
    res = res * invfac[n - m] % mod;
    res = res * invfac[m] % mod;
    return res;
}