OFFSET
1,2
COMMENTS
Except for 2, subsequence of nonprime numbers (A018252).
All odd terms are squares. - Michel Marcus, Dec 09 2014
Each term is either a square or twice a square. - Robert Israel, Dec 09 2014
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
8 is in the sequence because sigma(8)+2 = 15+2 = 17 (prime).
MAPLE
N:= 10^5: # to get all terms <= N
select(n -> isprime(numtheory:-sigma(n)+2), {seq(i^2, i=1..floor(sqrt(N))), seq(2*i^2, i=1..floor(sqrt(N/2)))} );
# if using Maple 11 or earlier, uncomment the next line
# sort(convert(%, list)); # Robert Israel, Dec 09 2014
MATHEMATICA
Select[Range[10000], PrimeQ[DivisorSigma[1, #] + 2] &] (* Vincenzo Librandi, Nov 14 2014 *)
PROG
(Magma) [n: n in[1..100000] | IsPrime(SumOfDivisors(n) + 2)]
(PARI) for(n=1, 10^5, if(isprime(sigma(n)+2), print1(n, ", "))) \\ Derek Orr, Nov 13 2014
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Jaroslav Krizek, Nov 13 2014
STATUS
approved