数论 - 裴属定理 + 扩展欧几里得算法 + 线性同余方程

数论 - 扩展欧几里得算法 + 线性同余方程

---

1、裴蜀定理

$对于任意正整数a,b，必有非零整数x,y，使得$$$ax+by=d$$

$其中d=gcd(a,b)。显然也有a{x}^{\prime }+b{y}^{\prime }=k\times d，即x^{\prime }=k\times x,y^{\prime }=k\times y，k\in Z。$

$通俗地讲，ax+by能够取到a和b的最大公约数的整数倍。$

重要推论：

$若a,b互质，即gcd(a,b)=1，则ax+by可取到任意正整数。$

裴属定理的通解:

$设(x_0,y_0)是方程ax+by=d的一组解，则通解为\{\begin{array}{l}x=x_0+k{b}^{\prime }\\ y=y_0-k{a}^{\prime }\end{array}，其中k\in Z，d=gcd(a,b)，a^{\prime }=\frac{a}{d}，b^{\prime }=\frac{b}{d}。$
$$\begin{gathered} 证明： \\ \{\begin{array}{l}a{x}_0+b{y}_0=d\\ ax+by=d\end{array},两式相减，得a(x-x_0)=b(y_0-y),则a^{\prime }(x-x_0)=b^{\prime }(y_0-y)，从而b^{\prime }|a^{\prime }(x-x_0)。 \end{gathered}$$

$又因为gcd(a^{\prime },b^{\prime })=1，所以b^{\prime }|(x-x_0)，所以x=x_0+k{b}^{\prime }。同理，y=y_0-k{a}^{\prime }。k=1,2,...。$

$所以知道一组特殊解，可以求得任意解。$

$$\begin{gathered} 事实上，就是给方程a{x}_0+b{y}_0=d添项。 \\ a{x}_0+b{y}_0=a(x_0+k{b}^{\prime })+b(y_0-k{a}^{\prime })=d，容易验证这个通解。 \end{gathered}$$

---

2、扩展欧几里得算法

$给定n对正整数a_i,b_i，对于每对数，求出一组x_i,y_i，使其满足$$$a_i{x}_i+b_i{y}_i=gcd(a_i,b_i)。$$

输入格式
第一行包含整数n。

接下来n行，每行包含两个整数ai,bi。

输出格式
输出共n行，对于每组ai,bi，求出一组满足条件的xi,yi，每组结果占一行。

本题答案不唯一，输出任意满足条件的xi,yi均可。

数据范围
1≤n≤10⁵,
1≤ai,bi≤2∗10⁹

输入样例：
2
4 6
8 18
输出样例：
-1 1
-2 1

---

分析：

$由欧几里得算法(辗转相除法)，我们知道gcd(a,b)=gcd(b,a\%b)$

$那么方程ax+by=gcd(a,b)<=>by+(a\%b)x=gcd(b,a\%b)。$

$$\begin{gathered} 即by+(a-\lfloor \frac{a}{b}\rfloor \cdot b)x=gcd(b,a\%b)， \\ 整理可得：ax+b(y-\lfloor \frac{a}{b}\rfloor \cdot x)=gcd(b,a\%b)。 \end{gathered}$$

$那么方程ax+by=gcd(a,b)的一组解(x,y)就可以通过先递归，再反推来确定。$

代码：

#include<iostream>
#include<cstdio>

using namespace std;

const int N=1e5+10;

int n;

int exgcd(int a,int b,int &x,int &y)
{
    if(!b)
    {
        x=1,y=0;
        return a;
    }
    
    int d=exgcd(b,a%b,y,x);
    y-=a/b*x;
    
    return d;
}

int main()
{
    cin>>n;
    while(n--)
    {
        int a,b,x,y;
        scanf("%d%d",&a,&b);
        exgcd(a,b,x,y);
        printf("%d %d\n",x,y);
    }
    
    return 0;
}

---

3、线性同余方程

$给定n组数据a_i,b_i,m_i，对于每组数求出一个x_i，使其满足a_i{x}_i\equiv b_i(mod{m}_i)，如果无解则输出impossible。$

输入格式
第一行包含整数n。

接下来n行，每行包含一组数据ai,bi,mi。

输出格式
输出共n行，每组数据输出一个整数表示一个满足条件的xi，如果无解则输出impossible。

每组数据结果占一行，结果可能不唯一，输出任意一个满足条件的结果均可。

输出答案必须在int范围之内。

数据范围
1≤n≤10⁵,
1≤ai,bi,mi≤2∗10⁹

输入样例：
2
2 3 6
4 3 5
输出样例：
impossible
7

---

分析：

$方程ax\equiv b(modm)<=>ax=my+b，其中y\in Z。$

$即求是否有整数解(x_0,y_0)满足a{x}_0+m{y}_0=b。$

$由裴属定理，上述方程有整数解，需保证b=k\times gcd(a,m)，k\in Z。$

$$\begin{gathered} 用扩展欧几里得算法容易解得ax+by=gcd(a,b)的一组解。 \\ 因此我们先求得a{x}_0+m{y}_0=gcd(a,m)的一组解的x_0，再将其扩大\frac{b}{gcd(a,m)}倍即可。 \end{gathered}$$

$$\begin{gathered} 输出答案需在int范围内，根据裴属定理， \\ 我们要将最后的答案对\frac{m}{gcd(a,m)}再取模，对m取模也行。 \end{gathered}$$

注意：$最后做乘法为避免int溢出，要强制类型转换。$

//ax+my=b
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>

using namespace std;

int n, a, b, x, y, m;

int exgcd(int a,int b,int &x,int &y)
{
    if(!b)
    {
        x=1, y=0;
        return a;
    }
    
    int d = exgcd(b,a%b,y,x);
    y-=(a/b)*x;
    return d;
}

int main()
{
    scanf("%d",&n);
    while(n--)
    {
        scanf("%d%d%d",&a,&b,&m);
        int t = exgcd(a,m,x,y);
        if(b%t!=0) puts("impossible");
        else
        {
            int k = b/t;
            x=((long long)x*k)%(m/t);
            printf("%d\n",x);
        }
    }
    
    return 0;
}

---