• <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>

            yuanyuelang

            常用鏈接

            統(tǒng)計(jì)

            最新評(píng)論

            數(shù)論(2)-------擴(kuò)展歐幾里得算法

            擴(kuò)展歐幾里得算法------求解線性方程ax+by=c

            1.應(yīng)用:
                  線性方程ax+by=c ,已知a,b,c,求解x,y.

            2.基本思路:
                
                  ax+by=c有解 => c=k*gcd(a,b)=kd(因?yàn)閐=gcd(a,b)=>d|(ax+by))
                  
                  我們先考慮求解  ax+by=d
                    由歐幾里得算法,d=bx'+(a mod b)y'=bx'+(a-[a/b]b)y'=ay'+b(x'-[a/b])y'    
                    則由上述兩式子,我們可以得出 x=y' ,y=x'-[a/b]y'
                    這樣子,在歐幾里得算法添加x,y變量,最后得到解。(可結(jié)合下面代碼源代碼進(jìn)行理解)

                  接下來我們來看看ax'+by'=d和ax+by=c之間的關(guān)系
                  
                    (c/d)ax'+(c/d)by'=(c/d)d 即 可以得到 x=(c/d)x',y=(c/d)y'
                     
                     所以可以得到ax+by的一組解
                     那么ax+by=c所有解的形式是什么呢?
                         a(x+qb)+b(y-qa)=c; q為任意整數(shù)
            (注意,當(dāng)要求y-qa的最小正整數(shù)min時(shí),由y-qa>=0, q取[y/a]最小,min=y-[y/a]y,但是,[y/a]可能為0,如果y是負(fù)數(shù),min此時(shí)也為負(fù)數(shù),不好,此時(shí)令min+=a就可以取得最小正整數(shù)值了([y/a]=0所以|y|<a),這段可以自己找個(gè)例子好好理解下啊)

            3.源代碼模板
             1int Extended_Euclid(int a,int b,int& x,int &y)
             2{
             3    if(b==0){
             4        x=1;
             5        y=0;
             6        return a;
             7    }

             8    int d=Extended_Euclid(b,a%b,x,y);
             9    int temp=x;x=y;y=temp-a/b*y;
            10    return d;
            11}

            12//用擴(kuò)展歐幾里得算法解線性方程ax+by=c;
            13bool linearEquation(int a,int b,int c,int& x,int &y)
            14{
            15    int d=Extended_Euclid(a,b,x,y);
            16    if(c%d) return false;
            17
            18    int k=c/d;
            19    x*=k;y*=k;//求的只是其中一個(gè)解
            20    return true;
            21}

            四.學(xué)習(xí)心得
                擴(kuò)展歐幾里得算法的妙用,值得好好體會(huì),對(duì)整個(gè)過程自己可以找個(gè)例子好好印證一下,肯定印象深刻。去Google一下POJ上的這方面的題目好好實(shí)踐一下吧!!!

            posted on 2009-09-05 16:51 原語餓狼 閱讀(4240) 評(píng)論(4)  編輯 收藏 引用 所屬分類: 數(shù)論

            評(píng)論

            # re: 數(shù)論(2)-------擴(kuò)展歐幾里得算法 2011-02-02 09:34 ding

            頂  回復(fù)  更多評(píng)論   

            # re: 數(shù)論(2)-------擴(kuò)展歐幾里得算法 2011-04-07 23:08 曉楠

            由歐幾里得算法,d=bx'+(a mod b)y'=bx'+(a-[a/b]b)y'=ay'+b(x'-[a/b])y'

            這句話應(yīng)該是
            d=bx'+(a mod b)y'=bx'+(a-[a/b]b)y'=ay'+b(x'-[a/b]y')
              回復(fù)  更多評(píng)論   

            # re: 數(shù)論(2)-------擴(kuò)展歐幾里得算法 2011-10-03 17:34 CrazyForAC

            min=y-[y/a]y錯(cuò)了.  回復(fù)  更多評(píng)論   

            # re: 數(shù)論(2)-------擴(kuò)展歐幾里得算法 2013-12-15 20:43 s

            定理:gcd(a,b) = gcd(b,a mod b) (a>b 且a mod b 不為0)  回復(fù)  更多評(píng)論   

            亚洲香蕉网久久综合影视| 国产精品久久久久国产A级| 91久久精一区二区三区大全| 久久久久高潮综合影院| 日产精品久久久一区二区| 亚洲国产精品成人久久| 久久精品国产福利国产秒| 久久国产一区二区| 麻豆精品久久久久久久99蜜桃| 污污内射久久一区二区欧美日韩| 久久成人国产精品免费软件| 国产精品久久波多野结衣| 综合久久一区二区三区| 国产麻豆精品久久一二三| 三级片免费观看久久| 国产欧美久久久精品| 久久99精品久久久久久久不卡| 精品久久久久久国产牛牛app| 国产精品亚洲综合久久| 久久久免费观成人影院| 亚洲国产精品一区二区久久| 久久人搡人人玩人妻精品首页 | 天堂无码久久综合东京热| 精品久久久噜噜噜久久久| 国产成人综合久久精品红| 久久丫忘忧草产品| 一级女性全黄久久生活片免费 | 久久久久久久97| 久久久精品人妻一区二区三区蜜桃| 久久无码一区二区三区少妇| 久久综合狠狠综合久久97色| 怡红院日本一道日本久久| 日韩人妻无码精品久久久不卡| 亚洲国产精品无码久久SM| 精品国产91久久久久久久a| 国产精品无码久久四虎| 久久人人爽人人爽人人片AV麻豆 | 久久精品国产99久久久| 国产精品一区二区久久国产| 99久久国语露脸精品国产| 精品一二三区久久aaa片|