欧几里得算法与逆元
(一)欧几里得与扩展欧几里得
(I)欧几里得辗转相除求gcd
描述
欧几里得辗转相除求最大公因数的表达式如下
gcd(a,b)={
agcd(b,amodb)b=0b≠1
代码
int gcd(int a,int b)
{
if(b==0)return a;
else return gcd(b, a%b);
}
(II)扩展欧几里得
描述
定理:对于任意整数a,b,都存在一组整数x、y使得ax+by=gcd(a,b)成立
于是设 bx’+(a%b)y’=gcd(b,a%b)
则ax+by=bx’+(a%b)y’
可以推导出x=y’,y=x’-[a/b]*y’
当b=0,a ≠ 0时,gcd(a,b)=a,则x=1,y=0,
由此便可以向上推出ax+by=gcd(a,b)中x,y的值
代码
#define ll long long
ll exgcd(ll a, ll b, ll &x, ll &y)
{
if(!b){
x=1;y=0;return a;}
ll t, d;
d=exgcd(b,a%b,x,y);
t=x;x=y;y=t-a/b*y;
return d;
}
求解方程ax+by=c
定理:当且仅当当gcd(a,b)|c时,方程有整数解,设d=gcd(a,b),令a’=a/d,b’=b/dc’=c/d,方程两边同时除以d就得到a’x+b’y=c’,我们可以先用扩展欧几里得求出a’x+b’y=1的一组解 x0和y0 ,再把 x0和y0 同时乘c’就得到了a’x+b’y=c’的一组解,记为xx和yy,即原方程的一组解。方程的所有解可以表示为x=xx+kb’,y=yy-ka’,其中 k∈Z
模板题 poj1061 青蛙的约会
http://poj.org/problem?id=1061
设青蛙的初始坐标分别为A、B,其他的字母与题中一样,设相遇时间为t
则有
A+