OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..5000
FORMULA
a(n) = 2*A064205(n) + 1.
EXAMPLE
sigma((17-1)/2) + tau((17-1)/2) = sigma(8) + tau(8) = 15 + 4 = 19; 19 is prime, so 17 is in the sequence.
MATHEMATICA
Select[Range[3, 10^7, 2], PrimeQ[DivisorSigma[1, #] + DivisorSigma[0, #]] &[(# - 1)/2] &] (* Michael De Vlieger, Apr 20 2016 *)
PROG
(Magma) [n: n in [3..1000000] | IsPrime(NumberOfDivisors((n-1) div 2) + SumOfDivisors((n-1) div 2)) and (n-1) mod 2 eq 0]
(PARI) isok(n) = isprime(sigma((n-1)/2) + numdiv((n-1)/2));
lista(nn) = forstep (n=3, nn, 2, if (isok(n), print1(n, ", "))); \\ Michel Marcus, Apr 19 2016
(PARI) is(n)=my(f=factor(n\2)); n>2 && isprime(sigma(f)+numdiv(f)) && isprime(n) \\ Charles R Greathouse IV, Apr 29 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
Jaroslav Krizek, Apr 19 2016
STATUS
approved