洛谷P1082 同余方程

链接:P1082

题目描述
求关于xx的同余方程 a x \equiv 1 \pmod {b}ax≡1(modb) 的最小正整数解。

输入格式
一行,包含两个正整数 a,ba,b,用一个空格隔开。

输出格式
一个正整数 x,即最小正整数解。输入数据保证一定有解。

输入输出样例
输入 #1
3 10
输出 #1
7
说明/提示
【数据范围】

对于 40%的数据,2 ≤b≤ 1,0002≤b≤1,000;

对于 60%的数据,2 ≤b≤ 50,000,0002≤b≤50,000,000;

对于 100%的数据,2 ≤a, b≤ 2,000,000,0002≤a,b≤2,000,000,000。

NOIP 2012 提高组 第二天 第一题


第一道数论?
对于同余方程ax ≡ 1 (mod b),如果转化为我们易懂的语言就是
求满足ax%b=1,1%b=1最小正整数解。

扩展欧几里得: 通常谈到最大公约数时,我们都会提到一个非常基本的事实:给予二个整数a、b,必存在整数x、y使得ax + by =
gcd(a,b)。
有两个数a,b,对它们进行辗转相除法,可得它们的最大公约数——这是众所周知的。然后,收集辗转相除法中产生的式子,倒回去,可以得到ax+by=gcd(a,b)的整数解。

但是我们在这里可以发现
对于任意a,b,x和y存在很多解
但是我们利用辗转相除法可以得到这个扩展欧几里得的方程的最小整数解·
而这个最小整数解就是这个答案所要求的

我们可以尝试解一下这个方程

ax+by=gcd(a,b)
根据辗转相除法
则得到ax+by=gcd(b,a%b);=bx1+(a%b)y1
ax+by=bx1+(a%b)y1
a(x-y1)=b[x1-y-(a/b)y1]
由此可以解得
x=y1,y=x1-[a/b]y1;

因此我们可以使用递归的方式进行解得辗转时当前的x1,x2
而当辗转到gcd(b,0),x=1,y=0,就可以从上进行调用了
则得到了最小解

但是最小正整数解一定是大于零的,这里这个方程可以能会解出负值
所以我们需要来把最后输出值给改为(x+q)%q (这里的q为输入的模值)

代码:

#include <cstdio>
#include <iostream>
using namespace std;
int a,q;
int x,y;
void exgcd(int a,int b)
{
	if (b==0)
	{
		x=1;
		y=0;
		return ;  //得到gcd(b,0)时到达边界值
	}  //
	else
	{
		exgcd(b,a%b);
		int k=x;
		x=y;
		y=k-(a/b)*y;  //根据上方推出的公式进行递归求出结果
	}
	return ;
}
int main()
{
	scanf("%d%d",&a,&q);  
	exgcd(a,q);
	printf("%lld",(x+q)%q);
	return 0;
}
  • 0
    点赞
  • 4
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值