数学知识：线性同余方程—扩展欧几里得算法

最新推荐文章于 2024-06-04 11:16:52 发布

奋斗吧！骚年！

最新推荐文章于 2024-06-04 11:16:52 发布

阅读量84

点赞数

分类专栏： ----------AcWing算法基础---------- # 数学知识文章标签：算法 c++ 数据结构扩展欧几里得

本文链接：https://blog.csdn.net/qq_46470984/article/details/125436833

版权

----------AcWing算法基础---------- 同时被 2 个专栏收录

95 篇文章 8 订阅

订阅专栏

数学知识

19 篇文章 0 订阅

订阅专栏

该博客介绍了如何使用欧几里得算法求解线性同余方程，并给出了C++代码实现。文章通过具体示例展示了如何在给定一组数据ai, bi, mi时找到满足条件的xi，当无解时输出'Impossible'。内容涉及算法原理和编程实践，适合计算机科学和技术领域的读者学习。

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

题目： AcWing 878. 线性同余方程
给定 n 组数据 a_i,b_i,m_i，对于每组数求出一个 x_i，使其满足 a_i×x_i≡b_i(mod m_i)，如果无解则输出 impossible。

输入格式
第一行包含整数 n。

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

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

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

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

数据范围
1≤n≤10⁵,
1≤a_i,b_i,m_i≤2×10⁹

输入样例：

2
2 3 6
4 3 5

输出样例：

impossible
-3

#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;
    int a,b,m;
    while(n--)
    {
        cin>>a>>b>>m;
        int x,y;
        int d=exgcd(a,m,x,y);
        if(b%d)cout<<"impossible"<<endl;
        else cout<<(LL)b/d*x%m<<endl;
    }
    return 0;
}