OFFSET
1,1
LINKS
Donovan Johnson, Table of n, a(n) for n = 1..1000
Eric W. Weisstein, Smith Number
EXAMPLE
a(4) = 141176320 = 2^9*5*55147;
3*1^6+2^6+3^6+4^6+6^6+7^6 = 1^6+9*2^6+4^6+3*5^6+7^6 = 169197
MATHEMATICA
fQ[n_] := Block[{id = Sort@ IntegerDigits@ n, fid = Sort@ Flatten[ IntegerDigits@ Table[ #[[1]], {#[[2]]}] & /@ FactorInteger@ n]}, While[ id[[1]] == 0, id = Drop[id, 1]]; While[ fid[[1]] == 0, fid = Drop[fid, 1]]; id != fid && Plus @@ (id^6) == Plus @@ (fid^6)]; k = 2; lst = {}; While[k < 50000001, If[fQ@ k, AppendTo[ lst, k]; Print@ k]; k++]; lst
CROSSREFS
KEYWORD
nonn,base,changed
AUTHOR
Paul Weisenhorn, Dec 19 2010
EXTENSIONS
Example corrected by Donovan Johnson, Jan 02 2013
STATUS
approved