noip2012 同余方程 exgcd

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

Loi_a

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

阅读量658

点赞数 2

分类专栏： ===数学=== 文章标签：扩展欧几里得

本文链接：https://blog.csdn.net/Loi_a/article/details/53101063

版权

===数学=== 专栏收录该内容

20 篇文章 0 订阅

订阅专栏

exgcd是求ax1+by1=gcd(a,b)的一组解。

ax1+by1 = gcd (a,b)
bx2+(a%b) y2 = gcd (b,a%b)

ax1+by1 = bx2+ (a%b) y2
∵ a%b 可以写成 a-[a/b]*b

ax1+by1=ay2+b ( x2-[a/b]*y2 )

x1 = y2 ; y1 = x2-[a/b]*y2 ;

这里求出的x可以为负数，所以要改为正数。

#include<iostream>
#include<cstdio>
using namespace std;
void exgcd(int a,int b,int &x,int &y)
{
    if(!b){x=1;y=0;return;}
    exgcd(b,a%b,x,y);
    int t=x;x=y;y=t-a/b*y;
}
int main()
{
    int a,b,x,y;
    scanf("%d%d",&a,&b);
    exgcd(a,b,x,y);
    x=(x%b+b)%b;
    cout<<x;
}