C++求组合数(再取模)

C(n,m) = n!/(n-m)!m!
C++在处理大数除法的时候好像会出现问题,所以用除法有时候会因为精度不够没法得到正确答案.需要把除法变成乘法.
知识储备:①逆元②快速幂
逆元:若a*x%M=1 我们称x是a的逆元,x写做a^-1(相当于整除中的倒数),那假设我们要求a/b%M的值,当b特别大时我们无法得到正确答案,这时候对式子进行一些操作:
a / b % M = a / b * 1 % M = a / b * (b * b^-1) % M --------------------- (b * b^-1 % M = 1)
= a * b^-1 % M
这时候只要求出b^-1的值就可以求出a/b%M的值了.
怎样求一个数的逆元呢?
首先,要求一个数的逆元,需要两个元素,一个是这个数,另一个是模数.
a * a^-1%M=1 若M是质数,根据费马小定理a^-1 = a^M-2 具体解释可以在网上搜到.a^M-2就要用到快速幂了.
既然会求逆元了,怎么求阶乘的逆元呢?
①可以用线性算法将1~n的所有逆元都求出来,设inv[i]是i的逆元.有性质:
inv[i] = (M- M/i) * inv[M%i] % M
然后将每个inv[i]乘起来就是i!的阶乘的逆元

inv[1]=1;
for