Três Métodos Rápidos de Elemento Inversor Fatorial Teorema de Fermat Recursão Linear
1. O pequeno teorema de Fermat
#include<cstdio>
typedef long long LL;
const LL MOD = 1e9 + 7;
LL fac[1000000+5]; //阶乘
LL inv[1000000+5]; //逆元
LL quickMod(LL a,LL b)
{
LL ans = 1;
while (b)
{
if (b&1)
ans = ans * a % MOD;
a = a*a % MOD;
b >>= 1;
}
return ans;
}
void getFac()
{
fac[0] = inv[0] = 1;
for (int i = 1 ; i <= 1000000 ; i++)
{
fac[i] = fac[i-1] * i % MOD;
inv[i] = quickMod(fac[i],MOD-2); //表示i的阶乘的逆元
}
}
LL getC(LL n,LL m) //C(n,m) = n!/((n-m)!*m!) % (1e9+7)
{
return fac[n] * inv[n-m] % MOD * inv[m] % MOD;
}
int main()
{
getFac();
int n,m;
while (~scanf ("%d %d",&n,&m))
printf ("%lld\n",getC((LL)n,(LL)m));
return 0;
}
2. Encontre recursivamente o fatorial inverso
Podemos considerar o uso do pequeno teorema de Fermat para encontrar primeiro o inverso do maior fatorial e depois empurrá-lo para trás, olhar diretamente para o código e explicar.
void init() {
fact[0] = 1;
for (int i = 1; i < maxn; i++) {
fact[i] = fact[i - 1] * i %mod;
}
inv[maxn - 1] = power(fact[maxn - 1], mod - 2);
for (int i = maxn - 2; i >= 0; i--) {
inv[i] = inv[i + 1] * (i + 1) %mod;
}
}
3. Inversão linear
Portanto, a recursão linear é
inv[i]=(mod-mod/i)*inv[mod%i]%mod;
Como usá-lo? Podemos usar isso para encontrar o valor do número combinado no sentido de mod, a função é a seguinte:
const ll mod=1e9+7;
const ll maxn=1e5+1;
ll f[maxn]={1,1};
ll f0[maxn]={1,1};
ll inv[maxn]={1,1};
void init()
{
for(int i=2;i<=maxn;i++)
{
//阶乘数
f[i]=f[i-1]*i%mod;
//i在mod意义下的逆元
f0[i]=(mod-mod/i)*f0[mod%i]%mod;
//阶乘逆元
inv[i]=inv[i-1]*f0[i]%mod;
}
}
//求阶乘数C(a,b)在mod意义下的值
ll C(ll a,ll b)
{
return f[b]*inv[a]%mod*inv[b-a]%mod;
}