OFFSET
1,2
COMMENTS
a(15) > 10^5. - Robert Price, Mar 17 2015
All terms are odd, since 4 * 10^(2*k) - 9 = (2 * 10^k - 3)*(2 * 10^k + 3). - Robert Israel, Mar 17 2015
LINKS
FORMULA
a(n) = A101848(n) + 1.
EXAMPLE
n = 1, 17, 19 are members since 31, 399999999999999991 and 39999999999999999991 are primes.
MAPLE
select(n -> n::odd and isprime(4*10^n-9), [$1..10000]); # Robert Israel, Mar 17 2015
MATHEMATICA
Do[ If[ PrimeQ[4*10^n - 9], Print[n]], {n, 0, 10000}]
PROG
(Magma) [n: n in [1..300] |IsPrime(4*10^n - 9)]; // Vincenzo Librandi, Mar 18 2015
(PARI) is(n)=ispseudoprime(4*10^n-9) \\ Charles R Greathouse IV, Jun 12 2017
CROSSREFS
KEYWORD
more,nonn
AUTHOR
Julien Peter Benney (jpbenney(AT)ftml.net), Jan 15 2005
EXTENSIONS
a(11) from Kamada link by Herman Jamke (hermanjamke(AT)fastmail.fm), Jan 02 2008
a(12)-a(14) from Kamada data by Robert Price, Mar 17 2015
STATUS
approved