数学知识：扩展欧几里得算法

本文档基于 AcWing 制作

什么是欧几里得算法？

定义

欧几里得算法$(gcd)$，又叫辗转相除法，是指用于计算两个非负整数a，b的最大公约数。

证明：$gcd(a,b)=gcd(b,amodb)$

证明:
设$a\le b,a=b\ast r+c,d=(a,b)$,则$c=a\%b=a-b\ast r$
(1).$gcd(a,b)<=gcd(b,amodb)$
令$d|a$, $d|b$, 则$d|c$, 即$a,b$的公约数同时也是$a,c$的公约数,
$a,b$的最大公约数也是$a,c$的公约数, $gcd(a,b)\le gcd(a,amodb)$
(2).$gcd(b,a)\le gcd(b,amodb)$
令$d|b$, $d|c$,则$d|a$ ,即$b,c$的公约数同时也是$a,b$的公约数, $b,c$的最大公约数同时也是$a,b$的公约数, $gcd(b,amodb)\le gcd(a,b)$
(3).综上可得,gcd(a,b)=gcd(a,a mod b),证毕

代码----欧几里得算法

int gcd(int a,int b)
{
   
    return b?gcd(b,a%b):a;
}

实例：

AcWing 872.最大公约数

原题链接：AcWing 872.最大公约数

题目描述

给定$n$对正整数 $a_i$,$b_i$，请你求出每对数的最大公约数。

输出格式

第一行包含整数$n$。
接下来$n$行,每行包含一个整数对$a_i,b_i$。

输出格式

输出共$n$行，每行输出一个整数对的最大公约数。

数据范围

$1$ $\le$ $n$ $\le$ $10^5$
$1$ $\le$ $a_i,b_i\le 2\times 10^9$

输入样例：

2
3 6
4 6

输出样例：

3
2

代码

• 时间复杂度为$O(\log n)$

#include<bits/stdc++.h>
#define int long long
using namespace std;

int gcd(int a,int b)
{
   
    return b?gcd(b,a%b):a;
}

signed main()
{
   
    ios::sync_with_stdio(false);
    cin.tie(0);cout.tie(0);
    int T;
    cin>>T;
    while(T--)
    {
   
        int a,b;
        cin>>a>>b;
        cout<<gcd(a,b)<<'\n';
    }
    return 0;
}

扩展欧几里得算法

算法基础——裴蜀定理

定义：

定义:若$a,b$是整数，且$gcd(a,b)=d$，那么对于任意整数$x,y$，$ax$