OFFSET
1,1
COMMENTS
An odd prime p is an optimus prime if (1 + sqrt(Legendre(-1, p)*p))^p - 1 = a + b*sqrt(Legendre(-1, p)*p), where gcd(a, b) = p.
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
Arkadii Slinko, Additive Representability of Finite Measurement Structures, 2007, 26 pp.
Arkadii Slinko, Additive Representability of Finite Measurement Structures, 2007, 26 pp. [Cached copy, permission requested]
PROG
(PARI) is(p)=if(p<3 || !isprime(p), return(0)); my(t=(2*quadgen(kronecker(-1, p)*p))^p); gcd(imag(t), real(t)-1)==p \\ Charles R Greathouse IV, Sep 26 2012
CROSSREFS
KEYWORD
nonn
AUTHOR
Arkadiusz Wesolowski, Sep 26 2012
STATUS
approved