题目描述
求关于 x 的同余方程 ax ≡ 1 (mod b)的最小正整数解。
输入格式:
输入只有一行,包含两个正整数 a, b,用一个空格隔开。
输出格式:
输出只有一行,包含一个正整数 x0,即最小正整数解。输入数据保证一定有解。
输入样例#1:
3 10
输出样例#1:
7
题解
对于
ax≡1(modb)
a
x
≡
1
(
mod
b
)
可以转化为:
ax−kb=1
a
x
−
k
b
=
1
,我们可以直接用
ExGcd()
E
x
G
c
d
(
)
,来求出一组解,然后套用同解公式不断的变小就可以。
注意开
longlong
l
o
n
g
l
o
n
g
不会扩展欧几里得的戳
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<cmath>
#define ll long long
using namespace std;
ll a,b;
ll Gcd(ll a,ll b)
{
if(!b) return a;
else return Gcd(b,a%b);
}
void Ex_Gcd(ll a,ll b,ll &x,ll &y)
{
if(!b)
{
x=1;y=0;
return ;
}
Ex_Gcd(b,a%b,x,y);
int temp=x;
x=y;
y=temp-(a/b)*y;
return ;
}
int main()
{
scanf("%lld%lld",&a,&b);
ll gcd=Gcd(a,b);
ll x,y;
Ex_Gcd(a,b,x,y);
ll lcm=a*b/gcd;//通解
while(x<0) x+=lcm/a;
while(x-(lcm/a)>0) x-=lcm/a;
printf("%lld",x);
return 0;
}