cp's OEIS Frontend

A007996 Primes that divide at least one term of Sylvester's sequence s = A000058: s(n+1) = s(n)^2 - s(n) + 1, s(0) = 2.

%I A007996 #54 Feb 16 2025 08:32:31
%S A007996 2,3,7,13,43,73,139,181,547,607,1033,1171,1459,1861,1987,2029,2287,
%T A007996 2437,4219,4519,6469,7603,8221,9829,12763,13147,13291,13999,15373,
%U A007996 17881,17977,19597,20161,20479,20641,20857,20929,21661,23689,23773,27031
%N A007996 Primes that divide at least one term of Sylvester's sequence s = A000058: s(n+1) = s(n)^2 - s(n) + 1, s(0) = 2.
%C A007996 Or, let S_1 = [2] and let S_{n+1} = list formed by sorting the union of S_n together with all prime factors of 1 + Product_i S_n(i) into increasing order; sequence is limit as n -> infinity of S_n.
%C A007996 Prime divisors of the terms of Sylvester's sequence A000058. - _Max Alekseyev_, Jan 03 2004. Also of A007018. - _N. J. A. Sloane_, Jan 27 2007
%C A007996 Because all terms of the sequence s(n) are coprime, a prime can divide at most one term. Odoni shows that primes p > 3 in this sequence must satisfy p = 1 (mod 6). - _T. D. Noe_, Sep 25 2010
%C A007996 See A180871(n) for the index of the first term of A000058 (this is one less than the index of the s-sequence) divisible by a(n). - _M. F. Hasler_, Apr 24 2014
%H A007996 Max Alekseyev, <a href="/A007996/b007996.txt">Table of n, a(n) for n = 1..12046</a> (first 8181 terms are also given at the Andersen link)
%H A007996 J. K. Andersen, <a href="http://primerecords.dk/sylvester-factors.htm">Factorization of Sylvester's sequence</a>
%H A007996 R. Mestrovic, <a href="http://arxiv.org/abs/1202.3670">Euclid's theorem on the infinitude of primes: a historical survey of its proofs (300 BC--2012) and another new proof</a>, arXiv preprint arXiv:1202.3670, 2012 - From N. J. A. Sloane, Jun 13 2012
%H A007996 R. W. K. Odoni, <a href="http://jlms.oxfordjournals.org/content/s2-32/1/1.extract">On the prime divisors of the sequence w_{n+1} = 1+w_1 ... w_n</a>, J. London Math. Soc. 32 (1985), 1-11.
%H A007996 Filip Saidak, <a href="http://www.jstor.org/stable/27642094">A New Proof of Euclid's Theorem</a>, Amer. Math. Monthly, Dec 2006
%H A007996 Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/SylvestersSequence.html">Sylvester's sequence</a>
%p A007996 n := 1; for p do if isprime(p) then x := 2 mod p; S := {}; while not member(x,S) do if x=0 then a[n] := p; n := n+1; break; fi; S := S union {x}; x := (x^2-x+1) mod p; od; fi; od;
%t A007996 t={}; p=1; While[Length[t]<100, p=NextPrime[p]; s=Mod[2,p]; k=0; modSet={}; While[s>0 && !MemberQ[modSet,s], AppendTo[modSet,s]; k++; s=Mod[s^2-s+1,p]]; If[s==0, AppendTo[t,{p,k}]]]; Transpose[t][[1]] (* _T. D. Noe_, Sep 25 2010 *)
%o A007996 (PARI) is(n)=my(k=Mod(2,n)); for(i=1, n, k=(k-1)*k+1; if(k==0, return(isprime(n)))); n==2 \\ _Charles R Greathouse IV_, Sep 30 2015
%Y A007996 The missing primes form A096264.
%Y A007996 Cf. A014546, A091335, A091336.
%Y A007996 Cf. A180871 (k such that a(n) divides A000058(k)).
%Y A007996 Cf. A323605 (smallest prime dividing A000058(n)).
%K A007996 nonn
%O A007996 1,1
%A A007996 Bennett Battaile (bennett.battaile(AT)autodesk.com)
%E A007996 More terms from _Max Alekseyev_, Jan 03 2004
%E A007996 Entry revised by _N. J. A. Sloane_, Jan 28 2007
%E A007996 Definition corrected (following a remark by _Don Reble_) by _M. F. Hasler_, Apr 24 2014