OFFSET
1,1
COMMENTS
All terms are either primes or powers of 2.
EXAMPLE
a(3) = 4 is a term because A002533(4) = 73 is prime.
MAPLE
B[0]:= 1: B[1]:= 1: R:= NULL: count:= 0:
for n from 2 while count < 18 do
B[n]:= 2*B[n-1]+5*B[n-2];
if (isprime(n) or n = 2^padic:-ordp(n, 2)) and isprime(A[n]) then
R:= R, n; count:= count+1
fi
od:
R;
MATHEMATICA
Flatten[Position[LinearRecurrence[{2, 5}, {1, 1}, 1000] , _Integer?PrimeQ]-1] (* James C. McMahon, May 06 2024 *)
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Robert Israel, May 03 2024
EXTENSIONS
a(21) from Michael S. Branicky, May 04 2024
STATUS
approved