exgcd——同余方程
题目描述
求关于 x 的同余方程 ax ≡ 1 (mod b) 的最小正整数解。
输入描述
输入只有一行,包含两个正整数a,b,用一个空格隔开。
输出描述
输出只有一行,包含一个正整数 x0,即最小正整数解。输入数据保证一定有解。
示例
输入
3 10
输出
7
备注
对于 40% 的数据,有 2 ≤ b ≤ 1000;
对于 60% 的数据,有 2 ≤ b ≤ 50000000;
对于 100% 的数据,有 2 ≤ a,b ≤ 2000000000。
分析
ax ≡ 1 (mod b) 可以转化成 ax + by = 1
这样就可以运用 exgcd 来求解,但所求的解可能为负数,所以应对 x 加上 b 再取余,即 (x + b) % b
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
void exgcd(ll a,ll b,ll &x,ll &y)
{
if(b==0){
x=1;
y=0;
return ;
}
exgcd(b,a%b,x,y);
ll k=y;
y=x-a/b*y;
x=k;
return ;
}
int main()
{
ll a,b,x,y;
cin>>a>>b;
exgcd(a,b,x,y);
cout<<(x+b)%b<<endl;
return 0;
}