OFFSET
1,1
COMMENTS
All terms, except a(2) = 5, are of the form (2^p + 1)/3 - the Wagstaff primes A000979 = {3, 11, 43, 683, 2731, 43691, 174763, ...}.
Indices of prime Jacobsthal numbers are listed in A107036 = {3, 4, 5, 7, 11, 13, 17, 19, 23, 31, 43, 61, ...}.
For n > 1, A107036(n) = A000978(n) (numbers m such that (2^m + 1)/3 is prime). - Alexander Adamchuk, Oct 10 2006
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..23
MATHEMATICA
Select[Table[(2^n + (-1)^(n - 1))/3, {n, 200}], PrimeQ] (* Vladimir Joseph Stephan Orlovsky, Mar 29 2011 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved