来看排列组合嘛

保证是你见过最简洁明了的排列组合文章

费马小定理

如果p是一个质数，而整数a不是p的倍数，则有$a^{p-1}\equiv 1$ (mod p）
证明：略

乘法逆元

简单的可以理解成一个数 对p取模 的倒数

证明：

若b、m互质，并且b | a，则存在一个整数x，使得a / b ≡ a * x (mod m)，则称x为b的模m乘法逆元，记为$b^{-1}$(mod m)。

因为$a/b=a\ast b^{-1}\equiv a/b\ast b\ast b^{-1}$，所以$b\ast b^{-1}\equiv 1$(mod m)

如果m是质数，（后文用p代表m）并且b < p,根据费马小定理$b^{p-1}\equiv 1$(mod p)，
即$b\ast b^{-2}\equiv 1$(mod p)。
因此，当模数p为质数的时候，$b^{p-2}$即为b的乘法逆元

结论：

• 当模数p与b互质的时候，$b^{p-2}$即为b的乘法逆元
• 通常用于除法中求模

如果p不是素数的情况，需要求逆元：
$a/b$ (mod m) = $x/d$(mod $d\ast m$)

组合的简单的写法

看别人的没看明白，顺着自己思路写了
也没什么好讲的，就是公式的写法。

组合的公式： $C_n^m=\frac{n!}{m!(n-m)!}$

根据上文的逆元
公式可以改写成：
$C_n^m=\frac{n!}{m!(n-m)!}=n!\ast inv[m]\%mod\ast inv[n-m]\%mod$
此处inv[m]已经是m阶乘的inv

问题如下：
给定n，m 求$C_n^m$，结果对$10^9+7$取模

思路是打表 c[n] = n!，以及对应的逆元inv[n]。然后输出结果

#include <iostream>
#define ll long long
using namespace std;
const int N = 1e5+5,mod= 1e9+7;
int a,b,t;
ll c[N], inv[N];

ll qpow(ll a,int b){
    ll ans = 1;
    while(b){
        if(b & 1) ans = ans * a % mod;
        a = a * a % mod;
        b >>= 1;
    }
    return ans;
}
inline ll C(ll a,ll b,ll p){
	return c[a] * inv[b] % p * inv[a-b] % p;
	 //  a / b  % mod == a * inv[b] % mod;
}
int main(){
    c[0] = 1;
    for(int i=1;i<=N;i++) c[i] = c[i-1] * i % mod;//  求排列数
    
    inv[N-1] = qpow(c[N-1], mod - 2);//  求排列数的逆元	b的逆元就是 b^(mod-2)
    for(int i = N-2;i >= 0;i--)	inv[i] = inv[i + 1] * (i + 1) % mod;
    
    cin >> t;
    while(t--){
        cin >> a >> b;
        cout << C(a, b, mod) << "\n";
    }
    return 0;
}

Lucas定理

求值 $C_a^b$ mod p
适用于：数论中求大组合数取模.
比如：
1 ≤ b ≤ a ≤ $10^{18}$
1≤ p ≤ $10^5$

Lucas定理完整内容：
$C_n^m=\prod _{i=0}^k{C}_{n_i}^{m_i}$(mod p)
$n=n_k{p}^k+n_{k-1}{p}^{k-1}+...+n_{1}p+n_0$
$m=m_k{p}^k+m_{k-1}{p}^{k-1}+...+m_{1}p+m_0$

简易写法：
$C_a^b\equiv C_{a\%p}^{b\%p}\ast C_{a/p}^{b/p}\%p$

其中$C_{a/p}^{b/p}$可以递归

那说起递归可能你比我还懂啊。

代码

ll lucas(ll a, ll b, ll p){
    if(a < p && b < p) return C(a, b, p);
    return C(a % p, b % p, p) * lucas(a / p, b / p, p) % p;
}