求组合数的几种方法

定义

从 n 个不同元素中取出 m 个元素的所有组合的数目

公式

性质

（1） C ( n , m ) = C ( n , n - m ) = C ( n - 1 , m - 1 ) + C ( n - 1 , m )
（2）

实现方法1（杨辉三角）

杨辉三角第 i 行第 j 列对应 C ( i ，j ) 的结果
用动态规划求从 n 个元素中取出 m 个元素的所有组合个数
c [ i ] [ j ] :在 i 个元素中选 j 个
对于当前状态 i j ，可以由两种情况转移而来：
1、前 i - 1 个元素选了 j 个
1、前 i - 1 个元素选了 j - 1 个
即 c [ i ] [ j ] = c [ i - 1 ] [ j ] + c [ i - 1 ] [ j - 1 ]
预处理 c [ i ] [ 0 ] = c [ i ] [ i ] = 1 ( 左右边界 1 )

代码：

int n,m,c[1010][1010];

int main( )
{
	cin>>n>>m;
	for(int i=0;i<=n;i++)
	{
		c[i][0]=1;
		c[i][i]=1;
	}
	for(int i=1;i<=n;i++)
		for(int j=1;j<i;j++)
			c[i][j]=c[i-1][j]+c[i-1][j-1];
	cout<<c[n][m];
	return 0;
} 

实现方法2（逆元）

乘法逆元 ：( a / b ) % mod = a * ( b ^ ( mod - 2 ) ) % mod， mod为质数
b 的逆元就是 b ^ ( mod - 2 )

代码

ll qpow(ll a,ll b)
{
    ll ans=1;
    while(b)
	{
        if(b&1) ans=ans*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return ans;
}
ll f[maxn],inv[maxn];
ll init()
{
    f[0]=1;
    for(int i=1;i<maxn;i++)
        f[i]=f[i-1]*i%mod;
    inv[maxn-1]=qpow(f[maxn-1],mod-2);
    for(int i=maxn-2;i>=0;i--)
        inv[i]=inv[i+1]*(i+1)%mod;
    return 0;
}
ll c(ll n,ll m)
{
    if (n<m) return 0;
    return f[n]*inv[m]%mod*inv[n-m]%mod;
}

待补。。。