同余方程(扩欧)

P1082 [NOIP2012 提高组] 同余方程

求关于 x 的同余方程 ax≡1(modb) 的最小正整数解。

可以将式子转化成ax + by = 1

那么可以用扩展欧几里得求得x和y的值,ax' + by' = d = gcd(a, b)

易得x = x' / d;

那么如何保证x为最小正整数解呢?

因为b是模数,所以x每次加b取模b即可

#include<bits/stdc++.h>

using namespace std;

int x, y;

int exgcd(int a, int b, int &x, int &y) {
	if(b == 0) {
		x = 1, y = 0;
		return a;
	}
	
	int d = exgcd(b, a % b, y, x);
	y -= a / b * x;
	
	return d;
}

int main() {
	int a, b;
	cin >> a >> b;
	
	int d = exgcd(a, b, x, y);
	while(x < 0) x += b;
	
	cout << x % b;
	
	return 0;
}

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值