欧几里得算法与逆元

欧几里得算法与逆元

（一）欧几里得与扩展欧几里得

（I）欧几里得辗转相除求gcd

描述

欧几里得辗转相除求最大公因数的表达式如下

$$gcd(a,b)= \begin{cases} a & b=0\\ gcd(b,a mod b) & b \neq 1 \end{cases}$$

代码

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$\neq$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的一组解$x_0和y_0$，再把$x_0和y_0$同时乘c’就得到了a’x+b’y=c’的一组解，记为xx和yy，即原方程的一组解。方程的所有解可以表示为x=xx+kb’，y=yy-ka’，其中$k\in Z$
模板题 poj1061 青蛙的约会
http://poj.org/problem?id=1061
设青蛙的初始坐标分别为A、B，其他的字母与题中一样，设相遇时间为t
则有

A+