扩展欧几里得算法

扩展欧几里得算法

数学推导

现有一二元一次方程：
$$a{x}_1+b{y}_1=m$$
$$a{x}_2+b{y}_2=gcd(a,b)$$
由数学中的某一定理可知（裴蜀定理，直接用结论），当二元一次方程的结果$m$是两系数的最大公因数或是最大公因数的倍数，则我们可以得至少有一组整数解可以满足方程。

且$x_1\%x_2==y_1\%y_2==0,x_1/x_2==y_1/y_2$。

现有一算法通过gcd的性质扩展可以求得一组$x_2,y_2$。

int exgcd(int a ,int b ,int &x , int &y){
    if(!b){
        x = 1 , y = 0 ;
        return a ; 
    }
    else {
        int d = exgcd(b , a%b , y , x) ;
        y = y - (a/b) * x ;
        return d ;
    }
}

同余方程（ACwing878）

原题

题目描述

给定 n 组数据 ai,bi,mi，对于每组数求出一个$x_i$，使其满足 $a_i\times x_i\equiv b_i(mod{m}_i)$ ，如果无解则输出 impossible。

输入格式

第一行包含整数 $n$。

接下来 $n$ 行，每行包含一组数据 $a_i,b_i,m_i$。

输出格式

输出共 $n$ 行，每组数据输出一个整数表示一个满足条件的 $x_i$，如果无解则输出 impossible。

每组数据结果占一行，结果可能不唯一，输出任意一个满足条件的结果均可。

输出答案必须在 int 范围之内。

数据范围

$1\le n\le 10^5,$
$1\le a_i,b_i\le 2\times 10^9$

样例

2
2 3 6
4 3 5

impossible
-3

思路

$$a_i\times x_i\equiv b_i(mod{m}_i)$$
可以转化为
$$a_i\times x_i-m_i\times y_i=b_i$$
令$y_i=-y_i$
得
$$a_i\times x_i+m_i\times y_i=b_i$$
之后便可利用扩展欧几里得解

#include <iostream>
using namespace std ;

typedef long long ll ;

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 ;
}


int main(){
    int n ; 
    cin >> n ;
    while(n -- ){
        int a ,b , m ;
        cin >> a >> b >> m ;
        int x ,y ;
        int d = exgcd(a , m , x , y ) ;
        
        if(b%d) puts("impossible");
        else cout <<  ((ll)x * b / d % m + m) % m << endl ;
        
    }
    return 0 ;
}