OFFSET
1,1
COMMENTS
Corresponding values of sigma(p) of primes p from A158913.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
Max Alekseyev, PARI/GP Scripts for Miscellaneous Math Problems (invphi.gp).
FORMULA
a(n) = A158913(n) + 1.
MATHEMATICA
Union@ Select[ DivisorSigma[1, Select[ Range@ 100, !PrimeQ@# &]], PrimeQ[# - 1] &] (* Robert G. Wilson v, Feb 09 2012 *)
PROG
(PARI) list(lim) = forprime(p = 1, lim, if(invsigmaNum(p+1) > 1, print1(p+1, ", "))); \\ Amiram Eldar, Dec 16 2024, using Max Alekseyev's invphi.gp
CROSSREFS
KEYWORD
nonn
AUTHOR
Jaroslav Krizek, Feb 07 2012
STATUS
approved