OFFSET
1,1
COMMENTS
All the prime numbers (A000040) are in this sequence. But there are other numbers (A023885). - Daniel Lignon, Feb 26 2015
LINKS
Daniel Lignon, Table of n, a(n) for n = 1..1000
MAPLE
q:=n->(irem(numtheory[sigma](n)-n, numtheory[tau](n)-1)=0):
select(q, [$2..159])[]; # Lorenzo Sauras Altuzarra, Feb 01 2023
MATHEMATICA
Select[Range[2000], IntegerQ[Mean[Most[Divisors[#]]]]&] (* Daniel Lignon, Feb 26 2015 *)
PROG
(PARI) isok(n) = (n != 1) && !((sigma(n) - n) % (numdiv(n)-1)); \\ Michel Marcus, Mar 01 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved