OFFSET
1,2
COMMENTS
Sophie Germain primes can alternatively be Lucasian primes, primes of the form 4k + 1, or, the individual prime 2.
FORMULA
EXAMPLE
There are five Sophie Germain Primes of the form 4k + 1 below 10^2: {5, 29, 41, 53, 89}, therefore a(2) = 5.
MATHEMATICA
nonLucSophies = Select[4Range[2500000] + 1, PrimeQ[#] && PrimeQ[2# + 1] &]; Table[Length[Select[nonLucSophies, # < 10^n &]], {n, 0, 7}]
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Rodolfo Ruiz-Huidobro, Mar 27 2019
STATUS
approved