原です。

某K川さんからも、直接メールいただきました。

require 'narray'
Max = 1000 # real maximum of this program is 1290
a = NArray.int(Max,Max).indgen!
x = a / Max + 1
y = a % Max + 1
pair = x * x * x + y * y * y
sorted = pair.sort
index  = pair.sort_index

f =  ( sorted[0...-4].eq sorted[3...-1] ).where
for i in (0...f.size)
  printf("%d = ",sorted[f[i]])
  printf("%d ^ 3 + %d ^ 3 = %d ^3 + %d ^3\n",
	x[index[f[i]]]  ,y[index[f[i]]],
	x[index[f[i]+1]],y[index[f[i]+1]])
end


Max は 170 ぐらいで十分なのですが、とにかく暴速です。
NArray のパワー恐るべし。Ruby プログラムに見えないと
ころも面白い。