（AcWing）线性同余方程

最新推荐文章于 2024-09-03 20:45:39 发布

Go Free919

最新推荐文章于 2024-09-03 20:45:39 发布

阅读量58

点赞数

文章标签：算法图论

本文链接：https://blog.csdn.net/GF0919/article/details/131792198

版权

文章介绍了一种利用扩展欧几里得算法求解线性同余方程ax≡b(modm)的方法。通过计算最大公约数d并检查b是否能被d整除来判断方程是否有解。如果有解，则可以通过算法得出x的值，并确保结果在int范围内。

摘要由CSDN通过智能技术生成

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

输入格式

第一行包含整数 n。

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

输出格式

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

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

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

数据范围

1≤n≤105,
1≤ai,bi,mi≤2×109

输入样例：

2
2 3 6
4 3 5

输出样例：

impossible
-3

/*题目：a*x%m = b，
则可化简为：a*x = m*y+b;
a*x+m*(-y) = b  (扩展欧几里得)
用扩展欧几里得求出x的值
*/

#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 t;
    cin>>t;
    while(t--)
    {
        int a,b,m,x,y;
        cin>>a>>b>>m;
        int d = exgcd(a,m,x,y);
        if(b%d) cout<<"impossible"<<endl;
        else cout<<(ll)b/d*x%m<<endl;
        /*考虑线性同余方程 ax ≡ b (mod m)，其中 a、b 和 m 都是整数。根据模运算的定义，两个整数 x 和 y 满足同余关系（即 x ≡ y (mod m)）当且仅当它们除以 m 的余数相等。
        求出的最大公约数d一定是b的倍数，d/b是整个a*x+m*(-y) = b扩大的倍数，x应该进行相应的扩大，%m是保证其在答案在int范围*/
    }
}