扩欧(解线性同余方程)+逆元

1.求解线性同余方程

定理1： $gcd(a,b)=gcd(b,a\%b)$

int gcd(int a,int b){
    if(b==0) return a;
    else return gcd(b,a%b);
}//欧几里得算法

定理2： 对于方程$a\ast x+b\ast y=c$,该方程等价于$a\ast x\equiv cmod\left(b\right)$，方程有整数解的充分必要条件是：$c\%gcd(a,b)=0$

1.求$ax+by=\gcd (a,b)$通解

因为欧几里得算法结束时变量a存放gcd，变量b存放0，因此此时显然有$$a\ast 1+b\ast 0=\gcd$$
我们采用递归求解，由数学知识得到递推公式：（若推不出来，望自学）
$$\{\begin{array}{l}{x}_1=y_2\\ y_1=x_2-\left(a/b\right)y_2\end{array}$$
之后通过x2,y2来反推x1,y2即可（就是目前已知递归边界成立时为x=1,y=0,需想办法反推出最初始的x和y）

int exgcd(ll a,ll b,ll &x,ll &y){
    if(!b){x=1, y=0; return a;}
    gc=exgcd(b,a%b,y,x);
    y-=a/b*x;
    return gc;//gc即是gcd(a,b)
}//递归求ax+by=gcd(a,b) 

通解为：
$$\begin{array}{r}\{\begin{array}{l}{x}^{\prime }=x+\frac{b}{\gcd }\cdot k\\ y^{\prime }=y-\frac{a}{\gcd }\cdot k\end{array}\\ \left(k\in Z\right)\end{array}$$

2.求$a\ast x+b\ast y=c$通解

根据上述定理，对于方程$a\ast x+b\ast y=c$，我们可以先用扩展欧几里得算法求出一组解$(x_0,y_0)$,使其满足$a{x}_0+b{y}_0=\gcd$，然后两边同时除以$gcd(a,b)$，再乘以$c$，这样就得到了方程:
$$\frac{a\cdot x_0\cdot c}{gcd\left(a,b\right)}+\frac{b\cdot y_0\cdot c}{gcd\left(a,b\right)}=c$$
我们也就找到了原方程的一组解$(x,y)$的表示形式:
$$\left(\frac{x_0\cdot c}{gcd\left(a,b\right)},\frac{y_0\cdot c}{gcd\left(a,b\right)}\right)$$，而且和全部解和$a{x}_0+b{y}_0=\gcd (a,b)$的通解类似，为:
$$\begin{array}{r}\{\begin{array}{l}{x}^{\prime }=x+\frac{b}{\gcd }\cdot k\\ y^{\prime }=y-\frac{a}{\gcd }\cdot k\end{array}\\ \left(k\in Z\right)\end{array}$$

特例： 若$gcd(a,b)=1$，且$x_0,y_0$为 $a\ast x+b\ast y=c$ 的一组解，则该方程的任一解可表示为：$x=x_0+bt，y=y_0-at$，且对任一整数$t$，皆成立。

观察全部解的形式，可以求出方程的所有解。但实际问题中，我们往往被要求去求最小非负整数解，也就是求一个特解$x$，方法为：$t=b/gcd(a,b)，x=(x\%t+t)$

if(c%gcd(a,b)==0){//判断是否有解
	ll gc=exgcd(a,b,x,y);
	x=c/gc*x;//x0转x
	ll t=b/gc;//求x的周期 
	x=(x%t+t)%t;//x的最小正整数解
	if(c-x*a>=0){//需要保证y也是正整数解 
	    printf("");
	    return 0;
	}
}

3.例题分析：斩杀线计算大师

例题：斩杀线计算大师
思路：枚举所有a卡可能的值，使问题转化为求解二元一次方程

 __gcd(a,b)函数用来求a和b的最大公因数，需要#include<algorithm>

ac代码：

#include<bits/stdc++.h>
#define ll long long
using namespace std;
ll gcd(ll a,ll b){
    if(!b)return a;
    return gcd(b,a%b);
}
ll exgcd(ll a,ll b,ll &x,ll &y){
    if(!b){x=1;y=0;return a;}
    int gc=exgcd(b,a%b,y,x);
    y-=a/b*x;
    return gc;
}//递归求ax+by=gcd(a,b)
int main(){
    ll a,b,c,k;
    scanf("%lld%lld%lld%lld",&a,&b,&c,&k);
    ll t=k/c,c1,x,y;
    for(ll i=0;i<=t;i++){
        c1=k-i*c;
        if(c1%gcd(a,b)==0){//判断是否有解
            ll gc=exgcd(a,b,x,y);
            x=c1/gc*x;//x0转x
            ll t=b/gc;//求x的周期 
            x=(x%t+t)%t;//x的最小正整数解
            if(c1-x*a>=0){//注意,x是正整数解但也要保证y是正整数解 
                printf("%lld %lld %lld",x,(c1-x*a)/b,i);
                return 0;
            }
        }
    }
    return 0;
}

2.逆元：

1.定义

假设 $a、b、m$ 是整数，$m>1$，且有 $ab\equiv 1(modm)$ 成立，那么就说 $a$ 和 $b$ 互为模 $m$ 的逆元，一般也记作 $a\equiv 1/b(modm)$ 或 $b\equiv 1/a(modm)$ 。通俗地说，如果两个整数的乘积模 $m$ 后等于 $1$ ，就称它们互为 $m$ 的逆元

逆元用处：

对乘法来说有$$（b\times a）\%m=（（b\%m）（a\%m）\%m$$
成立，但是对除法来说，$$（b/a）\%m=（（b\%m）/（a\%m））\%m$$
却不成立，$$（b/a）\%m=（（b\%m）/a）\%m$$也不成立，这时就需要逆元来计算$$（b/a）\%m$$
通过找到 $a$ 模 $m$ 的逆元，就有$$（b/a）\%m=（b\ast x）\%m$$成立 （证明过程略）

只考虑整数取模，也即假设$b\%a=0$，即b是a的整数倍，于是就把除法取模转化为乘法取模，这对于解决被除数b非常大（使得b已经取过模，不是原始值）的问题来说是非常实用的。

逆元的含义： 模$n$意义下，一个数$a$如果有逆元$x$，那么除以$a$相当于乘以$x$

由定义知，求$a$模$m$的逆元，就是求解同余式$ax\equiv 1(modm)$，并且在实际使用中，一般把$x$的最小正整数解称为$a$模$m$的逆元，因此逆元经常被默认为$x$的最小正整数解。

2.扩欧求解

显然，同余式$ax\equiv 1(modm)$是否有解取决于$1\%gcd(a,m)$是否为0，而这等价于$gcd(a,m)$是否为1：

①如果$gcd(a,m)\ne 1$，那么同余式$ax\equiv 1(modm)$无解，$a$不存在模$m$的逆元。
②如果$gcd(a,m)=1$，那么同余式$ax\equiv 1(modm)$在$(0,m)$上有唯一解，可以通过求解$ax+my=1$得到。

注意：由于$gcd(a,m)=1$，因此$ax+my=1=gcd(a,m)$，直接使用扩展欧几里得算法解出$x$之后就可以用$(x\%m+m)\%m$得到$(0,m)$范围内的解，也就是所需要的逆元。
下面的代码使用了扩展欧几里得算法来求解$a$模$m$的逆元，使用条件是$gcd(a,m)=1$
（我们不关心$m$是不是素数）

int exgcd(ll a,ll b,ll &x,ll &y){
    if(!b){x=1, y=0; return a;}
    gc=exgcd(b,a%b,y,x);
    y-=a/b*x;
    return gc;//gc即是gcd(a,b)
}//递归求ax+by=gcd(a,b) 
int inv(int a,int m){
    int x,y;
    int g=exgcd(a,m,x,y);//求ax+my=1 
    return (x%m+m)%m;//得到(0,m)范围内的解 
}//求逆元

3.费马小定理求解

费马小定理： 设$m$是素数，$a$是任意整数且$a\not\equiv 0\left(modm\right)$,则$a^{m-1}\equiv 1\left(modm\right)$

在此借鉴大佬的解释了（神崎兰子）

4.例题分析：抽卡

例题：抽卡

分析：

先求出全部抽完都没有想要的卡的概率，然后用1减去这个概率就可以了，全过程通过逆元把分数当成整数进行运算即可。

#include<bits/stdc++.h>
#define mod 1000000007
//int的最大值为2147483647 
using namespace std;
int n,a[100005],b[100005],ans;
int read(){
    int x=0;char ch=getchar();
    while(!isdigit(ch)) ch=getchar();
    while(isdigit(ch)) x=(x<<3)+(x<<1)+(ch^48),ch=getchar();
    return x;
}
int ksm(int x,int y){
    int res=1;
    while(y){
        if(y&1) res=1ll*res*x%mod;
        x=1ll*x*x%mod;
        y/=2;
    }
    return res;
}
int main(){
    n=read();
    for(int i=1;i<=n;i++) a[i]=read();//分母 
    for(int i=1;i<=n;i++) b[i]=read(),b[i]=a[i]-b[i];//分子 
    ans=1;
    for(int i=1;i<=n;i++) ans=1ll*ans*b[i]%mod*ksm(a[i],mod-2)%mod;
    printf("%d\n",(mod+1-ans)%mod);//确保非负 
    return 0;
}