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

O(1) 的小樂

Job Hunting

公告

記錄我的生活和工作。。。
<2025年11月>
2627282930311
2345678
9101112131415
16171819202122
23242526272829
30123456

統計

  • 隨筆 - 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>
            欧美午夜一区二区福利视频| 国产精品入口福利| 狠狠综合久久| 久久天堂成人| 老牛影视一区二区三区| 亚洲国产日韩欧美综合久久| 亚洲高清中文字幕| 欧美日韩国产一区二区三区地区| 一本色道久久综合狠狠躁的推荐| 99re66热这里只有精品3直播 | 91久久久久| 国产精品成人免费视频| 久久成人免费电影| 久久综合中文字幕| 国产精品99久久久久久宅男| 欧美亚洲视频在线看网址| 在线观看欧美激情| 日韩一级在线观看| 国产中文一区二区三区| 亚洲国产精品v| 国产精品任我爽爆在线播放| 美女网站在线免费欧美精品| 欧美日韩国产天堂| 久久精品免费观看| 欧美激情综合色综合啪啪| 欧美在线观看天堂一区二区三区| 久久久久久久性| 亚洲线精品一区二区三区八戒| 久久精品盗摄| 亚洲一区中文| 美女视频网站黄色亚洲| 亚洲欧美三级伦理| 欧美福利电影网| 久久国产毛片| 欧美午夜大胆人体| 欧美国产三级| 国产综合色在线| 中文亚洲欧美| 一本不卡影院| 免费观看国产成人| 久久露脸国产精品| 国产精品久久国产三级国电话系列| 欧美国产亚洲精品久久久8v| 国产午夜精品全部视频播放| 99综合电影在线视频| 亚洲经典在线| 久久精品免费观看| 久久精品国产综合| 国产精品免费aⅴ片在线观看| 亚洲第一精品影视| 极品尤物久久久av免费看| 亚洲午夜高清视频| 亚洲视频你懂的| 欧美激情导航| 亚洲人成啪啪网站| 最新亚洲电影| 欧美不卡视频一区| 欧美激情视频在线播放| 136国产福利精品导航| 久久成人免费日本黄色| 久久av资源网站| 国产午夜久久| 欧美影院在线| 老司机67194精品线观看| 国精品一区二区| 久久久久九九九| 免费欧美高清视频| 亚洲日本欧美| 欧美日韩免费高清一区色橹橹| 亚洲国产影院| 中文av一区特黄| 国产精品乱码人人做人人爱 | 欧美尤物巨大精品爽| 国产精品一区二区三区四区| 亚洲女人天堂av| 欧美在线视频观看| 激情另类综合| 欧美成人69av| 一级日韩一区在线观看| 亚洲欧美日韩国产中文在线| 国产美女搞久久| 久久精品亚洲精品| 欧美国内亚洲| 亚洲在线一区二区| 国产一区二区三区日韩| 久久综合电影一区| 91久久极品少妇xxxxⅹ软件| 亚洲五月婷婷| 国产精品私房写真福利视频| 久久精品国产第一区二区三区最新章节 | 亚洲精品孕妇| 国产精品激情电影| 久久精品国产免费| 亚洲国产一区二区在线| 亚洲欧美日本另类| 在线不卡中文字幕播放| 欧美日本中文字幕| 性欧美大战久久久久久久免费观看| 猫咪成人在线观看| 一区二区日韩欧美| 国产一区视频网站| 欧美精品少妇一区二区三区| 亚洲综合电影| 亚洲韩国精品一区| 久久精品日韩欧美| 99视频精品全部免费在线| 国产日韩欧美亚洲| 欧美日韩国产黄| 久久久久免费观看| 亚洲午夜激情在线| 亚洲国产精品久久久久秋霞蜜臀| 亚洲欧美中日韩| 亚洲精品视频免费观看| 国产女主播一区二区| 欧美好骚综合网| 久久久久国产一区二区三区| 一区二区三区免费看| 欧美国产视频一区二区| 久久高清一区| 午夜影院日韩| 亚洲视频一区二区免费在线观看| 在线观看国产精品淫| 国产欧美一区二区三区另类精品| 欧美韩日一区| 美女露胸一区二区三区| 欧美一区二区在线播放| 亚洲午夜视频在线观看| 99国产精品国产精品毛片| 欧美激情亚洲| 欧美国产大片| 免费欧美日韩| 免费永久网站黄欧美| 久久久999国产| 久久精品成人一区二区三区蜜臀| 亚洲免费一在线| 亚洲一区免费看| 亚洲专区一区| 亚洲欧美另类中文字幕| 制服丝袜激情欧洲亚洲| 99视频在线观看一区三区| 亚洲人久久久| 日韩视频免费观看高清在线视频 | 亚洲激情女人| 亚洲国产成人一区| 亚洲高清三级视频| 亚洲国产精品一区二区www| 在线观看日韩www视频免费 | 国产一区二区精品在线观看| 国产精品久久综合| 国产精品理论片| 国产乱码精品一区二区三区五月婷| 国产精品豆花视频| 国产精品视频999| 国产日韩专区在线| 精品成人在线视频| 亚洲国产精品久久久久秋霞蜜臀| 最新精品在线| 中文欧美字幕免费| 午夜精品亚洲一区二区三区嫩草| 欧美中文字幕第一页| 久久午夜电影| 91久久精品国产91久久性色| 亚洲精品永久免费| 亚洲欧美日韩国产综合| 久久―日本道色综合久久| 欧美成人免费播放| 国产精品久久久久久久免费软件| 国产亚洲精品激情久久| 亚洲国产成人在线| 亚洲网址在线| 久久亚洲综合色| 亚洲久久一区二区| 性欧美video另类hd性玩具| 美女尤物久久精品| 欧美性色视频在线| 黑人一区二区三区四区五区| 亚洲精品在线观看视频| 欧美诱惑福利视频| 亚洲激情自拍| 久久大逼视频| 欧美日韩日日骚| 一色屋精品亚洲香蕉网站| 中日韩男男gay无套| 可以看av的网站久久看| 日韩一区二区电影网| 久久久噜久噜久久综合| 欧美日韩一卡二卡| 1000部精品久久久久久久久| 亚洲字幕在线观看| 亚洲激情专区| 久久国产精品一区二区三区四区| 欧美激情一区二区三区 | 免费观看亚洲视频大全| 欧美图区在线视频| 亚洲第一中文字幕在线观看| 亚洲欧美日韩爽爽影院| 亚洲国产精品悠悠久久琪琪| 欧美一区二区三区四区视频| 欧美日韩在线观看一区二区| 亚洲高清在线播放|