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

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>
            久久久久成人精品| 欧美日韩妖精视频| 一区在线观看| 欧美aa国产视频| 欧美肥婆在线| 国产精品久久久久一区| 伊伊综合在线| 欧美黄色大片网站| 欧美日韩国产专区| 久久不射电影网| 久久久噜噜噜| 一区二区免费在线视频| 亚洲自拍电影| 伊人天天综合| 夜夜嗨av一区二区三区网站四季av | 国产日本欧美一区二区三区在线 | 欧美精品手机在线| 亚洲一二三级电影| 欧美在线视频免费播放| 欧美日韩中国免费专区在线看| 国产麻豆一精品一av一免费| 久久这里只有精品视频首页| 欧美日韩大片| 久久久久久久一区| 欧美日韩国产不卡| 久久久亚洲国产美女国产盗摄| 亚洲视频999| 在线免费观看日本欧美| 日韩视频精品在线| 精品91在线| 久久精品视频在线免费观看| 欧美成人免费在线视频| 久久精品99国产精品酒店日本| 亚洲黄页视频免费观看| 国产精品久久久久久久久久妞妞| 亚洲系列中文字幕| 亚洲免费影院| 日韩视频免费在线观看| 久久精品国语| 午夜精品久久99蜜桃的功能介绍| 99精品福利视频| 在线免费观看欧美| 亚洲欧美日本国产有色| 在线一区视频| 欧美va天堂| 亚洲免费精彩视频| 久久精品国产精品亚洲综合| 亚洲曰本av电影| 欧美日韩xxxxx| 欧美国产第一页| 永久免费视频成人| 久久精品国产精品亚洲精品| 欧美淫片网站| 国产日韩欧美综合精品| 久久五月天婷婷| 久久综合给合久久狠狠狠97色69| 精品动漫av| 欧美一级电影久久| 久久国内精品视频| 国产亚洲精品福利| 亚洲欧美日韩爽爽影院| 亚洲国产91| 9l视频自拍蝌蚪9l视频成人| 亚洲人成在线观看| 一本一本久久a久久精品综合妖精| 国产精品久久亚洲7777| 亚洲视频香蕉人妖| 小黄鸭视频精品导航| 国产精品海角社区在线观看| 一区二区三区**美女毛片| 一区二区免费看| 国产精品第三页| 午夜精品久久久久久久99樱桃| 在线观看91精品国产麻豆| 久久久美女艺术照精彩视频福利播放 | 一区二区三区欧美激情| 国产裸体写真av一区二区| 亚洲在线一区二区三区| 久久精品91久久香蕉加勒比| 激情丁香综合| 农夫在线精品视频免费观看| 91久久在线| 亚洲男人的天堂在线| 国产日韩欧美不卡| 最新日韩在线视频| 中日韩男男gay无套| 久久午夜av| 最新中文字幕一区二区三区| 国产美女精品免费电影| 久久成人免费日本黄色| 欧美激情亚洲国产| 亚洲午夜性刺激影院| 国产一区二区三区高清| 美日韩精品视频免费看| 99国产精品久久久久久久久久| 在线观看亚洲视频| 欧美日韩不卡| 羞羞漫画18久久大片| 亚洲第一精品电影| 亚洲欧美影音先锋| 亚洲国产一区二区在线| 国产精品免费视频xxxx| 亚洲永久免费精品| 亚洲激情视频在线播放| 欧美在线观看视频一区二区三区| 亚洲色图综合久久| 国产日韩精品在线观看| 欧美激情亚洲另类| 久久久xxx| 免费日韩av| 午夜精品久久久久久久久久久久| 欧美日韩aaaaa| 午夜激情综合网| 99精品欧美一区二区三区综合在线| 亚洲精品一区二区三区99| 久久全国免费视频| 久久亚洲影院| 亚洲自拍偷拍麻豆| 亚洲精品美女在线观看| 欧美fxxxxxx另类| 校园春色国产精品| 亚洲美女中文字幕| 欧美激情亚洲自拍| 亚洲美女在线观看| 激情综合电影网| 国产视频自拍一区| 国产精品一区二区在线观看| 亚洲专区在线视频| 日韩天堂av| 亚洲人成亚洲人成在线观看| 麻豆国产精品777777在线| 久久久久国产精品厨房| 久久精品99国产精品| 狠狠色狠狠色综合日日五| 欧美在线精品免播放器视频| 亚洲一区二区三区中文字幕| 99v久久综合狠狠综合久久| 91久久精品www人人做人人爽| 亚洲手机视频| 一本色道久久综合亚洲精品小说 | 老司机精品视频一区二区三区| 欧美成人一区二区在线| 久久久精品国产免大香伊| 永久域名在线精品| 欧美美女bbbb| 欧美激情亚洲综合一区| 欧美激情偷拍| 欧美日韩精品三区| 性欧美1819性猛交| 欧美在线电影| 久久中文字幕一区| 欧美成人性生活| 欧美精品三级日韩久久| 欧美午夜久久| 国产午夜亚洲精品理论片色戒| 猫咪成人在线观看| 欧美成人自拍视频| 欧美激情日韩| 国产精品乱子久久久久| 国产午夜精品一区理论片飘花| 欧美激情一区二区在线| 欧美91福利在线观看| 欧美中文字幕| 久久三级视频| 亚洲国产中文字幕在线观看| 91久久久久| 亚洲欧美视频一区| 久久中文欧美| 国产精品成人免费视频 | 国产精品亚洲综合一区在线观看| 一区二区免费看| 你懂的视频欧美| 亚洲欧洲另类| 亚洲一区二区黄色| 久久久亚洲国产美女国产盗摄| 一本久久a久久精品亚洲| 亚洲自拍偷拍色片视频| 亚洲黄色影院| 亚洲欧美激情一区二区| 欧美激情一区二区三级高清视频| 亚洲欧洲99久久| 蜜桃av久久久亚洲精品| 亚洲毛片网站| 久久久久9999亚洲精品| 欧美三区视频| 亚洲成色999久久网站| 亚洲私人影院在线观看| 日韩午夜黄色| 久久精品国产清高在天天线| 亚洲综合精品自拍| 麻豆精品精华液| 久热精品视频| 亚洲综合精品四区| 亚洲无人区一区| 美女黄色成人网| 国产亚洲第一区| 亚洲无线视频| 亚洲国产精品国自产拍av秋霞| 亚洲国产99精品国自产| 久久精彩视频|