青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品

O(1) 的小樂

Job Hunting

公告

記錄我的生活和工作。。。
<2010年9月>
2930311234
567891011
12131415161718
19202122232425
262728293012
3456789

統計

  • 隨筆 - 182
  • 文章 - 1
  • 評論 - 41
  • 引用 - 0

留言簿(10)

隨筆分類(70)

隨筆檔案(182)

文章檔案(1)

如影隨形

搜索

  •  

最新隨筆

最新評論

閱讀排行榜

評論排行榜

Project Euler 140

Need Patience!!calm down..

We get A(x) = (x + 3x2)/(1 - x - x2) = n is positive integer.
Then:
(n + 3)x2 + (n + 1)x - n = 0
D = 5n2 + 14n + 1 = m2 for some integer m.
Again we get:
5n2 + 14n - (m2 - 1) = 0
D = 5m2 + 44 = t2 for some integer t.
We get more general Pell's equation in the form t2 - 5m2 = 44.
I used http://mathworld.wolfram.com/PellEquation.html and .
First solve "unit" form of this equation r2 - 5s2 = 1. First solution is (9, 4), then use recurence formula r(i + 1) = r(1)r(i) + 5s(1)s(i) and s(i + 1) = r(1)s(i) + s(1)r(i) to get all other solutions.
Then build set of soutions of t2 - 5m2 = 44 using it's set of fundamental solutions:{(t, m)} = {(7, 1), (8, 2), (13, 5), (17, 7), (32, 14), (43, 19)} and using the identity t2 - 5m2 = (t2 - 5m2)(r2 - 5s2) = (tr + 5ms)2 - 5(ts + mr)2 = 44. This gives tr + 5ms for the new value of t.
Finally, if t % 5 == 2 then n = (t - 7)/5 :)
My very simple program is:

Java

[hide code]

public class Problem140 {
    public static void main(String[] args) {
        long g, t, n;
        long[] p = {7, 8, 13, 17, 32, 43};
        long[] q = {1, 2, 5, 7, 14, 19};
        long[] r = {9, 9, 9, 9, 9, 9};
        long[] s = {4, 4, 4, 4, 4, 4};   int count = 5;
        // First 5 values of n
        long sum = 2 + 5 + 21 + 42 + 152;   while (true) {
            for (int k = 0; k < 6; k++) {
                g = r[k];
                r[k] = 9 * r[k] + 20 * s[k];
                s[k] = 9 * s[k] + 4 * g;
                t = p[k] * r[k] + 5 * q[k] * s[k];   if (t % 5 == 2) {
                    n = (t - 7) / 5;
                    sum += n;
                    count++;
                    System.out.println("n(" + count + ")=" + n);
                    if (count == 30) {
                        System.out.println("Sum=" + sum);
                        return;
                    }
                }
            }
        }
    }
}

 

 

Answer:

5673835352990

posted on 2010-09-17 21:15 Sosi 閱讀(225) 評論(0)  編輯 收藏 引用


只有注冊用戶登錄后才能發表評論。
網站導航: 博客園   IT新聞   BlogJava   博問   Chat2DB   管理


統計系統
青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <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>
            亚洲激情第一页| 久久综合久久综合久久综合| 亚洲视频在线观看免费| 亚洲无毛电影| 亚洲欧美日韩另类| 国产乱码精品一区二区三区不卡| 羞羞答答国产精品www一本| 亚洲欧美一区二区三区久久 | 久久久噜久噜久久综合| 99riav国产精品| 国产精品99久久久久久www| 国内精品久久久| 久久福利精品| 欧美成年视频| 久久久亚洲午夜电影| 欧美激情a∨在线视频播放| 久久久久国产一区二区| 亚洲一区二区三区高清 | 欧美精品久久久久久久| 国产精品成人aaaaa网站| 久久综合中文字幕| 欧美自拍丝袜亚洲| 欧美精选午夜久久久乱码6080| 久久久青草婷婷精品综合日韩| 欧美区在线播放| 欧美日韩免费一区| 欧美日产在线观看| 极品尤物久久久av免费看| 国产亚洲成av人片在线观看桃| 久久精品女人| 欧美黄色影院| 美日韩丰满少妇在线观看| 亚洲精品美女在线观看播放| 美女脱光内衣内裤视频久久网站| 久久在线免费视频| 亚洲一区二区三区777| 亚洲欧美日韩天堂一区二区| 亚洲国产成人久久综合一区| 欧美亚洲综合在线| 欧美专区日韩专区| 国产精品嫩草久久久久| 国产日韩精品入口| 国产一区二区久久| 欧美一级视频一区二区| 久久久久一区| 午夜综合激情| 免费成人毛片| 欧美日本韩国一区| 亚洲日本va午夜在线影院| 日韩一本二本av| 亚洲欧美激情诱惑| 欧美激情aaaa| 国产精品国产三级国产普通话三级| 亚洲蜜桃精久久久久久久| 在线亚洲自拍| 欧美在线精品一区| 欧美在线视频免费| 久久久水蜜桃av免费网站| 国产日韩欧美精品一区| 亚洲成人在线视频网站| 99精品福利视频| 亚洲精品四区| 欧美一区二视频| 性欧美1819sex性高清| 麻豆国产va免费精品高清在线| 欧美成人一区二区| 欧美成人在线免费观看| 国产精品久久精品日日| 欧美一区在线看| 欧美成人免费在线观看| 另类av一区二区| 你懂的网址国产 欧美| 欧美午夜视频在线| 亚洲香蕉成视频在线观看| 麻豆国产精品一区二区三区| 久久在精品线影院精品国产| 老牛影视一区二区三区| 国产精品亚洲综合一区在线观看| 欧美在线播放一区二区| 亚洲国产你懂的| 欧美在线你懂的| 亚洲乱码精品一二三四区日韩在线| 久久福利电影| 免费看的黄色欧美网站| 国产一区二区三区在线观看免费视频| 国产精品久久久免费| 久久精品国产精品亚洲综合| 亚洲天堂网站在线观看视频| 夜夜爽夜夜爽精品视频| 美女日韩欧美| 欧美亚州一区二区三区| 久久午夜电影网| 亚洲午夜在线观看视频在线| 一区视频在线播放| 久久精品女人| 免费日韩av电影| 亚洲免费在线观看视频| 最新69国产成人精品视频免费| 国产精品乱子久久久久| 亚洲视频一二| 久久全国免费视频| 欧美精品日韩| 日韩亚洲欧美一区| 久久精品卡一| 亚洲国产婷婷| 在线播放一区| 亚洲无限av看| 国产精品视频免费| 亚洲精品日产精品乱码不卡| 欧美激情视频在线免费观看 欧美视频免费一| 亚洲一区观看| 亚洲娇小video精品| 亚洲第一在线综合在线| 国产精品视频你懂的| 欧美亚洲免费高清在线观看| 欧美1区免费| 国产真实乱子伦精品视频| 久久免费国产精品1| 国产精品精品视频| 久久精品av麻豆的观看方式| 亚洲精品视频一区| 亚洲精品午夜| 一本色道久久综合亚洲精品婷婷 | 久久精品女人天堂| 欧美日韩精品免费看| 亚洲社区在线观看| 亚洲一二三区在线观看| 一区二区三区毛片| 亚洲一区二区三区色| 一个色综合av| 亚洲欧美日韩一区二区三区在线| 国产视频久久久久| 性色av一区二区三区在线观看| 黄色精品网站| 久久国产精品72免费观看| 亚洲国产乱码最新视频| 亚洲精品永久免费精品| 99精品视频一区| 亚洲综合色视频| 欧美中文字幕视频| 久久九九热re6这里有精品| 久久久久久亚洲精品中文字幕| 久久精品一区二区三区四区 | 亚洲自拍偷拍福利| 一区二区三区久久精品| 亚洲图色在线| 欧美一区二区三区婷婷月色 | 国产亚洲美州欧州综合国| 亚洲欧美日韩一区| 亚洲免费观看在线视频| 欧美日韩一区二区三区| 久热精品在线视频| 亚洲在线一区| 久久精品亚洲乱码伦伦中文| 免费在线成人av| 亚洲人www| 精品成人一区二区| 免费人成精品欧美精品| 亚洲精品孕妇| 欧美一区2区视频在线观看 | 欧美在线观看你懂的| 久热这里只精品99re8久| 亚洲激情不卡| 亚洲香蕉在线观看| 欧美xart系列在线观看| 欧美精品www| 亚洲经典三级| 国产精品xxxxx| 亚洲国产精品毛片| 在线综合亚洲欧美在线视频| 久久综合久久综合久久综合| 国内久久视频| 欧美一级专区| 欧美在线观看一二区| 亚洲日本中文字幕区| 欧美一区二区三区视频在线| 欧美刺激性大交免费视频| 国产精品一国产精品k频道56| 一区二区三区四区国产| 国产嫩草一区二区三区在线观看| 91久久亚洲| 欧美在线首页| 国产精品欧美在线| 欧美福利视频在线| 另类综合日韩欧美亚洲| 日韩一区二区精品葵司在线| 另类尿喷潮videofree| 又紧又大又爽精品一区二区| 欧美日韩一区在线观看| 亚洲美女av网站| 乱中年女人伦av一区二区| 欧美高清视频一区| 亚洲欧美一区二区三区久久| 亚洲国产欧美日韩另类综合| 在线观看国产欧美| 欧美三区视频| 欧美一区二区播放| 亚洲黄色免费网站| 狠狠久久综合婷婷不卡| 香港成人在线视频|