求关于 x 的同余方程 ax ≡1(modb)的最小正整数解。
输入格式
输入只有一行,包含两个正整数 a,b,用一个空格隔开。
输出格式
输出只有一行,包含一个正整数 x0,即最小正整数解。输入数据保证一定有解。
数据范围
对于 40% 的数据,2≤b≤1,000;
对于 60% 的数据,2≤b≤50,000,000;
对于 100% 的数据,2≤a,b≤2,000,000,000。
Sample Input
3 10
Sample Output
7
---------------------------------------------------------------------------------------------------------------------------------
本质就是求a的逆元,这里使用拓展欧几里得方法求解,因为费马小定理要求b是质数,输入的b不一定符合要求
#include<stdio.h>
#include<iostream>
#include<string.h>
#include<math.h>
#define ll long long
using namespace std;
ll exgcd(ll a,ll b,ll &x,ll &y)
{
if(b==0)
{
x=1,y=0;
return a;
}
ll d=exgcd(b,a%b,x,y);
ll t=x;
x=y;
y=t-a/b*y;
return d;
}
int main()
{
ll a,b,x,y,d;
cin>>a>>b;
d=exgcd(a,b,x,y);
x=(x%b+b)%b;
cout<<x;
}