NOIP2012 同余方程的两种解法

最新推荐文章于 2024-08-03 12:01:10 发布

ngncmh

最新推荐文章于 2024-08-03 12:01:10 发布

阅读量7.2k

点赞数

分类专栏：分块数学扩展欧几里得

本文链接：https://blog.csdn.net/KQZXCMH/article/details/8199423

版权

这篇博客介绍了如何求解同余方程ax mod b = 1的最小正整数解，重点讨论了扩展欧几里得算法和利用预处理进行快速查找的方法。通过将数值区间划分并利用已知的倍数块内的模运算结果，可以有效地找到满足条件的最小正整数解。

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

求 ax mod b = c 的最小正整数x解题目中的c=1

解法一：经典的扩展欧几里得算法解ax+by=1的一组特解，然后算出最小正整数解

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>

using namespace std;
long long a,b,x,y;
long long extended_gcd(long long  a,long long b,long long &x,long long &y)
{
    if (!b) { x=1; y=0; return a; }
    else { 
         long long ans=extended_gcd(b,a%b,x,y); 
         long long t=x; x=y; y=t-(a/b)*y;
         return ans;
    }
}
int main()
{
    cin>>a>>b;
    long long ans = extended_gcd(a, b, x, y);
    //cout << a*x+b*y<<endl;

    if (ans!=1) cout<<"No answer"<<endl;
    else
    {
        x=x%b;
        while (x<0) x+=b;
        cout <<x<<endl;
    }
    system("pause");
    return 0;
}