码迷,mamicode.com
首页 > 编程语言 > 详细

扩展欧几里得算法+推论

时间:2019-01-15 10:41:35      阅读:142      评论:0      收藏:0      [点我收藏+]

标签:代码   time   算法   相减   等于   取值   long   上下   return   

什么是扩展欧几里得?

扩展欧几里得算法是建立在欧几里得算法(gcd)之上。

首先,我们知道有\(a*x+b*y=gcd(a,b)\)

我们怎么求这个\(x,y\)呢?

这时候我们就得使用exgcd算法,我们来推导一下吧!

\(a*x+b*y=gcd(a,b)\)
\(a*x+b*y=gcd(b,a\% b)\)
\(a*x+b*y=b*x'+(a- \left \lfloor \frac{a}{b} \right \rfloor*b)*y'\)
\(a*x+b*y= a*y'+b*(x'-\left \lfloor \frac{a}{b} \right \rfloor*y')\)

因此,根据系数对应,我们得到了

\(x=y'\),\(y=x'-\left \lfloor \frac{a}{b} \right \rfloor*y'\)

那这个式子我们显然可以在递归里面顺带算出来嘛。

再想一想,我们gcd的递归出口为\(b=0\),即\(a*x=gcd(a,b)\),所以说我们的\(x,y\)的递归出口也为\(x=1,y=0\)

代码大概长这样qwq:

long long exgcd(long long A,long long B,long long &x,long long &y)
{
    if(B==0) 
    {
        x=1,y=0;
        return A;
    }
    long long t=exgcd(B,A%B,x,y),tx=x;
    x=y;
    y=tx-(A/B)*y;
    return t;
}

酱紫,我们就求出了\(x,y\)的一组解\(x_0,y_0\)


进一步推导

如果我们要求x的最小正整数解,那不免要求x的通项公式。

首先我们的推导建立在已经求出了一组\((x_0,y_0)\)使得\(a\times x+b\times y=gcd(a,b)\)

我们要求的\(x\)的通项公式是建立在\(x_0\)之上的,我们假设\(x=x_0+p\)\(y=y_0-q\)

现在问题就变为了如何求这个\(p\)

原式变为:
\[a(x_0+p)+b(y_0-q)=gcd(a,b)\]

展开得:
\[a*x_0+a*p+b*y_0-q*b=gcd(a,b)\]

与原式\(a*x_0+b*y_0=gcd(a,b)\)相减得
\[ap=bq\]

\[p=\frac{b*q}{a}\]

我们设\(d=gcd(a,b)\),有\(a=a'*d\),\(b=b'*d\)

将上一个式子上下同时除以\(d\)
\[p=\frac{b'*q}{a'}\]

因为\(p\)为正整数,因此我们的\(q\)至少等于\(a'\)才能使得\(p\)的取值最小
\[p=b*\frac{a'}{a}=b*\frac{1}{d}=b*\frac{1}{gcd(a,b)}\]

解毕,我们得到了\(x\)的通项公式\(x=x_0+k*\frac{b}{gcd(a,b)}\)

.

完结撒花(*′???`)?

扩展欧几里得算法+推论

标签:代码   time   算法   相减   等于   取值   long   上下   return   

原文地址:https://www.cnblogs.com/GoldenPotato/p/10269979.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!