OFFSET
1,2
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
FORMULA
Probably lim n -> infinity a(n)/n = 10.
MAPLE
filter:= n -> numtheory:-sigma(n)-n mod 10 = 0:
select(filter, [$1..1000]); # Robert Israel, Feb 28 2020
MATHEMATICA
Select[Range[600], Mod[#, 10]==Mod[DivisorSigma[1, #], 10]&] (* Harvey P. Dale, Mar 07 2024 *)
PROG
(PARI) isok(n) = (n % 10) == (sigma(n) % 10); \\ Michel Marcus, Nov 26 2013
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Benoit Cloitre, Aug 06 2002
STATUS
approved