OFFSET
0,1
LINKS
Michael S. Branicky, Table of n, a(n) for n = 0..2391
MATHEMATICA
a = {}; Do[k = 0; While[ !PrimeQ[k 4^n + (4^n - 1)/3], k++ ]; AppendTo[a, k], {n, 0, 50}]; a (*Artur Jasinski*)
lnk[n_]:=Module[{k=0, n4=4^n}, While[!PrimeQ[k*n4+(n4-1)/3], k++]; k]; Array[ lnk, 60, 0] (* Harvey P. Dale, May 28 2018 *)
PROG
(Python)
from sympy import isprime
def a(n):
k, fourn = 0, 4**n
while not isprime(k*fourn + (fourn-1)//3): k += 1
return k
print([a(n) for n in range(72)]) # Michael S. Branicky, May 18 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Artur Jasinski, Jan 19 2007
EXTENSIONS
Offset corrected and a(51) and beyond from Michael S. Branicky, May 18 2022
STATUS
approved