链接: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;
}