|
MATHEMATICA
|
(* complete sequence *) tx = {}; Do[d = {i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20}; If[Total[d^3] == Total[d]^2, n = FromDigits[d]; AppendTo[tx, n]], {i1, 0, 9}, {i2, i1, 9}, {i3, i2, 9}, {i4, i3, 9}, {i5, i4, 9}, {i6, i5, 9}, {i7, i6, 9}, {i8, i7, 9}, {i9, i8, 9}, {i10, i9, 9}, {i11, i10, 9}, {i12, i11, 9}, {i13, i12, 9}, {i14, i13, 9}, {i15, i14, 9}, {i16, i15, 9}, {i17, i16, 9}, {i18, i17, 9}, {i19, i18, 9}, {i20, i19, 9}]; tx = Rest[tx]
(* partial sequence *) nddQ[n_] := Module[{idn=IntegerDigits[n]}, Min[Differences[idn]] >= 0 && Total[idn^3] == Total[idn]^2]; Select[Range[2000000], nddQ] (* Harvey P. Dale, Sep 01 2013 *)
|