n <- 5000
s <- cumsum((n - 1) : 1)
d <- dist((1:n) ^ 3)
r <- rle(sort(d[duplicated(d)]))
for (x in r$value[r$length >= 4]) {
i <- which(d == x)
b <- findInterval(i, s, left.open = TRUE) + 1
a <- i - s[b - 1] + b
cat
(x
, ": ", toString
(sprintf("(%d, %d)", a
, b
)), "\n", sep
= "")}
biA8LSA1MDAwCnMgPC0gY3Vtc3VtKChuIC0gMSkgOiAxKQpkIDwtIGRpc3QoKDE6bikgXiAzKQpyIDwtIHJsZShzb3J0KGRbZHVwbGljYXRlZChkKV0pKQoKZm9yICh4IGluIHIkdmFsdWVbciRsZW5ndGggPj0gNF0pIHsKICAgIGkgPC0gd2hpY2goZCA9PSB4KQogICAgYiA8LSBmaW5kSW50ZXJ2YWwoaSwgcywgbGVmdC5vcGVuID0gVFJVRSkgKyAxCiAgICBhIDwtIGkgLSBzW2IgLSAxXSArIGIKICAgIGNhdCh4LCAiOiAiLCB0b1N0cmluZyhzcHJpbnRmKCIoJWQsICVkKSIsIGEsIGIpKSwgIlxuIiwgc2VwID0gIiIpCn0=
1412774811: (1134, 357), (1155, 504), (1246, 805), (2115, 2004), (4746, 4725)