扩展欧几里德
来源:互联网 发布:南方大数据100指数雪球 编辑:程序博客网 时间:2024/05/24 01:45
http://hihocoder.com/problemset/problem/1297
假设经过t时间过后,速度快的人刚好超过了速度慢的人k圈,且到达同一个位置。
s1+v1t=s2+v2t−km >注:不需要考虑s1
将这个式子进行变换得到:
(v1−v2)t+km=(s2−s1)
令
求解该式子的算法我们称为扩展欧几里德算法。
该算法分为两个部分:
(1) 判定是否存在解
对于形如”Ax+By=C”的式子,其存在解的条件为C为A和B最大公约数的整数倍。
我们将A和B的最大公约数记为gcd(A,B)。因此其有解的条件是C=n*gcd(A,B)。
那么我们应该如何来求解gcd(A,B)呢?
一个朴素的算法是枚举1~min(A,B),最大的一个能同时被A,B整除的数即gcd(A,B)。显然这个算法是非常没有效率的。
为了求解gcd(A,B),欧几里德提出了一个辗转相除法(欧几里得算法):
首先要证明一个定理:gcd(A,B) = gcd(B, A mod B)
证明:
假设A = k*B+r,有r = A mod B。不妨设d为A和B的一个任意一个公约数,则有A = pd, B = qd。
由r = A - k*B = pd - k*qd = (p - kq)*d,所以有d也为r的约数,因此d是B和A mod B的公约数。
由于对任意一个A和B的公约数都满足这个性质,gcd(A,B)也满足,因此有gcd(A,B)=gcd(B,A mod B)。
利用这个性质,我们可以得到算法:
A mod B = 0, 则B为gcd(A,B)
A mod B ≠ 0, 则gcd(A,B) = gcd(B, A mod B)
通过不断的模运算,数据的规模也越来越小,因此能够快速得收敛到一个解。将其写成伪代码为:
gcd(A, B): If (A mod B == 0) Then Return B End If Return gcd(B, A mod B)
(2) 求解
先求解
Ax+By=gcd(A,B)
此时,我们可以先求解出上式的解
(x,y)=(x′∗Cgcd(A,B),y′∗Cgcd(A,B)) 。
那么接下来我们来研究如何求解
假设A>B>0,同时我们设:
A * x[1] + B * y[1] = gcd(A, B)
B * x[2] + (A mod B) * y[2] = gcd(B, A mod B)
已知gcd(A,B)=gcd(B, A mod B),因此有:
A * x[1] + B * y[1] = B * x[2] + (A mod B) * y[2]
A * x[1] + B * y[1] = B * x[2] + (A - kB) * y[2] // A = kB + r
=> A * x[1] + B * y[1] = A * y[2] + B * x[2] - kB * y[2]
=> A * x[1] + B * y[1] = A * y[2] + B * (x[2] - ky[2])
=> x[1] = y[2], y[1] = (x[2] - ky[2])
利用这个性质,我们可以递归的去求解(x,y)。
其终止条件为gcd(A, B)=B,此时对应的(x,y)=(0,1)
或者终止条件为B=0,此时对应的(x,y)=(1,0)
将这个过程写成伪代码为:
extend_gcd(A, B): If (A mod B == 0) Then Return (0, 1) End If (tempX, tempY) = extend_gcd(B, A mod B) x = tempY y = tempX - (A / B) * tempY Return (x, y)
那么我只需要把A=(v1-v2),B=m,C=(s2-s1),x=t,y=k代入就可以得到t了么?
是的,在已知A,B,C的情况下,我们的确能够顺利求解出一组合法的(x,y)。
但是在求解过程中,我们并没有保证x是最小的非负整数,它不能直接作为我们的解。
求最小的非负整数
我们上面求出了AX+BY=C 的一个解
设另一组解为
Ax1+By1=C
Ax2+By2=C
两式相减得
A(x1−x2)=B(y2−y1)
同时除以gcd(A,B),得
A′(x1−x2)=B′(y2−y1)
故,
设A,B,C为任意整数,若方程Ax+By=C的一组整数解为
(x0,y0) ,则它的任意整数解都可以写成(x0+kB′,y0−kA′) 。其中
A′=A/gcd(A,B),B′=B/gcd(A,B) 。
AC代码为:
#include <iostream>#include <cmath>using namespace std;long gcd(long a, long b) { long r; while (b) { r = a % b; a = b; b = r; } return a;}void ExtendGcd(long a, long b, long& x, long& y) { if (b == 0) { x = 1; y = 0; return; } ExtendGcd(b, a % b, x, y); long t = x; x = y; y = t - (a / b) * y;}int main() { long s1, s2, v1, v2, m; cin >> s1 >> s2 >> v1 >> v2 >> m; long x = 0, y = 0; long a = v1 - v2, b = m, c = s2 - s1; long g = gcd(a, b); if (c % g != 0) { cout << "-1" << endl; return 0; } ExtendGcd(a, b, x, y);//ax+by=gcd(a,b) x = x * (c / g);// x/g * c b = abs(b / g); x = (x % b + b) % b; cout << x << endl; return 0;}/*0 1 1 2 6 //ans = 5 1 2 3 4 5 //ans = 4*/
- 欧几里德和扩展欧几里德
- 欧几里德, 与 扩展欧几里德
- 欧几里德|扩展欧几里德
- 欧几里德&&扩展欧几里德详解
- 欧几里德&&扩展欧几里德
- 扩展欧几里德
- 扩展欧几里德
- 扩展欧几里德
- 扩展欧几里德
- 扩展欧几里德
- 扩展欧几里德
- 扩展欧几里德
- 扩展欧几里德
- 扩展欧几里德
- 扩展欧几里德!!!!!
- 扩展欧几里德
- 扩展欧几里德
- 扩展欧几里德
- MySql中varchar(10)和varchar(100)的区别==>>以及char的利弊
- Android工具类系列-ToastUtil
- 链表后半部分反转(2016亚信实习生笔试题)
- Log4j.properties配置详解
- 批处理脚本编写
- 扩展欧几里德
- 验证当前日期是否为一年中最后一天,否则,默认返回当前日期前一年的最后一天日期
- 百度鹰眼Android SDK 运行配置
- SDO处理(二)
- 解决php网页运行超时问题:Maximum execution time of 30 seconds exceeded
- 转发和重定向的本质区别 简洁明了
- 指针
- JAVA除法保留小数点后两位的方法
- 第13周 阅读程序2-(1)交通工具类