逆元 入门?

零、共识

$a\mid b$：a能整除b( $b\%a=0$ )
$a\nmid b$：a不能整除b( $b\%a\ne 0$ )
$a\equiv b(modm)$：a与b模n同余($a\%n=b\%n$)
$(a,b)=gcd(a,b)$
$[a,b]=lcm(a,b)$
$inv(a)=\frac{1}{a}$
充分必要条件： $a\iff b$

一、定义及性质

1 什么是逆元

如果一个线性同余方程 $ax\equiv 1(modb)$，则称 $x$ 为 $amodb$的逆元，记作 $a^{-1}$ 。

2 整除的性质

1. $a\mid b\iff -a\mid b\iff a\mid -b\iff |a|\mid |b|$
2. $a\mid b,b\mid c\Rightarrow a\mid c$
3. $a\mid b,a\mid c\Rightarrow a\mid (bx+cy)$
4. $a\mid b,m\ne 0\Rightarrow ma\mid mb$
5. $m\mid ab,(a,b)=1\Rightarrow m\mid b$
6. $a_i\mid b(1\le i\le n)\iff [a_1,a_2\dots a_n\mid b$

3 同余的性质

1. 当 $a\equiv b(modb)$ 时：

$a\equiv a(modb)$ 自反

$b\equiv a(modb)$ 对称

$b\equiv c(modb)$ $\Rightarrow a\equiv c(modb)$ 传递
$a+c\equiv b+c(modb)$ 同加
$a\times c\equiv b\times c(modb)$ 同乘
$a^c\equiv b^c(modb)$ 同幂

2. $amodp=x,amodq=x,(p,q)=1\Rightarrow amodpq=x$

3. $d\ge 1,d\mid m,a=b(modm)\Rightarrow a\equiv b(modd)$

4. $d\ne 0\Rightarrow a\equiv b(modm)\iff da\equiv db(mod|d|m)$

5. $ca\equiv cb(modm)\iff a\equiv b(mod\frac{m}{(c,m)})$

4 费马小定理

如果$p$ 是一个质数，而整数 $a$ 不是 $p$ 的倍数，则$a^{p-1}\equiv 1(modp)$。

5 逆元

根据费马小定理：
$\because a^{p-1}\%p=1$
$\therefore (a\times a^{p-2})\%p=1$
$\because 1\%p=1$
$\therefore (a\times a^{-1})\%p=1$
$\therefore a^{-1}=a^{p-2}$

二、实践

1 求逆元

1.1 题目

求 $b=10^9+7$ 时 $amodb$ 的逆元（ $a^{-1}$ ）。

1.2 思路

直接用 $a^{-1}=a^{p-2}$ 即可。

1.3 代码

#include<bits/stdc++.h>
#pragma GCC optimize(2)
#define ll long long
#define bug printf("---OK---")
#define pa printf("A: ")
#define pr printf("\n")
#define pi acos(-1.0)
using namespace std;
ll a;
const ll mod=1e9+7;
ll qpow(ll a,ll b){
    ll ans=1;
    while(b){
        if(b&1){
            ans=ans*a%mod;
        }
        b>>=1;
        a=a*a%mod;
    }
    return ans;
}
int main(){
    cin>>a;
    cout<<qpow(a,mod-2);
    return 0;
}

2 线性求逆元

2.1 题目

给定 $n,p$ ，求 $1$ 到 $n$ 所有整数在模 $p$ 意义下的乘法逆元。

2.2 思路

以下除法都是向下取整。设:
$$k=\frac{p}{i},j=pmodi$$
能够得到：
$$p=ki+j$$
带入同余式中：
$$p\equiv ki+j\equiv 0(modp)$$
$$ki+j\equiv 0(modp)$$
两边同时乘以i的逆元,得到:
$$k{j}^{-1}+i^{-1}\equiv 0(modp)$$
移项：
$$i^{-1}\equiv -k{j}^{-1}(modp)$$
也就是：
$$i^{-1}\equiv -\frac{p}{i}\equiv (p-\frac{p}{i})j^{-1}\equiv (p-\frac{p}{j})(pmodi{)}^{-1}(modp)$$

2.3 代码

#include<bits/stdc++.h>
#pragma GCC optimize(2)
#define ll long long
#define bug printf("---OK---")
#define pa printf("A: ")
#define pr printf("\n")
#define pi acos(-1.0)
using namespace std;
ll n,p;
ll inv[3000002];
int main(){
	cin>>n>>p;
	inv[1]=1;
	cout<<1;pr;
	for(int i=2;i<=n;i++){
		inv[i]=(p-p/i)*inv[p%i]%p;
		cout<<inv[i];pr;
	}
	return 0;
}

3 逆元求组合数

3.1 题目

求组合数 $C_n^m$ 对 $10^9+7$ 取模的值。

3.2 思路

首先初始化 $1\%mod$ ~ $10^6\%mod$ 的阶乘。由于组合数 $C_n^m=\frac{n!}{((n-m)!\times m!)}$ ，每次处理时使用公式$a^{-1}=a^{p-2}$ 即可，注意取模。

3.3 代码

#include<bits/stdc++.h>
#pragma GCC optimize(2)
#define ll long long
#define bug printf("---OK---")
#define pa printf("A: ")
#define pr printf("\n")
#define pi acos(-1.0)
using namespace std;
ll fac[1000005];
const ll mod=1e9+7;
ll qpow(ll a,ll b){
    ll ans=1;
    while(b){
        if(b&1){
            ans=ans*a%mod;
        }
        b>>=1;
        a=a*a%mod;
    }
    return ans;
}
void init() {
    fac[0]=1;
    for(int i=1;i<=1000000;i++) {
        fac[i]=fac[i-1]*i%mod;
    }
}
ll C(ll n,ll m) {
    return fac[n]*qpow(fac[m]*fac[n-m]%mod,mod-2)%mod;
}
int main(){
    init();
    int t;
    cin>>t;
    while(t--) {
        ll n,m;
        cin>>n>>m;
        cout<<C(n, m);pr;
    }
    return 0;
}