同余方程 (扩展欧几里得) (洛P1082)

题目地址

题目化简 a x − y b = 1 ax -yb = 1 axyb=1 要求最小x.典型扩展欧几里得方程式。

显然 1 ∣ g c d ( a , b ) 1 | gcd(a,b) 1gcd(a,b) 故方程始终有解。

用模板求出 x.

对方程 a x + b y = c , x = x + k b 且 b = b g c d ( a , b ) ax+by = c,x = x + kb 且 b= \frac{b}{gcd(a,b)} ax+by=cx=x+kbb=gcd(a,b)b

博客,性质1

所以把答案范围控制再 (0, b g c d ( a , b ) \frac{b}{gcd(a,b)} gcd(a,b)b-1)

代码:

#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#include<bits/stdc++.h>
#define int long long
using namespace std;
typedef pair<int,int> pii;
typedef long long ll;

const int INF = 0x3f3f3f3f;
const double eps = 1e-5;
const int mod = 1e9+7;
const int N = 2e5+10;

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;
}

signed main(){
	IOS
	#ifdef ddgo
		freopen("C:\\Users\\asus\\Desktop\\ddgoin.txt","r",stdin);
    #endif
	int a,b; cin>>a>>b;
	int x,y;
	int d = exgcd(a,b,x,y);
    int p = b/d;
	cout<<(x % p + p) % p<<endl;
	return 0;
}

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

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值