n^2+(n+1)^2为平方数

medie2005

请求出$10^100$内所有满足$n^2+(n+1)^2=c^2$的n.

无心人

是毕达哥拉斯三角形的三个边的解
$n, n + 1, c$
设$X = n, Y = n + 1, Z = c$
则有$X = M^2 - N^2, Y = 2MN, Z = M^2 + N^2$
同样有$Y = M^2 - N^2, X = 2MN, Z = M^2 + N^2$
即$M^2 - N^2 - 2MN = +- 1$
变换得到$(M - N)^2 -2N^2 = +- 1$
是佩尔方程
有通解
设$M - N = x, N = y$
则有
$x = [(1 + sqrt(2))^r + (1-sqrt(2))^r] / 2$
$y = [(1 + sqrt(2))^r - (1-sqrt(2))^r] / {2sqrt(2)}$
设
$a_0 = 1, b_0 = 1, a_n = a_{n-1} + 2 b_{n-1}, b_n = a_{n-1} + b_{n-1}$
则有
$x = a_r, y = b_r$
$M = x + y, N = y$
已经得到结果
以下省略
呵呵

medie2005

$X_{0}=1,Y_{0}=0$
$X_{k+1 }= 3*X_{k} + 4*Y_{k} + 2$
$Y_{k+1} = 2* X_{k} + 3*Y_{k} + 1$
那么，$Y_{k}(k>=1)$就是所求的n的所有解.

无心人

那就

    takeWhile (\(x, y) -> y <= 10^100)  \$ iterate (\(x, y) -> (3*x + 4*y + 2, 2*x + 3*y + 1))  (5, 3)

得到你的答案

zYr

楼上的数学公式是不是没写好？

无心人

随便写的
哪里错了？
请指正

zYr

没显示出来呀！

zYr

颜色不是蓝的啊
你没看出来么

gxqcn

估计不是数学公式，而是某个外星人用的程序代码。

无心人

呵呵
是计算这个问题的程序代码啊

和数学公式差不多吧