暑假集训之gcd&exgcd简介

gcd算法就是平时所说的欧几里得算法，主要是来算最大公约数的，其原理是前面所提到的辗转相除法，这里就不再赘述了。除此之外，求最小公倍数的问题也是借助这个算法实现，具体的实现方法是：

a，b最小公倍数=（a*b）/gcd(a,b);通常写为a/gcd(a,b)*b(防止运算中间溢出)；

gcd算法的模板

1、递归算法

int gcd(int x,int y)
{
    return y?gcd(y,x%y):x;
}

2，非递归算法

在求解多个数的最小公倍数的问题是，通常采用两辆比较的方式

模板

int work(a[],x)
{
    int ans=a[0];
    for(int i=1;i<x;i++)
    {
        ans=gcd(a[i],ans);
    }
    return ans;
}

求多个数的最大公约数同理！

--------------------------扩展欧几里得算法--------------------------------------

扩展欧几里得求的是一个方程的解
方程如下：
ax + by = gcd(a,b)(类似这样的方程)
现在我们假设 a > b,接下来就要有两种情况
1) 那么当b == 0的时候gcd(a,b)=a;
那么我们现在要解的这个方程就是
ax=a
==> x = 1, y = 0;

2)当b!=0的时候，我们可以设：
ax1+by1=gcd(a,b)
bx2+(a%b)y2=gcd(b,a%b)
又因为：gcd(a,b)=gcd(b,a%b)
由所以我们可以得到：==>ax1+by1=bx2+(a%b)y2
那么现在a%b又可以写成a-floor(a/b)*b ---[这里面floor()是向下取整的意思]
==>ax1+by1=bx2+[a-floor(a/b)*b]*y2

==>x1*a+y1*b=y2*a+[x2-floor(a/b)*y2]*b
我们现在让系数相等也就是a和b当作未知数
==> x1 = y2
==> y1 = x2-floor(a/b)*y2
那么我们就得到了扩展欧几里得的递归算法
方程a*x+b*y=c
有解的情况是:c mod gcd(a,b)==0
模板：

void Ex_gcd(int a, int b, int &x, int &y)
{
    if(b == 0)//递归出口
    {
        x = 1;
        y = 0;
        return;
    }
    int x1, y1;
    Ex_gcd(b, a%b, x1, y1);
    x = y1;
    y = x1-(a/b)*y1;
}（模板这东西每个人都对问题有自己的见解，所以模板也不尽相同，网上的gcd模板也是多种多样，可以随便找一个顺手的，当然最好事根据自己的理解编一个自己的模板）