整点直角三角形

无心人

确实是37， 41有解

无心人

下面是51， 53， 61有解

无心人

63似乎无解

wayne

63是无解的，因为含有因子3,7，均是模4余3的

无心人

s(2)应该是1，本质上只有1组解

无心人

大的n情况下，用我原来的想法暴力搜索不行了
根据wayne思想，俺考虑如下算法

固定一个定点在原点，考虑在正整数坐标中，是否存在点A
使得$l_{AO} = n$，至少有两组平凡解的

然后，再确定第三点

无心人

用这个定义一下距离原点恰好是n的点集
let pn n = [(x, y)| x <- [0..n], y <- [0..n], x^2 + y^2 == n^2, x<=y]

无心人

测试了下，在Haskell中，计算pn 1000速度是可以接受的

无心人

完整的修改代码如下
let p n = [(x, y)|x<-[0..n], y<- [0..n], (x/=0)||(y/=0)]
let pn n =filter (\(x, y) ->(x <= y) && (x^2 + y^2 == n^2)) (p n)
let p3 n = [((0, 0), (x1, y1), (x2, y2))|(x1, y1)<- pn n, (x2, y2)<-p n, (x1, y1)/=(x2, y2), x2^2 + y2^2 + (x1-x2)^2 + (y1-y2)^2 == n^2]

wayne

63似乎无解
无心人 发表于 2010-6-10 14:56

63刚才被我漏掉了，现在补全了