欧几里得定理及其扩展应用

先来介绍一下欧几里得的应用

可用于求两个数的最大公约数

核心等式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)成立

　　得证！

 

给定 nn 对正整数 ai,biai,bi，请你求出每对数的最大公约数。

输入格式

第一行包含整数 nn。

接下来 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;

 

给定 nn 对正整数 ai,biai,bi，对于每对数，求出一组 xi,yixi,yi，使其满足 ai×xi+bi×yi=gcd(ai,bi)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;
}

希望大家能够喜欢！

标签：b%,gcd,int,欧几里得,扩展,bi,ai,x2,定理
来源： https://www.cnblogs.com/AC--Dream/p/14730433.html