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

O(1) 的小樂

Job Hunting

公告

記錄我的生活和工作。。。
<2025年9月>
31123456
78910111213
14151617181920
21222324252627
2829301234
567891011

統計

  • 隨筆 - 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 閱讀(221) 評論(0)  編輯 收藏 引用

統計系統
青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <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>
            能在线观看的日韩av| 美女黄网久久| 欧美成人日本| 免费视频一区| 亚洲国产高潮在线观看| 久久精品日韩| 欧美成人精品一区二区三区| 免费观看一区| 亚洲国产cao| 亚洲人成人99网站| 99国产精品久久久久久久成人热| 亚洲美女视频网| 中国成人黄色视屏| 欧美在线免费视频| 看片网站欧美日韩| 欧美三区在线视频| 国产永久精品大片wwwapp| 亚洲福利电影| 中文亚洲免费| 久久免费视频网| 欧美国产三级| 亚洲午夜精品一区二区三区他趣| 老司机免费视频久久| 亚洲乱码精品一二三四区日韩在线 | 国产一区二区在线观看免费播放| 国产日韩欧美在线观看| 亚洲激情视频在线播放| 欧美一区成人| 亚洲日本aⅴ片在线观看香蕉| 亚洲一区二区三区国产| 久久综合狠狠综合久久综合88 | 狠狠干综合网| 国产精品激情av在线播放| 国产日韩欧美综合| 亚洲免费精品| 老司机午夜精品视频| 妖精视频成人观看www| 久久久www成人免费无遮挡大片| 欧美人成免费网站| 91久久精品国产91久久性色tv| 性欧美8khd高清极品| 亚洲精品欧美专区| 欧美高清视频一区| 亚洲国产精品一区二区第一页| 亚洲欧美国产高清va在线播| 最新成人av网站| 久久久五月天| 国模 一区 二区 三区| 午夜精品福利一区二区蜜股av| 亚洲精品国产精品乱码不99 | 午夜欧美大尺度福利影院在线看| 欧美成人亚洲| 亚洲高清在线视频| 久久青草久久| 久久久福利视频| 国产一区二区三区四区在线观看 | 欧美国产精品| 亚洲国产欧美日韩精品| 美女日韩欧美| 久久综合色播五月| 亚洲国产精品国自产拍av秋霞 | 久久九九99视频| 性做久久久久久免费观看欧美| 国产精品美女久久福利网站| 一区二区三区.www| 日韩网站免费观看| 国产精品成人国产乱一区| 亚洲女同精品视频| 亚洲天堂av图片| 欧美日韩视频在线观看一区二区三区| 亚洲免费高清| 99在线精品视频在线观看| 国产精品久久久久久久久久ktv| 亚洲中无吗在线| 亚洲主播在线| 狠狠色狠狠色综合人人| 欧美国产精品va在线观看| 欧美国产综合视频| 中文亚洲免费| 久久精品国产91精品亚洲| 亚洲激情不卡| 国产精品99久久久久久久vr | 国产精品99久久久久久www| 欧美一区成人| 亚洲一二三区视频在线观看| 欧美激情一区二区| 欧美激情自拍| 亚洲欧美日韩国产一区二区三区 | 国产精品日韩精品| 久久久综合香蕉尹人综合网| 久久夜色精品国产欧美乱| 一区二区欧美在线| 午夜欧美精品| 91久久久一线二线三线品牌| 一区二区三区视频在线| 一色屋精品视频在线看| 一区二区三区产品免费精品久久75| 国产一区在线视频| 一区二区三区av| 亚洲欧洲一区二区三区久久| 亚洲深夜福利视频| 亚洲激情影视| 欧美一区在线直播| 亚洲视频自拍偷拍| 久久精品毛片| 亚洲一区中文| 欧美精品电影在线| 久久久久久亚洲综合影院红桃| 欧美人交a欧美精品| 久久亚洲私人国产精品va| 欧美啪啪成人vr| 麻豆视频一区二区| 国产深夜精品| 亚洲一区成人| 日韩视频永久免费| 久久久精品午夜少妇| 久久国产精品久久久久久电车| 欧美激情综合色| 老司机午夜精品视频在线观看| 国产精品你懂的在线欣赏| 亚洲国产成人午夜在线一区| 很黄很黄激情成人| 国产精品99久久久久久久vr| 亚洲第一久久影院| 久久精品国产一区二区电影| 香蕉久久夜色精品国产使用方法| 欧美日韩大陆在线| 亚洲精品国产视频| 日韩亚洲在线观看| 女人香蕉久久**毛片精品| 免播放器亚洲| 亚洲第一网站| 久久综合九色欧美综合狠狠| 久久精品亚洲一区二区三区浴池| 国产精品久久久久久久午夜| 一区二区日本视频| 亚洲一品av免费观看| 欧美三区美女| 日韩一级片网址| 亚洲午夜av| 国产精品观看| 亚洲一区二区三区久久| 亚洲综合社区| 国产免费一区二区三区香蕉精| 中国日韩欧美久久久久久久久| 亚洲视频免费在线| 9久re热视频在线精品| 夜夜狂射影院欧美极品| 欧美日韩在线高清| 一区二区三区 在线观看视频| 亚洲女性裸体视频| 国产亚洲a∨片在线观看| 久久av红桃一区二区小说| 乱人伦精品视频在线观看| 亚洲国产精品视频一区| 欧美理论电影在线播放| 亚洲视频福利| 欧美一区二区三区男人的天堂| 国产日韩av高清| 久久精品色图| 亚洲日韩欧美视频| 亚洲综合色网站| 好吊色欧美一区二区三区四区| 免费视频一区| 99在线观看免费视频精品观看| 先锋影音网一区二区| 尤物精品国产第一福利三区 | 亚洲最新中文字幕| 久久久久免费| 中文av一区二区| 黄色av日韩| 欧美日韩在线观看一区二区| 欧美一区=区| 亚洲精品一区久久久久久| 久久精品国产成人| 一区二区三区你懂的| 狠色狠色综合久久| 国产精品成人国产乱一区| 久久久久青草大香线综合精品| 9色精品在线| 久久蜜臀精品av| 亚洲午夜一级| 亚洲国产成人久久综合一区| 国产精品成人v| 欧美电影在线播放| 久久精品国产清高在天天线| 一本久久a久久免费精品不卡| 久久亚洲风情| 午夜视频在线观看一区二区| 亚洲精品免费在线播放| 国产日韩亚洲欧美精品| 欧美日韩国产精品专区| 久久爱www.| 亚洲综合国产精品| 亚洲乱码视频| 亚洲精品视频免费观看| 亚洲第一免费播放区| 久久婷婷一区| 久久av一区二区三区| 亚洲午夜黄色|