exgcd

exgcd学习笔记

exgcd的应用

exgcd即扩展欧几里得算法，用来解决形如 $ax+by=gcd(a,b)$ ($a,b$ 为常数) 的方程的一组整数解 。

exgcd推导过程

当 $b=0$ 时，容易发现方程的解即为 $x=1,y=0$

当 $b\ne 0$ 时，考虑到 $gcd(a,b)=gcd(b,amodb)$ ，假设函数 $exgcd(a,b,x,y)$ 中的 $x,y$为解

那么 $exgcd(a,b,x,y)=exgcd(b,amodb,x_0,y_0)$

解这个方程

$ax+by=b{x}_0+(amodb)y_0$

而 $amodb=a-\lfloor \frac{a}{b}\rfloor b$

所以 $ax+by=b{x}_0+a{y}_0-\lfloor \frac{a}{b}\rfloor b{y}_0$

即 $ax+by=a{y}_0+b(x_0-\lfloor \frac{a}{b}\rfloor y_0)$

则方程 $ax+by=gcd(a,b)$ 的一个解为 $x=y_0,y=x_0-\lfloor \frac{a}{b}\rfloor y_0$

通俗一点的代码

int exgcd(int a,int b,int &x,int &y){
    if(!b){
        x=1;y=0;
        return a;
    }
    int x0,y0;
    int d=exgcd(b,a%b,x0,y0);//先求解x0,y0
    x=y0,y=x0-a/b*y0;//给目前的解赋值
    return d;
}

简洁的代码

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);//先求解gcd(b,a%b)的解
    y-=a/b*x;
    return d;
}

返回值为 $gcd(a,b)$

exgcd的通解

假设存在 $x_0=x+i,y_0$ 也满足 $a{x}_0+b{y}_0=gcd(a,b)=ax+by$

那么 $ax+ai+b{y}_0=ax+by$

化简得 $y_0=\frac{by-ai}{b}$

即$y_0=y-\frac{ai}{b}$

需要满足 $y_0$ 为整数 即需找出最小的 $i$ 满足 $a\ast imodb=0$

首先 $i$ 为 $b$ 一定是可行的 考虑怎么缩小 $i$

容易发现 $a\ast i$ 一定会包含 $a$ 的因数 那么我们就可以将 $i$ 中 $a$ 的因数删去

即 $i=b/gcd(a,b)$

那么

$$\begin{gathered} x_0=x+i=x+\frac{b}{gcd(a,b)} \\ {y}_0=y-\frac{ai}{b}=y-\frac{a}{gcd(a,b)} \end{gathered}$$

即$\{\begin{array}{l}{x}_0=x+\frac{b}{gcd(a,b)}\\ y_0=y-\frac{a}{gcd(a,b)}\end{array}$

当然 这样表示不合理 我们用一般用 $x_0,y_0$表示特解 $x,y$ 表示通解。

于是方程 $ax+by=gcd(a,b)$ 的通解即为
$$\{\begin{array}{l}x=x_0+k\frac{b}{gcd(a,b)}\\ y=y_0-k\frac{a}{gcd(a,b)}\end{array},k\in Z$$

exgcd的最小正整数解

令 $d_x$ 代表通解中 $x$ 的变化量，$d_y$ 代表通解中 $y$ 的变化量

$k_x$代表将特解 $x_0$ 变为最小正整数解的 $k$ 值

$k_y$代表将特解 $y_0$ 变为最小正整数解的 $k$ 值

需要注意通解中 $x,y$ 变化量符号是相反的

那么满足 $x,y$ 均为正整数的 $k$ 的区间即为 $[k_x,k_y]$

if(x<=0){
    Kx+=-x/dx+1; //向下取整+1
}
else {
    Kx-=(x-1)/dx; //减到1(因为x>0) 向下取整
}
if(y<=0){
    Ky-=-y/dy+1;
}
else {
    Ky+=(y-1)/dy;
}

看这道题目 P5656

题目描述

给定不定方程

$ax+by=c$

若该方程无整数解，输出 $-1$。
若该方程有整数解，且有正整数解，则输出其正整数解的数量，所有正整数解中 $x$ 的最小值，所有正整数解中 $y$ 的最小值，所有正整数解中 $x$ 的最大值，以及所有正整数解中 $y$ 的最大值。
若方程有整数解，但没有正整数解，你需要输出所有整数解中 $x$ 的最小正整数值， $y$ 的最小正整数值。

正整数解即为 $x,y$ 均为正整数的解，$0$ 不是正整数。
整数解即为 $x,y$ 均为整数的解。
$x$ 的最小正整数值即所有 $x$ 为正整数的整数解中 $x$ 的最小值，$y$ 同理。

输入格式

第一行一个正整数 $T$，代表数据组数。

接下来 $T$ 行，每行三个由空格隔开的正整数 $a,b,c$。

输出格式

$T$ 行。

若该行对应的询问无整数解，一个数字 $-1$。
若该行对应的询问有整数解但无正整数解，包含 $2$ 个由空格隔开的数字，依次代表整数解中，$x$ 的最小正整数值，$y$ 的最小正整数值。
否则包含 $5$ 个由空格隔开的数字，依次代表正整数解的数量，正整数解中，$x$ 的最小值，$y$ 的最小值，$x$ 的最大值，$y$ 的最大值。

代码(开long long)

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<iomanip>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<stack>
//#include<random>
//#include<bits/stdc++.h>
using namespace std;
#define int long long
#define pii pair<int,int>
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;
}
void solve(){
    int a,b,c;
    cin>>a>>b>>c;
    int g=__gcd(a,b);
    if(c%g){
        cout<<-1<<'\n';
        return ;
    }
    int x,y;
    exgcd(a,b,x,y);
    x*=c/g,y*=c/g;
    int dx=b/g,dy=a/g;
    int Kx=0,Ky=0;
    if(x<=0){
        Kx+=-x/dx+1;
    }
    else {
        Kx-=(x-1)/dx;
    }
    if(y<=0){
        Ky-=-y/dy+1;
    }
    else {
        Ky+=(y-1)/dy;
    }
    if(Ky-Kx+1<=0){
        cout<<x+Kx*dx<<' '<<y-Ky*dy<<'\n';
    }
    else {
        cout<<Ky-Kx+1<<' ';
        cout<<x+Kx*dx<<' '<<y-Ky*dy<<' ';
        cout<<x+Ky*dx<<' '<<y-Kx*dy<<'\n';
    }
}
signed main(){
    ios::sync_with_stdio(false);
    cin.tie(0),cout.tie(0);
    int t=1;
    cin>>t;
    while(t--)solve();
    return 0;
}