Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A073918
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
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

page 1

Search completed in 0.002 seconds

Lookup | Welcome | Find friends | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
More pages | Superseeker | Maintained by N. J. A. Sloane (njas@research.att.com)

Last modified September 6 16:04 EDT 2008. Contains 143483 sequences.


AT&T Labs Research