扩展欧几里得算法-线性同余方程

扩展欧几里得算法-线性同余方程

概念：给定a, b, m，求整数x,使得： ax = b(mod m)；

举例：2x = 3 (mod 6) 无解

          4x = 3 (mod 5) :4*某一个数 同于于 4*某一个数/5的余数是3. x可以等于 2，7等。

等式变形：

ax = b(mod m) 等价于 存在一个整数y属于z，使得 ax = m*y + b  --> ax - my = b --> y` = y, ax + my` = b;

注意：该式子能成立的条件是b能整除(a,m)最大公约数，(a, m) | b, 那么就有解，否则一定无解。

Example

给定nn组数据ai,bi,miai,bi,mi，对于每组数求出一个xixi，使其满足ai∗xi≡bi(mod mi)ai∗xi≡bi(mod mi)，如果无解则输出impossible。

输入格式

第一行包含整数nn。

接下来nn行，每行包含一组数据ai,bi,miai,bi,mi。

输出格式

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

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

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

数据范围

1≤n≤1051≤n≤105,
1≤ai,bi,mi≤2∗1091≤ai,bi,mi≤2∗109

输入样例：

2
2 3 6
4 3 5

//扩展欧几里得算法
#include <iostream>
 
using namespace std;

typedef long long LL; 

int exgcd(int a, int b, int &x, int &y)
{
    //return b ? gcd(b, b % a) : a;
    if (!b) //如果b = 0,ax + by = a
    {
        //得到一组解
        x = 1, y = 0;
        return a;//如果b = 0,返回a
    }
    //by + (a mod b)x = (a, b)最大公约数d
    int d = exgcd(b, a % b, y, x);//存下最大公约数,翻转了一下
    //a mod b = a - (a / b)下取整 * b带入
    //by + (a - (a / b)下取整 * b)x = d;
    //展开整理：ax + b (y - (a/b)*x)
    //得到a的系数不变为a,y的系数变成：y - (a/b)*x
    y -= a / b * x;
    
    return d;
    
}
 
int main()
{
    int n;
    scanf("%d", &n);
    
    while ( n --)
    {
        int a, b, m;
        scanf ("%d%d%d", &a, &b, &m);
        int x, y;
        int d = exgcd(a, m, x, y);//存下结果
        //使用x的系数
        if (b % d) puts("impossible");//如果b不是d的倍数，输出无解
        else printf("%d\n", (LL)x * (b / d) % m);//否则输出x扩大若干倍，mod到m范围内
    }
    return 0;
}