• <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>
            隨筆 - 68  文章 - 57  trackbacks - 0
            <2010年3月>
            28123456
            78910111213
            14151617181920
            21222324252627
            28293031123
            45678910

            常用鏈接

            留言簿(8)

            隨筆分類(74)

            隨筆檔案(68)

            搜索

            •  

            最新評論

            閱讀排行榜

            評論排行榜

              是個經典問題了,求滿足a ^ 2 + b ^ 2 = c ^ 2且a、b、c兩兩互質的三元組個數,其中a、b、c都小于等于n,同時還要求出n以內不屬于三元組(不必互質)的數的個數,其中n <= 1000000。
              數據范圍很大,單純的n ^ 2的枚舉肯定是不行的。我發現互質的三元組不多,想找一個辦法求出互質的,然后篩出其他的三元組(乘以倍數),但是仍然想不出怎樣很快找到互質的三元組。到網上查了一下,找到了一個很好的網站(http://xserve.math.nctu.edu.tw/people/cpai/carnival/fraction/05.htm)講解這個問題。
              首先只需找到互質的a和b就可以了,其余的可以通過乘以一定的倍數得到。首先a和b必須不能都是奇數,否則的話,不妨設a = 2p + 1, b = 2q + 1, c = 2r。帶入a ^ 2 + b ^ 2 = c ^ 2有:2 * (r ^ 2 - p ^ 2 - q ^ 2 - p - q) = 1,出現了矛盾。這樣必然a和b一奇一偶,設a是偶數。恒等變換一下:(這里是x ^ 2 + y ^ 2 = z ^ 2,x是偶數)


              令u = (z - y) / 2,v = (z + y) / 2。因為y、z互質,那么u、v互質。如果不互質,有v = ku,變換后得到:z / y = (k + 1) / (k - 1)。因為z和y互質,這樣這個等式的解只能是z = k + 1且y = k - 1。這樣的話u = 1。可以理解為這也是互質的特殊情況(gcd = 1)。這樣(x / 2) ^ 2 = u * v,u和v沒有公共的質因數,因此必然u和v都是完全平方數。接下來的事情就比較簡單了,令u = m ^ 2,v = n ^ 2,然后利用m和n表示x、y、z有:y = n ^ 2 - m ^ 2,x = 2 * m * n,z = n ^ 2 + m ^ 2。這樣題目中數據范圍1000000,而枚舉m和n的范圍最大1000,這個復雜度就可以接受了。枚舉的時候不用擔心出現重復的互質的x與y,因為這里x一定是偶數,如果x和y已經互質,那么y一定是奇數,解是唯一的。
              很囧的是這個巧妙的方法居然公元前250年就有大牛想到了(丟番圖),實在orz!

             

            posted on 2009-06-09 09:09 sdfond 閱讀(333) 評論(0)  編輯 收藏 引用 所屬分類: Algorithm - Number Theory
            亚洲AV无码久久寂寞少妇| 久久精品国产亚洲Aⅴ蜜臀色欲 | 国产综合久久久久| 国产成人久久精品激情| 99久久精品免费看国产| 久久久久精品国产亚洲AV无码| 韩国免费A级毛片久久| 久久精品国产WWW456C0M| 亚洲精品tv久久久久久久久| 久久精品成人影院| 91精品国产高清久久久久久io| 久久久久亚洲精品男人的天堂| 国产产无码乱码精品久久鸭| 亚洲欧美日韩中文久久| 久久久久久久亚洲精品| 久久亚洲国产精品一区二区| 久久精品国产亚洲av麻豆蜜芽| 久久精品无码一区二区三区日韩| 国产亚洲欧美精品久久久| 国产精品久久久久久久久久影院| 久久综合九色综合97_久久久| 中文字幕乱码人妻无码久久| 免费一级做a爰片久久毛片潮| 久久综合九色综合久99 | 成人国内精品久久久久影院VR| 久久天天躁狠狠躁夜夜不卡 | 99re久久精品国产首页2020| 日韩欧美亚洲综合久久 | 久久精品人人做人人爽电影蜜月| 亚州日韩精品专区久久久| 精品久久久久久久中文字幕| 91精品免费久久久久久久久| 精品久久久久久亚洲| 97精品久久天干天天天按摩 | 亚洲精品无码久久毛片| 久久国产香蕉一区精品| 久久综合久久鬼色| 久久人人爽人人爽人人片av麻烦| 中文精品久久久久人妻| 久久久久亚洲AV无码专区首JN| 久久人人爽人人爽人人片av麻烦|