先来介绍一下欧几里得的应用
欧几里得定理主要用于求两个数的最大公约数
核心等式 gcd(a,b)=gcd(b%a,a) 前提是a不等于0
证明:
证明两个正数大小相等可以转换为证明这两个整数可以相互整除
先设d=gcd(a,b) 现在证明d | gcd(b%a,a)
b%a=b-[b/a]*a,由于b%a是a与b的线性组合,故有d | (b%a)又因为d | a,所以gcd(a,b) | gcd(b%a,a)成立
下证gcd(b%a,a) | gcd(a,b)成立
设gcd(b%a,a)=d, 则有d | (b%a)以及d | a;
又因为b%a=b-[b/a]*a,所以 d | b-[b/a]*a+[b/a]*a,也就是说d | b,所以有gcd(b%a,a) | gcd(a,b)成立
得证!
给定 n 对正整数 ai,bi,请你求出每对数的最大公约数。
输入格式
第一行包含整数 n。
接下来 n 行,每行包含一个整数对 ai,bi。
输出格式
输出共 n 行,每行输出一个整数对的最大公约数。
数据范围
1≤n≤10^5
1≤ai,bi≤2×10^9
输入样例:
2
3 6
4 6
输出样例:
3
2
#include<bits/stdc++.h>
using namespace std;
int gcd(int x,int y)
{
return x==0?y:gcd(y%x,x);
}
int main()
{
int n;
cin>>n;
int a,b;
while(n--)
{
cin>>a>>b;
printf("%d\n",gcd(a,b));
}
return 0;
}
下面来介绍一下扩展欧几里得及其应用
先引入一下裴蜀定理:对于任何一对整数a,b,对于a*x+b*y=gcd(a,b)一定有解
接下来我将通过构造解的方法来证明这个定理
a*x1+b*y1=gcd(a,b) (b%a)*x2+a*y2=gcd(b%a,a)
由欧几里得定理我们可以知道gcd(a,b)=gcd(b%a,a);
则a*x1+b*y1=(b%a)*x2+a*y2①
(b%a)*x2+a*y2 = (b-[b/a]*a)*x2+a*y2 = a*(y2-[b/a]*x2)+b*x2②
①②a与b对应相等,可知 x1=y2-[b/a]*x2 y1=x2
由此我们可以通过递归求出x,y;
给定 n 对正整数 ai,bi,对于每对数,求出一组 xi,yi,使其满足 ai*xi+bi*yi=gcd(ai,bi) 。
输入格式
第一行包含整数 n。
接下来 n 行,每行包含两个整数 ai,bi。
输出格式
输出共 n行,对于每组 ai,bi求出一组满足条件的 xi,yi每组结果占一行。
本题答案不唯一,输出任意满足条件的 xi,yi均可。
数据范围
1≤n≤10^5
1≤ai,bi≤2×10^9
输入样例:
2
4 6
8 18
输出样例:
-1 1
-2 1
#include<bits/stdc++.h>
using namespace std;
int exgcd(int a,int b,int &x,int &y)
{
if(a==0)
{
x=0;y=1;
return b;
}
int d=exgcd(b%a,a,y,x);
x-=b/a*y;
return d;
}
int main()
{
int n,a,b;
cin>>n;
while(n--)
{
cin>>a>>b;
int x,y;
int t=exgcd(a,b,x,y);
printf("%d %d\n",x,y);
}
return 0;
}
下面引入一道欧几里得扩展的变形题目
给定 n 组数据 ai,bi,mi,对于每组数求出一个 xi,使其满足 ai×xi≡bi(modmi),如果无解则输出 impossible
。
输入格式
第一行包含整数 n。
接下来 n 行,每行包含一组数据 ai,bi,mi。
输出格式
输出共 n 行,每组数据输出一个整数表示一个满足条件的 xi,如果无解则输出 impossible
。
每组数据结果占一行,结果可能不唯一,输出任意一个满足条件的结果均可。
输出答案必须在 int 范围之内。
数据范围
1≤n≤10^5,
1≤ai,bi,mi≤2×10^9
输入样例:
2
2 3 6
4 3 5
输出样例:
impossible
-3
/*题目分析:
a*x≡b (mod m) 有解等价于存在y使得a*x-m*y=b
由此题目转化为扩展欧几里得定理的应用
我们可以用扩展欧几里得定理求出a*x-m*y=gcd(a,m)的系数
也可以看出方程有解的条件是gcd(a,m)|b
最后别忘了将解乘以gcd(a,m)|b
*/
#include<bits/stdc++.h>
using namespace std;
int exgcd(int a,int b,int &x,int &y)
{
if(a==0)
{
x=0;y=1;
return b;
}
int t=exgcd(b%a,a,x,y);
int tx=x,ty=y;
x=ty-b/a*tx;y=tx;
return t;
}
int main()
{
int n;
cin>>n;
int a,b,x,y,m;
while(n--)
{
cin>>a>>b>>m;
int t=exgcd(a,m,x,y);
if(b/t*t==b) printf("%lld\n",(long long)x*(b/t)%m);//求得是关于gcd(a,b)的等式,所以最后要乘以 b/t
else printf("impossible\n");
}
return 0;
}
拓展:
若关于二元方程 AX+BY=gcd(A,B)已经通过拓展欧几里得找到一组解X0,Y0,则可得到二元方程AX+BY=C的一组解X1=X0*(C/gcd(A,B)),Y1=Y0*(C/gcd(A,B)),则方程AX+BY=C的通解为X=X1+n*(B/gcd(A,B)),Y=Y1-n*(A/gcd(A,B));
若求一个数x在模n下的最小解,则可利用取余的方法来求:x=(x%n+n)%n;(前提n是正数)