login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A167604 A variant of Euclid-Mullin (A000945) : a(1)=2, a(n+1) is the least prime dividing [prod_{i in I} a(i) + prod_{i not in I} a(i)], minimized over all I \subset {1,..,n}. 3
2, 3, 5, 11, 37, 13, 7, 29, 17, 19, 43, 23, 47, 41, 53, 31, 61, 59, 67, 79, 83, 73, 97, 71, 101, 89, 103, 127, 107, 113, 137, 131, 139, 109, 149, 151, 163, 157, 167, 173, 193, 211, 179, 191, 181, 223, 199, 197, 233, 227, 229, 239, 241, 251, 257, 307, 281, 269, 271, 293 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

COMMENTS

By Euclid's argument, the a(i) are distinct.

One can ask whether all primes occur in this sequence.

LINKS

Table of n, a(n) for n=1..60.

Andrew R. Booker, A variant of the Euclid-Mullin sequence containing every prime, arXiv preprint arXiv:1605.08929, 2016.

FORMULA

For any n, we have Legendre symbol (-a(1)*a(2)*...*a(n-1) / a(n)) = 1. If p is the smallest prime such that (-a(1)*a(2)*...*a(n-1) / p) = 1, then a(n) >= p. Conjecture: For all n, a(n) = p. Note that if b is such that b^2 == -a(1)*a(2)*...*a(n-1) (mod p) and for some I, b == prod_{i in I} a(i) (mod p), then a(n) = p. Heuristically, I must exist for large enough n, since the number of possible subsets I is much larger than p. - Max Alekseyev, Nov 11 2009, May 20 2015

EXAMPLE

a(4)=11 which is the smallest prime dividing the 4 partitions 2+3*5=17, 3+2*5=13, 5+2*3=11, 1+2*3*5=31.

PROG

(PARI) { A167604() = my(a, A, p, b, q, z); a = []; A=1; while(1, p=2; while( kronecker(-A, p)!=1, p=nextprime(p+1) ); b=lift(sqrt(-A+O(p))); z=znprimroot(p); q=lift(prod(i=1, #a, Mod(1+x^znlog(Mod(a[i], p), z, p-1), 1-x^(p-1)) )); if( polcoeff(q, znlog(Mod(b, p), z, p-1), x)==0, error("conjecture failed")); a=concat(a, [p]); A*=p; print1(p, ", ") ) } /* Max Alekseyev, May 20 2015 */

CROSSREFS

A167605 lists such n that the first n terms of a(n) is a permutation of the first n primes.

A000945 is the original Euclid-Mullin sequence (where I is restricted to the empty set).

Sequence in context: A061935 A067078 A124561 * A065510 A006721 A111289

Adjacent sequences:  A167601 A167602 A167603 * A167605 A167606 A167607

KEYWORD

nonn

AUTHOR

Kok Seng Chua (chuakokseng(AT)hotmail.com), Nov 07 2009

EXTENSIONS

Edited and extended by Max Alekseyev, Nov 11 2009

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified October 17 13:57 EDT 2019. Contains 328113 sequences. (Running on oeis4.)