OFFSET
1,2
LINKS
Giovanni Resta, Table of n, a(n) for n = 1..1000
EXAMPLE
sigma(541443) = 921600 = 5^2 * 4^2 * 1^2 *4^2* 4^2 *3^2, so 541443 is a term of the sequence.
MATHEMATICA
f[n_] := Module[{a, l}, a = IntegerDigits[n]; l = Length[a]; Product[a[[i]], {i, 1, l}]^2]; Do[If[f[n] == DivisorSigma[1, n], Print[n]], {n, 2, 10^6}]
Select[Range[25*10^6], DivisorSigma[1, #]==(Times@@IntegerDigits[#])^2&] (* Harvey P. Dale, Aug 29 2014 *)
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Joseph L. Pe, Mar 26 2002
EXTENSIONS
Added missing first term and a(14)-a(48) from Donovan Johnson, Nov 20 2010
STATUS
approved