OFFSET
1,1
COMMENTS
A prime p is in this sequence iff all prime divisors of ord_p(2)/2 are in this sequence, where ord_p(2) is the order of 2 modulo p. - Max Alekseyev, Jul 30 2006
LINKS
Joerg Arndt, Table of n, a(n) for n = 1..220 (terms up to 10^9, terms for n = 1..100 from T. D. Noe)
Alexander Kalmynin, On Novák numbers, arXiv:1611.00417 [math.NT], 2016. See Chapter 4 p. 7 Novák primes.
C. Smyth, The terms in Lucas Sequences divisible by their indices, JIS 13 (2010) #10.2.4.
EXAMPLE
2^171 + 1 == 0 (mod 171), 171 = 3^2*19, 2^13203+1 == 0 (mod 13203), 13203 = 3^4*163.
MATHEMATICA
S = {2}; Reap[For[p = 3, p < 2 10^6, p = NextPrime[p], f = FactorInteger[ MultiplicativeOrder[2, p]]; If[f[[1, 1]] != 2 || f[[1, 2]] != 1, Continue[]]; f = f[[All, 1]]; If[Length[Intersection[S, f]] == Length[f], S = Union[S, {p}]; Print[p]; Sow[p]]]][[2, 1]] (* Jean-François Alcover, Nov 11 2018, from PARI *)
PROG
(PARI) { A057719() = local(S, f); S=Set([2]); forprime(p=3, 10^7, f=factorint(znorder(Mod(2, p))); if(f[1, 1]!=2||f[1, 2]!=1, next); f=f[, 1]; if(length(setintersect(S, Set(f)))==length(f), S=setunion(S, [p]); print1(p, ", "))) }
CROSSREFS
KEYWORD
nonn
AUTHOR
Ignacio Larrosa Cañestro, Oct 26 2000
EXTENSIONS
Edited by Max Alekseyev, Jul 30 2006
STATUS
approved