OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
f[n_] := f[n] = Prime[n] - n * DivisorSigma[0, n]; Select[Range[2, 6500], f[#+1] == f[#-1] &] (* Amiram Eldar, Apr 08 2024 *)
PROG
(PARI) f(n) = prime(n) - n * numdiv(n);
is(n) = n > 1 && f(n+1) == f(n-1); \\ Amiram Eldar, Apr 08 2024
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, Feb 24 2002
STATUS
approved