OFFSET
1,1
COMMENTS
EXAMPLE
a(4)=33, because 3+3 is the perfect number 6.
MATHEMATICA
Select[Range[1500], MemberQ[{6, 28, 496, 8128}, Total[IntegerDigits[#]]]&] (* Harvey P. Dale, Mar 14 2023 *)
PROG
(R) perf=c(6, 28, 496, 8128, 33550336, 8589869056, 137438691328, 2305843008139952128); i=6; psod=c(); while(length(psod)<1000) {x=which(sumod(i)==perf); if(length(x)) psod=c(psod, i+x/10); i=i+1}
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Kevin L. Schwartz and Christian N. K. Anderson, Mar 23 2013
STATUS
approved