• <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>

            The Fourth Dimension Space

            枯葉北風寒,忽然年以殘,念往昔,語默心酸。二十光陰無一物,韶光賤,寐難安; 不畏形影單,道途阻且慢,哪曲折,如渡飛湍。斬浪劈波酬壯志,同把酒,共言歡! -如夢令

            歐幾里德算法及其擴展算法的迭代版本(轉)

            整除和最大公約數。

            如果m能整除n,記m|n

            對于不全為0的兩個整數ab,能同時整除他們倆的最大整數為它們的最大公約數,記為gcd(a,b)。如果gcd(a,b)=1,則ab互質。

            gcd(a,b)可用歐幾里德算法:

            a = q1 * b + r1

            b = q2 * r1 + r2

            r1 = q3 * r2 + r3

            r2 = q4 * r3 + r4

            ……

            rn-3 = qn-1 * rn-2 + rn-1

            rn-2 = qn * rn-1 + rn => rn就是gcd

            rn-1 = qn+1 * rn + 0

            歐幾里德算法非常好寫。以pascal語言為例:

            function gcd(a,b:longint):longint;

            begin

            if b=0 then gcd:=a

            else gcd:=gcd(b,a mod b);

            end;

            線性方程ax+by=c的解法。

            對于abax+byab的線性組合。這里xy可以是任何整數。

            g=gcd(a,b),則g整除所有的ax+by。特別地,所有ax+by的取值中,最小的正值為g

            因此,方程ax+by=c,僅當gcd(a,b)|c時有解。因此只考慮方程ax+by=gcd(a,b)的解法。

            回憶上一章中講的歐幾里德算法,在第i步,我們有

            ri-2 = qi * ri-1 + ri 。若ri+1=0,ri就是gcd

            假設ri可表示為ab的線性組合(xi,yi),即ri=axi+byi,則有(xi,yi)=(xi-2,yi-2)-qi* (xi-1,yi-1)

            初始時, 由于a=q1*b+r1,即r1=a-b*q1,由序列{Rn}的通項公式Rn=R(n-2)-R(n-1)*q1,我們可以把a看做r(-1),

            b看做,r0,那么對應于表達式ri=a*xi+b*yi,r(-1)=a=a*1+b*0=(1,0),r0=b=a*0+b*1=(0,1);

            即r-1=a=(1,0)r0=b=(0,1)。這樣便可一步步推出rn=(xn,yn),即為原方程的一組解。

            因此有如下定理:

            ab為非零整數,g=gcd(a,b),則方程

            ax+by=g 總可以通過歐幾里德算法找出一組可行整解,記為(x1,y1)。它的通解可以表示為:

            (x,y)=(x1+k*b/g , y1-k*a/g) ,其中k為任意整數。

            競賽中常用的歐幾里德擴展算法采用的是倒推的方法,相比而言,本章提供的方法需要的變量稍多,常數稍大(只大了一點點)。但由于易于用迭代實現,而且實際測試中兩種算法時間差別十分微小,因此本算法仍有應用價值。

            posted on 2009-07-10 17:59 abilitytao 閱讀(591) 評論(0)  編輯 收藏 引用

            久久www免费人成看国产片| 久久久九九有精品国产| 国产精品99精品久久免费| 人妻无码αv中文字幕久久琪琪布 人妻无码久久一区二区三区免费 人妻无码中文久久久久专区 | 免费观看久久精彩视频| 久久国产精品成人免费| 久久精品国产精品亚洲| 五月丁香综合激情六月久久| 国内精品久久久久久中文字幕| 一97日本道伊人久久综合影院| 一本色道久久HEZYO无码| 97精品久久天干天天天按摩| 国产精品日韩欧美久久综合| 精品久久久久久久久午夜福利| 香港aa三级久久三级老师2021国产三级精品三级在 | 国产91色综合久久免费| 国产精品一区二区久久精品涩爱| 久久久精品午夜免费不卡| 狠狠色噜噜色狠狠狠综合久久| 久久www免费人成精品香蕉| 69久久精品无码一区二区| 日韩人妻无码精品久久久不卡| 日产久久强奸免费的看| 94久久国产乱子伦精品免费| 久久久九九有精品国产| 国产精品一区二区久久精品| 天天躁日日躁狠狠久久| 久久亚洲AV成人无码国产| 国产A级毛片久久久精品毛片| 亚洲午夜无码久久久久小说| 久久亚洲国产成人精品无码区| 久久线看观看精品香蕉国产| 成人资源影音先锋久久资源网| 亚洲精品国产字幕久久不卡| 五月丁香综合激情六月久久| 久久综合精品国产二区无码| 日产精品久久久一区二区| AV无码久久久久不卡蜜桃| 久久久精品国产免大香伊| 久久精品麻豆日日躁夜夜躁| 久久亚洲精品人成综合网|