OFFSET
1,1
LINKS
Chai Wah Wu, Table of n, a(n) for n = 1..10000
EXAMPLE
PROG
(PARI) { e=6; forprime(p=2, 1000, if(p%4==1, k=lift(sqrt(-1+O(p^e))); if(k>p^e/2, k=p^e-k); print1(k, ", "))) }
(Python)
from itertools import islice
from sympy import nextprime, sqrt_mod_iter
def A145299_gen(): # generator of terms
p = 1
while (p:=nextprime(p)):
if p&3==1:
yield min(sqrt_mod_iter(-1, p**6))
CROSSREFS
KEYWORD
nonn
AUTHOR
Klaus Brockhaus, Oct 17 2008
EXTENSIONS
More terms and efficient PARI program from. - Max Alekseyev, Oct 28 2008
STATUS
approved