OFFSET
1,2
COMMENTS
a(40) > 2.5*10^5. - Robert Price, Oct 17 2015
Since 4^(6k) + 25 = 4096^k + 25 == (1^k + 25) mod 13 = 26 mod 13 == 0 mod 13, no multiple of 6 will be in this sequence. - Timothy L. Tiffin, Aug 07 2016
FORMULA
a(n) = A157006(n)/2. - Robert Price, Oct 17 2015
EXAMPLE
From Timothy L. Tiffin, Aug 09 2016: (Start)
a(1) = 1, since 4^1 + 25 = 4 + 25 = 29, which is prime.
a(2) = 2, since 4^2 + 25 = 16 + 25 = 41, which is prime.
a(3) = 3, since 4^3 + 25 = 64 + 25 = 89, which is prime.
a(4) = 4, since 4^4 + 25 = 256 + 25 = 281, which is prime.
a(5) = 5, since 4^5 + 25 = 1024 + 25 = 1049, which is prime.
a(6) = 10, since 4^10 + 25 = 1048576 + 25 = 1048601, which is prime. (End)
MATHEMATICA
Select[Range[2000], PrimeQ[4^# + 25] &] (* T. D. Noe, Feb 03 2012 *)
PROG
(PARI) for(n=1, 1e5, if(isprime(4^n + 25), print1(n", "))) \\ Altug Alkan, Oct 17 2015
(Magma) [n: n in [1..1000] | IsPrime(4^n+25)]; // Vincenzo Librandi, Aug 07 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Jan 15 2012
EXTENSIONS
a(22)-a(39) derived from A157006 by Robert Price, Oct 17 2015
a(40)-a(48) derived from A157006 by Elmo R. Oliveira, Nov 28 2023
STATUS
approved