算法|数学与数论|逆元

数学与数论|逆元

1.费马小定理
2.线性递推
3.扩展欧几里得（exgcd）

心有猛虎，细嗅蔷薇。你好朋友，这里是锅巴的C\C++学习笔记，常言道，不积跬步无以至千里，希望有朝一日我们积累的滴水可以击穿顽石。

线性同余方程(模线性方程)：
线性即未知数都是一次，是最基本的方程，形如ax ≡ b(mod m),其中a,b为整数,m为正整数，x为未知数。该方程可以变形为：ax +by = m,当gcd(a,b)|m时方程有解。

逆元：
如果一个线性同余方程ax ≡ 1(mod b),则x称为a在模b意义下的逆元，记作a^-1。一个数x的逆元通常表示为inv(x)。
逆元存在充分必要条件：a和p（模数）互质：gcd(a,p)=1,此时逆元唯一存在。

费马小定理

定义：
若p为素数，a，p互质：gcd(a,p)=1,则a^(p-1)≡1(mod p)。也有表示为a^p≡a(mod p)
(注：费马小定理必须保证模数p为素数才能使用)

推导：
设x=a^-1，因为ax≡1(mod p)，且a^(p-1)≡1(mod p)
则a^(p-1)≡ax(mod p)
两边同乘x(a^-1),x≡a^(p-2)(mod p)
即 a的逆元在模p意义下与a^(p-2)是同模的，所以我们可以直接用快速幂求 a 的逆元 a^-1。

实践代码：

const int mod = 1e9+7;//模
int qpow(int a,int b){//快速幂
    int ans=1;
    while(b){
        if(b&1) ans*=a,ans%=mod;
        a*=a,a%=mod;
        b>>=1;
    }
    return ans%mod;
}
int inv(int x,int mod){//逆元
    return qpow(x,mod-2);
}
void solve(){
   int x;cin>>x;
   cout<<inv(x,mod);
}

线性递推

费马小定理和exgcd都可以快速地在O(logn)的复杂度内计算出单个值的逆元，但在求解一大段连续值的逆元时，很轻易就会超时，因此我们需要找一种方法线性计算逆元。

递推公式：inv[i]=-p/i*inv[p%i] (p为模数)
证明：
设t=p/i(向下取整)，k=p%i -> p= i * t + k -> i * t + k ≡ 0(mod p) -> i * t ≡ - k(mod p)
i^-1和k^-1分别为i和k对于模p的逆元，上式两边同时乘以i^-1k^-1得 i * t * i^-1 k^-1 ≡ -k * i^-1 * k^-1(mod p) 即 t * k^-1 ≡ -i^-1(mod p)
用inv来表达上式，即inv[i] ≡ -t * inv[k](mod p)*
将t和k代换，得 inv[i] = -p/i * inv[p%i] %p

实践代码：

const int N =1e5+10;
const int p = 1e9+7;
int inv[N];
void solve(){
    inv[0]=inv[1]=1;
    //通常计算的时候都是正数，所以再在递推公式中加一个p
    for(int i=2;i<=1e5;i++) inv[i]=(p-p/i*inv[p%i]%p)%p;
}

扩展欧几里得（exgcd）

裴蜀定理
$对任意整数a,b,存在一堆整数x，y满足ax+by=gcd(a,b)$
扩展欧几里得算法（EXGCD）
常用于求$ax+by=gcd(a,b)的一组可行解。$

在欧几里得算法最后一步，有$b=0，x=1，y=0，$使得$a\ast 1+0\ast 0=gcd(a,0)。$若b>0，则$gcd(a,b)=gcd(b,amodb)$,
假设存在x,y满足 $b\ast x+(amodb)\ast y=gcd(b,amodb)$
$bx+(amodb)y=bx+(a-b\lfloor$a/b$\rfloor )y=ay+b(x-\lfloor$a/b$\rfloor y)$，
所以令$x^{\prime }=y,y^{\prime }=x-\lfloor$a/b$\rfloor y$，就得到了$a{x}^{\prime }+b{y}^{\prime }=gcd(a,b)$的解。

注意
exgcd求逆元适用于模数不为质数的情况

实践代码：

int exgcd(int a,int b,int &x,int &y)
{
    if(!b)
    {
        x=1,y=0;
        return a;
    }
    int d=exgcd(b,a%b,y,x);
    y-=(a/b*x);
    return d;
}
void solve(){
    int a,b,x,y;//a b取任意值
    cin>>a>>b;
    cout<<(exgcd(a,b,x,y)==1?(x+b)%b:-1)<<endl;
}

心有猛虎，细嗅蔷薇。再见了朋友~