同余方程

求关于 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;
} 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值