RosettaCodeData/Task/Pythagorean-quadruples/R/pythagorean-quadruples.r
2023-07-01 13:44:08 -04:00

5 lines
188 B
R

squares <- d <- seq_len(2200)^2
aAndb <- outer(squares, squares, '+')
aAndb <- aAndb[upper.tri(aAndb, diag = TRUE)]
sapply(squares, function(c) d <<- setdiff(d, aAndb + c))
print(sqrt(d))