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

            poj 1006 Biorhythms 中國(guó)剩余定理

               此題本來(lái)模擬即可,但是注意有容易出錯(cuò)的地方。
               這題主要是可以用中國(guó)剩余定理來(lái)做。
               根據(jù)題意可以抽象出這樣的模型。給出三個(gè)數(shù)A,B,C分別是模23,28,33后的余數(shù),求最小的數(shù)字
            使得其模23,28,33分別為A,B,C,并且要大于給定的數(shù)字D。
               中國(guó)剩余定理很好的解決了這種余數(shù)問(wèn)題。令模數(shù)為Ni,余數(shù)為Ai,設(shè)Mi = N1*N2*...*Ni-1*Ni+1*...*Nn,
            那么答案一定滿足形式ans = ΣAi*Mi*(Mi對(duì)Ni的乘法逆) % N。(N為所有Ni的乘積)。
               很明顯,由于ans的第i項(xiàng)有Mi因子,所以模N1-Ni-1和Ni+1-Nn肯定是0,而Ai*Mi*(Mi對(duì)Ni的乘法逆) %Ni
            就是Ai。這樣就滿足了要求。
               代碼如下:
            #include <stdio.h>
            #include <algorithm>
            #include <string.h>
            #include <vector>
            using namespace std;

            int Egcd(int nN, int nM, int& nX, int& nY)
            {
                if (nM == 0)
                {
                    nX = 1, nY = 0;
                    return nN;
                }
                int nRet = Egcd(nM, nN % nM, nX, nY);
                int nT = nX;
                nX = nY;
                nY = nT - (nN / nM) * nY;
                return nRet;
            }

            int main()
            {
                int nA, nB, nC, nD;
                int nDays = 21252;
                int nCase = 1;
                
                while (scanf("%d%d%d%d", &nA, &nB, &nC, &nD),
                       nA != -1 || nB != -1 || nC != -1 || nD != -1)
                {
                    int nFirst = 0;
                    nA %= 23;
                    nB %= 28;
                    nC %= 33;
                    int nM1= 28 * 33, nM2 = 23 * 33, nM3 = 23 * 28;
                    int nN1, nN2, nN3, nTemp;
                    Egcd(23, nM1, nTemp, nN1);
                    Egcd(28, nM2, nTemp, nN2);
                    Egcd(33, nM3, nTemp, nN3);
                    nFirst = (nA * nM1 * nN1 + nB * nM2 * nN2 + nC * nM3 * nN3) % nDays;
                    while (nFirst <= nD)nFirst += nDays;
                    printf("Case %d: the next triple peak occurs in %d days.\n",
                           nCase++, nFirst - nD);
                }
                
                return 0;
            }

            posted on 2012-10-03 23:11 yx 閱讀(883) 評(píng)論(0)  編輯 收藏 引用 所屬分類(lèi): 數(shù)論

            <2012年9月>
            2627282930311
            2345678
            9101112131415
            16171819202122
            23242526272829
            30123456

            導(dǎo)航

            統(tǒng)計(jì)

            公告

            常用鏈接

            留言簿(3)

            隨筆分類(lèi)

            隨筆檔案

            me

            好友

            同學(xué)

            網(wǎng)友

            搜索

            最新評(píng)論

            閱讀排行榜

            評(píng)論排行榜

            久久亚洲国产精品五月天婷| 精品综合久久久久久98| 亚洲国产精品久久久久婷婷软件 | 亚洲欧美一级久久精品| 久久婷婷人人澡人人爽人人爱| 色诱久久久久综合网ywww| 93精91精品国产综合久久香蕉 | 香蕉久久一区二区不卡无毒影院| 国产综合精品久久亚洲| 人妻少妇久久中文字幕一区二区| 久久国产精品免费| 69SEX久久精品国产麻豆| 亚洲午夜无码AV毛片久久| 999久久久国产精品| 精品乱码久久久久久久| 精品久久久一二三区| 国产午夜精品久久久久九九| 久久99国产精品久久99果冻传媒 | 久久中文骚妇内射| 狠狠色丁香婷婷久久综合| 久久久久亚洲?V成人无码| 一本伊大人香蕉久久网手机| 97r久久精品国产99国产精| 伊人久久精品无码av一区| 一本一本久久a久久精品综合麻豆| 亚洲国产精久久久久久久| 色综合久久最新中文字幕| 久久久一本精品99久久精品66| 久久国语露脸国产精品电影| 思思久久好好热精品国产| 久久精品国产欧美日韩99热| 久久久久久久久久免免费精品| 久久精品亚洲乱码伦伦中文| 久久国产精品免费| 深夜久久AAAAA级毛片免费看| 三级片免费观看久久| 久久人妻少妇嫩草AV蜜桃| 色青青草原桃花久久综合| 久久无码AV一区二区三区| 97久久国产综合精品女不卡| 久久人妻少妇嫩草AV无码专区 |