OFFSET
1,1
EXAMPLE
Do(179) = 25665020 = 25236484 + 428536 = Do(178) + Do(46);
Do(184) = 27880600 = 15086400 + 12794200 = Do(150) + Do(142).
MATHEMATICA
(* This is just a recomputation of j values, given i values. *)
do[n_] := n*(3*n - 1)*(3*n - 2)/2;
triples = Reap[Module[{s, i, j, k, n, ijk}, s[i_] := Solve[j >= k > 0 && do[i] == do[j] + do[k], {j, k}, Integers]; For[n = 1, n <= Length[A053017], n++, i = A053017[[n]]; ijk = {i, j, k} /. s[i] // First; Print[ijk]; Sow[ijk]]]][[2, 1]];
A053018 = triples[[All, 2]] (* Jean-François Alcover, Feb 17 2015, updated Jul 09 2022 *)
CROSSREFS
KEYWORD
nice,nonn
AUTHOR
Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de), Feb 24 2000
EXTENSIONS
More terms from Jon E. Schoenfield, Aug 13 2007
a(27)-a(32) from Donovan Johnson, Aug 15 2010
STATUS
approved