login
Septic artiads: primes p congruent to 1 mod 14 for which all solutions of the congruence x^3 + x^2 - 2x - 1 == 0 (mod p) are 7th power residues.
7

%I #16 Mar 06 2020 15:47:03

%S 14197,21617,23801,24977,25999,34763,37549,41959,42407,45053,45599,

%T 54713,55987,56099,60271,61657,63463,66067,72577,75307,76343,76777,

%U 79283,83357,88397,90469,91309,99611,107927,111217,111301,111791,124699,126127,131251,132287

%N Septic artiads: primes p congruent to 1 mod 14 for which all solutions of the congruence x^3 + x^2 - 2x - 1 == 0 (mod p) are 7th power residues.

%H Eric M. Schmidt, <a href="/A270800/b270800.txt">Table of n, a(n) for n = 1..1000</a>

%H E. Lehmer, <a href="http://dx.doi.org/10.1016/0022-247X(66)90145-4">Artiads characterized</a>, J. Math. Anal. Appl. 15 1966 118-131. See page 126 (but beware errors).

%H E. Lehmer, <a href="/A001583/a001583b.pdf">Artiads characterized</a>, J. Math. Anal. Appl. 15 1966 118-131 [annotated and corrected scanned copy]

%o (Sage)

%o def is_septic_artiad(n) :

%o if not (n % 14 == 1 and is_prime(n)) : return False

%o R.<t> = PolynomialRing(GF(n))

%o return all(r[0]^((n-1)//7) == 1 for r in (t^3 + t^2 - 2*t - 1).roots())

%o # _Eric M. Schmidt_, Apr 02 2016

%Y Cf. A001583.

%K nonn

%O 1,1

%A _N. J. A. Sloane_, Apr 01 2016

%E Definition added and sequence extended and corrected by _Eric M. Schmidt_, Apr 02 2016