OFFSET
1,2
COMMENTS
Terms are odd.
The next prime term is 807281. - Michel Marcus, Jun 01 2024
LINKS
Michel Marcus, Table of n, a(n) for n = 1..74
MAPLE
select(n -> A373224(n) = 0, [seq(1..300)]);
PROG
(SageMath)
def is_a(n):
p = nth_prime(n)
r = (p - 1)//2
return 0 == sum((-1)^(r*((q - 1)//2)) for q in prime_range(p))
def aList(upto): return [n for n in range(1, upto, 2) if is_a(n)]
print(aList(222))
(PARI) KP(p, q) = kronecker(p, q);
isok(k) = my(vp=primes(k)); vecsum(vector(k, j, KP(vp[k], vp[j])*KP(vp[j], vp[k]))) == 0; \\ Michel Marcus, Jun 01 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Peter Luschny, May 31 2024
STATUS
approved