定义:设a , b为常数,对于一个表达式 a \* x + b \* y =GCD( a , b ) ,一定存在解( x , y)使之成立。
我们就可以通过扩展原来的辗转相除法来求解。
求解的过程如下(**最好自己手推一下**):
> 初始表达式
> `$ a * x + b * y =GCD( a , b ) $`
> 由之前的知识可得
> `$ GCD( a , b ) =GCD( b , a \% b ) $`
> 因此
> `$ a*x_1+b*y_1=GCD( a , b ) =GCD( b , a\%b ) =b*x_2+(a\%b)*y_2 $`
> `$ a\%b=a-[a/b]*b $`
> 其中\[a/b\]表示整除。带入化简可得:
> `$ a*x_1+b*y_1=a*x_2+b*(x_2-[a/b])*y_2 $`
> 由恒等关系可得
> `$ x_1=y_2 $`
> `$ y_1=(x_2-[a/b]) $`
> 因此我们只要求出 x2 和 y2 的值就可以求解 x1 和 y1。而 x2 , y2 可通过同种方法求解。
>
> 特别的,当 b=0时,表达式为 a \* x + b \* y = GCD( a , b ) =GCD( a , 0 ) = 0
>
> 此时,可求得 x = 1 。y 的值对表达式的值没有影响。
>
> 上面就是求解初始表达式的方法。
代码如下:
~~~
int extGcd(int a,int b,int &x,int &y)//这里取 x , y的地址
{
if(b==0)
{
x=1;y=0;
return a;//函数的返回值一直是 gcd(a,b)
}
int r=exGcd(b,a%b,x,y);
int t=x;
x=y; //x1=y2
y=t-a/b*y; //y1=x2-[a/b]*b
return r;
}
~~~
- 简介
- 零基础快速入门ACM
- C语言快速入门
- C语言快速入门
- C/C++基础
- 输入输出
- 数组和字符串
- 数组
- 字符数组
- 函数和递归
- C++标准容器库(STL)
- Vector
- Map
- Set
- String
- Stack
- Queue
- Priority_queue
- 贪心
- 硬币问题
- 区间调度问题
- 字典序最小问题
- 独木舟问题
- 搜索
- DFS
- BFS
- 剪枝
- 记忆化搜索
- 常用技巧
- 倍增
- 高精度
- 大数加法
- 大数减法
- 大数乘法
- 大数除法
- 大数阶乘
- 输入挂
- 前缀和
- 二分
- 本地预处理
- 尺取
- 枚举
- 坐标离散化
- 分治
- 数列分治
- 树上分治
- 平面分治
- 计算几何
- 凸包
- 点积
- 叉积
- 线段关系
- 皮克定理
- 最近点对
- 数据结构
- 线段树
- 树状数组
- 并查集
- 动态规划
- 基础知识
- 信息学竞赛中的动态规划
- 动态规划初步
- 最长上升子序列(LIS)
- 最长公共子序列(LCS)
- 最大子段和
- 背包问题
- 部分背包
- 0 1 背包
- 完全背包
- 多重背包
- 背包的可行性问题
- 线性DP
- 树形DP
- 区间DP
- 数位DP
- 动态规划优化
- 字符串
- KMP算法
- 拓展KMP
- 字符串Hash
- Manacher算法
- 后缀数组
- Trie树
- 博弈论
- Nim博弈
- Bash博弈
- 斐波那契博弈
- 威佐夫博弈
- SG函数
- 数论
- 整数研究
- 素数判断
- 素数筛选
- 快速幂
- 算数基本定理
- 最大公约数(Gcd)
- 费马小定理
- 扩展欧几里得
- 逆元
- 同余定理
- 组合数学
- 容斥原理
- 抽屉原理
- 卢卡斯(Lucas)
- 卡特兰(Catalan)
- 著名函数
- 欧拉函数
- 莫比乌斯函数
- 线性代数
- 矩阵运算
- 矩阵快速幂
- 图论
- 存图方法
- 邻接矩阵
- 邻接表
- Vector存图
- 链式前向星
- 图的遍历
- 拓扑排序
- 最短路
- Dijkstra算法
- SPFA算法
- Floyed 算法
- 最小生成树
- Prim算法
- Kruskal算法
- 最近公共祖先 (LCA)
- 二分图匹配
- 网络流
- 其他
- 莫队