组合数计算与模运算优化技巧

1. 题目背景与核心考察点

1656题"Combination"是《信息学奥赛一本通》中组合数学章节的经典题目。这道题主要考察选手对组合数计算的理解以及处理大数取模问题的能力。在实际比赛中，这类题目经常作为中等难度题出现，需要选手掌握模运算下的组合数计算技巧。

组合数C(n,m)表示从n个不同元素中取出m个元素的方案数，其计算公式为：
C(n,m) = n! / (m! * (n-m)! )

在算法竞赛中，由于n和m的取值可能很大（比如1≤m≤n≤10^6），直接计算阶乘会导致数值溢出，因此必须采用模运算下的组合数计算方法。

2. 常用解法与算法选择

2.1 暴力计算法的局限性

最直观的方法是直接套用组合数公式，预先计算各个数的阶乘，然后进行除法运算。这种方法在小数据范围时可行，但当n超过20时，阶乘值会迅速超过普通数据类型的表示范围（如C++中long long类型最大约9×10^18，而21! ≈ 5.1×10^19）。

即使采用高精度算法，时间复杂度也会变得不可接受。因此，在算法竞赛中，我们通常需要考虑模运算下的解决方案。

2.2 模运算下的组合数计算

题目通常会给出一个质数P（常见的有1e9+7、998244353等），要求输出C(n,m) mod P的结果。这种情况下，我们可以利用数论中的逆元概念来避免直接除法运算。

核心思路是将除法转换为乘法逆元：
a / b ≡ a * b^(-1) (mod P)

其中b^(-1)表示b在模P下的乘法逆元，即满足b * b^(-1) ≡ 1 (mod P)的整数。

3. 预处理阶乘与逆元的优化方法

3.1 费马小定理求逆元

当P是质数时，根据费马小定理：
b^(P-1) ≡ 1 (mod P)

因此，b的逆元可以表示为：
b^(-1) ≡ b^(P-2) (mod P)

这样我们就可以用快速幂算法在O(logP)时间内计算单个数的逆元。

3.2 线性预处理阶乘和逆元

为了优化多次查询的效率，我们可以预先计算0到n的阶乘数组fact[]及其逆元数组inv_fact[]：

1. 计算阶乘数组：
   fact[i] = fact[i-1] * i % P

2. 计算逆元数组：
   利用性质：(i!)^(-1) ≡ ((i+1)!)^(-1)