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

            ArcTan

            dfs
            隨筆 - 16, 文章 - 117, 評論 - 6, 引用 - 0
            數據加載中……

            poj 1061(線性同余)--青蛙的約會

            http://poj.org/problem?id=1061

            剛剛學了點數論,也好幾天沒有寫代碼了,生疏了不少,總是有錯誤,編譯個程序怎么都那么糾結啊,看來水平實在是差啊!!
            線性同余,方程還是很好就構造出來了,extended_euclid,中國剩余定理,可是后面卻不指導怎么求最小的整數了:
            看了看大牛的題解,原來這樣啊,自己數論還得多學學才行,多想想啊:

            分析:設青蛙跳了k次,那么就有(x+mk)-(y+nk)=p*L.

            即x-y+(m-n)k=p*L,即(m-n)*k≡(y-x) (mod L).

            這個線性同余方程有解當且僅當gcd(m-n,L)|(y-x).

            令a=m-n,b=L,c=y-x.用擴展歐幾里得解方程ax+by=c.

            可以求出原方程的一個解.如何求最小正整數解呢?

            假設我們已經得到一個x0,令d=gcd(m-n,L),

            那么所有解可以表示為x=x0+k*L/d.

            設L'=L/d.

            Xmin=(x0 mod L'+L') mod L'.


            WA兩次,0Ms,囧,還有一次編譯錯誤!!!

            #include<stdio.h>
            #include
            <string.h>
            #include
            <math.h>
            long long  c,d;
            long long gcd_ext(long long a,long long b)
            {
                
            long long gcd,t;
                
            if (!b)
                {
                    c
            =1;d=0;
                    
            return a;
                }
                gcd
            =gcd_ext(b,a%b);
                t
            =c;c=d;d=t-a/b*d;
                
            return gcd;
            }
            int main()
            {
                
            long long x,y,m,n,L,a,b,gcd;
                
            while (scanf("%I64d%I64d",&x,&y)==2)
                {
                    scanf(
            "%I64d%I64d%I64d",&m,&n,&L);
                    a
            =m>n?m-n:n-m;
                    b
            =m>n?y-x:x-y;
                    gcd
            =gcd_ext(a,L);
                    L
            =L/gcd;
                    
            if (b%gcd==0)
                        printf(
            "%I64d\n",((c*b/gcd)%L+L)%L);
                    
            else
                        printf(
            "Impossible\n");
                }
                
            return 0;
            }


            總結:代碼,還是得天天寫,三日不練手生。自己多想想,多思考思考才能提升能力哈。

                     不要總是去看別人的題解,要有自己的思路哈。

                     數論,還得繼續看,繼續學。要吃透才行。


            posted on 2012-04-12 00:18 wangs 閱讀(582) 評論(0)  編輯 收藏 引用 所屬分類: ACM-模擬

            国产精品99久久久久久宅男| 久久精品国产精品亜洲毛片 | 久久人人妻人人爽人人爽| 久久久SS麻豆欧美国产日韩| 久久久无码精品亚洲日韩按摩 | 久久丫精品国产亚洲av不卡| 伊人久久大香线焦综合四虎| 一极黄色视频久久网站| 精品久久久久久无码专区不卡| 久久久青草青青国产亚洲免观| 嫩草伊人久久精品少妇AV| 久久精品国产亚洲精品| 久久夜色精品国产噜噜噜亚洲AV| 久久精品中文字幕有码| 国内精品久久久久久99| 久久精品国产亚洲AV蜜臀色欲| 久久这里只精品国产99热| 久久精品国产亚洲精品2020 | 人妻精品久久久久中文字幕69 | 亚洲人成电影网站久久| 99久久精品国产一区二区| 国内精品久久久久伊人av| 2021最新久久久视精品爱 | 欧美精品国产综合久久| 久久久久99精品成人片三人毛片 | 国产亚洲精品久久久久秋霞| 久久精品国产黑森林| 亚洲国产成人久久精品影视| 国产精品一区二区久久精品| 久久久久99精品成人片欧美| 日本久久久久亚洲中字幕| 亚洲香蕉网久久综合影视| 99久久国产综合精品女同图片| 久久久久久av无码免费看大片| 国产高潮国产高潮久久久91| 国内精品久久久久久久影视麻豆 | 久久精品国产亚洲AV高清热| 韩国免费A级毛片久久| 97久久国产亚洲精品超碰热| 国产精品视频久久| 久久国产精品一区|