|
Search: id:A073918
|
|
|
| A073918 |
|
Smallest prime which is 1 more than a product of n distinct primes: a(n) is a prime and a(n) - 1 is a square-free number with n prime factors. |
|
+0 7
|
|
| 2, 3, 7, 31, 211, 2311, 43891, 870871, 13123111, 300690391, 6915878971, 200560490131, 11406069164491, 386480064480511, 18826412648012971, 693386350578511591, 37508276737897976011, 3087649419126112110271
(list; graph; listen)
|
|
|
OFFSET
|
0,1
|
|
|
COMMENT
|
Apparently the same as record values of A055734: least k such that phi(k) has n distinct prime factors, where phi is Euler's totient function. If the Mathematica program is used for large n, then "fact" should be reduced to, say, 1.1 in order to increase the search speed. - T. D. Noe (noe(AT)sspectra.com), Dec 17 2003
|
|
LINKS
|
M. F. Hasler (Maximilian.Hasler(AT)gmail.com) Jun 16 2007, Table of n, a(n) for n = 0..24
|
|
FORMULA
|
Theorem: For any m>0 there is K>0 such that for all k>K, a(k)-1 is divisible by the first m primes. - M. F. Hasler (Maximilian.Hasler(AT)gmail.com) Jun 16 2007
Corollary: For any m>1 there is K>0 such that for all k>K, a(k) = 1 (mod m). (M. F. Hasler, Jun 16 2007)
Let K(m) be the smallest possible K satisfying the Theorem. Conjecture: K(m) ~ m, i.e. a(k) ~ A002110(k), only very few of the last factors will be (insignificantly) larger. - M. F. Hasler (Maximilian.Hasler(AT)gmail.com) Jun 16 2007
|
|
EXAMPLE
|
a(6) = 43891 = 2*3*5*7*11*19+1.
a(0) = 1+1 = 2 (empty product of zero primes).
|
|
MATHEMATICA
|
Generate[pIndex_, i_] := Module[{p2, t}, p2=pIndex; While[p2[[i]]++; Do[p2[[j]]=p2[[i]]+j-i, {j, i+1, Length[p2]}]; t=Times@@Prime[p2]; t<fact*base, AppendTo[s, t]; If[i<Length[p2], Generate[p2, i+1]]]]; fact=2; Table[pin=Range[n]; base=Times@@Prime[pin]; s={base}; Do[Generate[pin, j], {j, n}]; s=Sort[s]; noPrime=True; i=0; While[noPrime&&i<Length[s], i++; noPrime=!PrimeQ[1+s[[i]]]]; If[noPrime, -1, 1+s[[i]]], {n, 20}] - from T. D. Noe
|
|
PROGRAM
|
(PARI) A073918(n, b=0 /*best*/, p=1 /*product*/, f=[]/*factors*/)={ if( #f<n, f=vector( n, i, prime(i)); p=prod( i=1, n-1, f[i] ); b=prime(n); /* get upper limit by incrementing last factor until prime is found */ while( !isprime( 1+p*b), b=nextprime(b+1)); b=1+p*b; p*=f[n] ); if( isprime( 1+p ), return( 1+p )); /* always p < b */ /* increase the n-th factor to recursively explore all solutions < b */ p /= f[n]; until( b <= 1+p*f[n] || ( n < #f && f[n] >= f[n+1] ) || !b = A073918( n-1, b, p*f[n], f), f[n]= nextprime( f[n]+1 ) ); b } /* then vector(30, n, A073918(n-1)) gives the first 30 terms */ - M. F. Hasler (Maximilian.Hasler(AT)gmail.com) Jun 16 2007
|
|
CROSSREFS
|
Cf. A055734 (number of distinct prime factors of phi(n)).
Cf. A073917, A098026.
Cf. A000040 (primes), A002110 (primorial), A081545 (same with composite instead of primes).
Sequence in context: A046972 A006862 A038710 this_sequence A096350 A018239 A066279
Adjacent sequences: A073915 A073916 A073917 this_sequence A073919 A073920 A073921
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
Amarnath Murthy (amarnath_murthy(AT)yahoo.com), Aug 18 2002
|
|
EXTENSIONS
|
More terms from Vladeta Jovovic (vladeta(AT)Eunet.yu), Aug 20 2002
|
|
|
Search completed in 0.002 seconds
|