#求组合数

求组合数（一）

1. 适用思路：范围1e3左右
2. 时间复杂度：n^2

算法思路：

• 首先引入一个关于组合排列的数学公式：Cnm = Cn-1,m + Cn-1,m-1
• 通过这个组合数的公式(杨辉三角系数)来线性递推出所求组合数的值

模板CODE

---

#include<iostream>
#include<algorithm>
using namespace std;
const int mod=1e9+7;
int dp[2080][2080];//存组合数的数值
int n,m;
int t;
void init()
{
    for(int i=0;i<2080;i++)
    {
        for(int j=0;j<=i;j++)
        {
            if(!j) dp[i][j]=1;//初始化Cn0
            else
            {
                dp[i][j]=(dp[i-1][j-1]+dp[i-1][j])%mod;//组合数公式
            }
        }
    }
}
int main()
{
    init();
    cin>>t;
    while(t--)
    {
        cin>>n>>m;
        cout<<dp[n][m]<<endl;//回答查询
    }
    return 0;
}

---

求组合数（二）

1. 适用思路：范围1e6左右
2. 时间复杂度：O(n)

算法思路：

• 在数学表达式上n的阶乘的逆元=（n+1)的阶乘的逆元乘以(n+1)

模板CODE

---

#include<iostream>
#include<cstring>
#define endl '\n'
const int mod = 1e9 + 7, N = 1e5;
using namespace std;
int fact[N+10], invf[N+10];
// C(n,m) = n! /  m!*(n-m)!

int qmi(int m, int n = mod-2)
{
	int res = 1;
	while(n)
	{
		if(n & 1) res = 1ll*res*m%mod;
		m = 1ll*m*m%mod;
		n >>= 1;
	}
	return res;
}

int comb(int n, int m)
{
	return 1ll*fact[n]*invf[m]%mod*invf[n-m]%mod;
}

void solve()
{
	int n, m; cin >> n >> m;
	cout << comb(n, m) << endl;
}

int main()
{
	int t; cin >> t;
	fact[0] = invf[0] = 1;
	for(int i = 1; i <= N; i ++)
	{
		fact[i] = 1ll*fact[i-1]*i%mod;
	}
	invf[N] = qmi(fact[N]);
	cout << invf[N] << endl;
	for(int i = N-1; i >= 1; i --)
	{
		invf[i] = 1ll*invf[i+1]*(i+1)%mod;
	}

	while(t --)
	{
		solve();
	}
	return 0;
}

---