18/11/08 16:05:00.26 wKTjJ6Fa.net
いつもの顰蹙解w
今回はダンプリストではなくてRのスクリプト(HaskellやPythonは独学中w)
a^2+b^2の平方根が整数の組み合わせを考えればいいんだから、簡単にプログラムが組めた。
A=100
pita=NULL
for(a in 1:A){
B=floor(a^2/2-1/2)
for(b in a:B){
c=a^2+b^2
if(floor(sqrt(c)) == sqrt(c) ){
pita=rbind(re,c(a,b,sqrt(c)))
}
}
}
> pita[7,]
[1] 9 40 41
> pita[77,]
[1] 42 56 70
> pita[100,]
[1] 50 120 130