OFFSET
1,2
EXAMPLE
The sum of the digits of 531^sigma(531) is 9558 and 9558 is divisible by 531, so 531 is in the sequence.
MATHEMATICA
Do[s = n^DivisorSigma[1, n]; k = Plus @@ IntegerDigits[s]; If[Mod[k, n] == 0, Print[n]], {n, 1, 10000}]
Select[Range[150000], Divisible[Total[IntegerDigits[#^DivisorSigma[ 1, #]]], #]&] (* Harvey P. Dale, Jul 19 2013 *)
CROSSREFS
KEYWORD
base,nonn,more
AUTHOR
Ryan Propper, Aug 06 2005
EXTENSIONS
More terms from Ryan Propper, Oct 10 2005
More terms from Harvey P. Dale, Jul 19 2013
STATUS
approved