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

O(1) 的小樂

Job Hunting

公告

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

統(tǒng)計

  • 隨筆 - 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)  編輯 收藏 引用


只有注冊用戶登錄后才能發(fā)表評論。
網(wǎng)站導(dǎo)航: 博客園   IT新聞   BlogJava   博問   Chat2DB   管理


統(tǒng)計系統(tǒng)
青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <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>
            一区二区三区欧美在线| 亚洲午夜日本在线观看| 欧美在线视频观看免费网站| 亚洲日本成人网| 欧美成人r级一区二区三区| 亚洲欧洲一区二区在线观看| 亚洲激情社区| 国产精品久久久久久影视 | 欧美一区二区在线视频| 国产乱人伦精品一区二区 | 欧美在线免费观看| 欧美一级大片在线观看| 影音先锋久久资源网| 欧美激情一区二区三区不卡| 欧美精品国产一区二区| 午夜视频在线观看一区| 久久精品99| 亚洲三级观看| 亚洲综合第一| 亚洲第一黄网| 亚洲视频在线观看三级| 红桃视频欧美| 99日韩精品| 18成人免费观看视频| 亚洲精品中文字幕有码专区| 国产欧美精品一区| 亚洲国产91| 国产亚洲成精品久久| 亚洲狠狠丁香婷婷综合久久久| 国产精品电影观看| 久久综合色播五月| 欧美色一级片| 欧美国产日本韩| 国产视频在线观看一区二区三区| 欧美77777| 国产麻豆日韩| 夜夜嗨av一区二区三区四区 | 国产丝袜一区二区| 亚洲国产精品一区制服丝袜| 国产精品专区一| 亚洲蜜桃精久久久久久久| 国产一区二区三区在线观看精品| 亚洲精品一区在线| 一区二区视频免费在线观看| 一区二区三区久久网| 亚洲国产一区二区三区a毛片| 亚洲尤物在线| 亚洲伊人一本大道中文字幕| 蜜臀99久久精品久久久久久软件| 久久精品夜色噜噜亚洲aⅴ| 欧美日韩在线影院| 亚洲国产精品久久久久| 在线欧美小视频| 欧美一区成人| 久久国产成人| 国产精品系列在线| 亚洲深爱激情| 亚洲视频第一页| 欧美伦理在线观看| 亚洲国产成人精品视频| 亚洲第一天堂av| 久久精品一区中文字幕| 久久久噜噜噜久噜久久| 国产精品影片在线观看| 亚洲一区二区在| 午夜精品剧场| 国产精品有限公司| 午夜精品久久久久久久男人的天堂 | 欧美久久久久免费| 久久精品视频导航| 国产偷国产偷亚洲高清97cao| 亚洲无限av看| 欧美怡红院视频| 国产亚洲精品福利| 欧美综合77777色婷婷| 久久精品电影| 激情亚洲网站| 欧美高清在线| 夜夜狂射影院欧美极品| 亚洲免费视频网站| 国产精品网站在线| 久久超碰97人人做人人爱| 久久一二三四| 日韩视频精品| 国产精品卡一卡二卡三| 午夜在线视频一区二区区别| 久久精视频免费在线久久完整在线看| 国产一区二区三区在线观看视频| 久久精品亚洲精品国产欧美kt∨| 欧美成人嫩草网站| 中文在线不卡视频| 国产精品亚洲综合色区韩国| 久久精品视频免费播放| 亚洲国产成人av| 亚洲一区二区高清| 国内精品久久久久久久影视蜜臀| 免费一级欧美片在线观看| 日韩图片一区| 久久久之久亚州精品露出| 亚洲欧洲精品一区二区| 国产精品扒开腿爽爽爽视频| 久久国产福利国产秒拍| 亚洲黑丝在线| 久久久国产精品一区二区中文| 亚洲国产精品成人va在线观看| 欧美日韩在线综合| 久久久久久电影| 亚洲精品视频啊美女在线直播| 欧美亚洲综合在线| 亚洲国产精品一区二区尤物区| 国产精品久久久久久久久动漫| 久久久久国产免费免费| 亚洲视频你懂的| 亚洲成色www久久网站| 午夜精品理论片| 日韩天天综合| 激情成人av在线| 国产精品日韩专区| 欧美精品三级在线观看| 久久精品国产第一区二区三区最新章节 | 久久久久高清| 亚洲欧美日韩国产综合精品二区| 在线观看欧美日韩| 国产亚洲欧美一级| 国产精品福利在线| 欧美人与性动交cc0o| 久久九九免费| 亚洲影院免费| 一区二区三区四区精品| 亚洲国产成人tv| 欧美国产成人在线| 久久人人九九| 欧美在线一二三四区| 亚洲欧美一区二区激情| 一区二区三区黄色| 99视频国产精品免费观看| 精品二区视频| 精品福利免费观看| 国模私拍一区二区三区| 国产女人精品视频| 国产精品亚洲综合久久| 国产精品网站在线播放| 国产精品久久久久aaaa九色| 欧美日韩不卡合集视频| 欧美精品久久久久久久免费观看 | 亚洲国产综合视频在线观看| 美女主播一区| 欧美成人免费在线| 欧美黄色小视频| 亚洲激情欧美| 亚洲毛片在线观看| 一区二区毛片| 亚洲在线不卡| 久久精品盗摄| 久久亚洲春色中文字幕久久久| 久久一区二区三区四区| 久久综合久久久久88| 欧美大片国产精品| 欧美日韩国产成人在线观看 | 亚洲欧美日韩高清| 欧美一区二区三区男人的天堂| 欧美一区二区精美| 久久五月激情| 欧美激情小视频| 国产精品久久久久久久9999| 国产精品一区二区在线观看| 国产一区二区无遮挡| 亚洲国产一区在线| 亚洲深夜福利在线| 久久精品日韩一区二区三区| 欧美h视频在线| 亚洲精品日日夜夜| 亚洲欧美清纯在线制服| 久久免费一区| 欧美日韩免费观看一区| 国产日本欧美一区二区| 亚洲国产一区二区三区青草影视| 日韩视频免费在线| 性色av一区二区三区在线观看| 久久影视精品| 一本综合久久| 美女视频一区免费观看| 国产精品ⅴa在线观看h| 国内在线观看一区二区三区| 亚洲精品久久久一区二区三区| 午夜精品久久久久影视 | 亚洲欧美日韩久久精品| 久久亚洲综合| 中文欧美日韩| 免费中文字幕日韩欧美| 国产女人aaa级久久久级| 日韩一区二区精品在线观看| 久久国产精品久久久久久| 亚洲高清不卡在线| 欧美在线免费观看亚洲| 欧美日韩一区自拍| 激情成人中文字幕| 亚洲一区激情| 亚洲精品国产精品国自产观看| 久久精品九九|