

A073899


a(1) = 1; then kth prime prime(k) followed by prime(k) consecutive composite numbers not occurring earlier.


1



1, 2, 4, 6, 3, 8, 9, 10, 5, 12, 14, 15, 16, 18, 7, 20, 21, 22, 24, 25, 26, 27, 11, 28, 30, 32, 33, 34, 35, 36, 38, 39, 40, 42, 13, 44, 45, 46, 48, 49, 50, 51, 52, 54, 55, 56, 57, 58, 17, 60, 62, 63, 64, 65, 66, 68, 69, 70, 72, 74, 75, 76
(list;
graph;
refs;
listen;
history;
text;
internal format)



OFFSET

1,2


COMMENTS

Primes grow faster than composites. Question: For what (the smallest) value of m is a(m) prime and bigger than the previous term which is obviously composite? Answer: There is no such m. Proof: Assume a(m) is the nth prime and a(m1) < a(m). Checking manually gives n>10. Then a(m) < 2*n*log(n). The number of composite numbers appearing before a(m) is apparently the sum of the first n1 primes, which is bigger than (n1)^2. This means that a(m1) is definitely bigger than (n1)^2. Therefore we have a(m) < 2*n*log(n) as well as a(m1) > (n1)^2. Therefore a(m1) > a(m).


LINKS

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


MATHEMATICA

a = {1}; For[n = 1, n < 9, n++, AppendTo[a, Prime[n]]; For[j = 1, j < Prime[n] + 1, j++, i = 4; While[PrimeQ[i]  Length[Intersection[a, {i}]] == 1, i++ ]; AppendTo[a, i]]]; a


PROG

(PARI) { nonprim = listcreate(50000) ; for(n=2, 50000, if( !isprime(n), listput(nonprim, n)) ; ) ; print("1, ") ; k=2 ; indxn = 1 ; for (n = 2, 80, pr=prime(k1); print1(pr, ", ") ; for(i=1, pr, print1(nonprim[indxn], ", ") ; indxn++ ; ); print("") ; k++ ; ) }  R. J. Mathar, Mar 14 2006


CROSSREFS

Cf. A073900.
Sequence in context: A115316 A089088 A251622 * A232846 A101543 A073900
Adjacent sequences: A073896 A073897 A073898 * A073900 A073901 A073902


KEYWORD

nonn


AUTHOR

Amarnath Murthy, Aug 18 2002


EXTENSIONS

More terms from R. J. Mathar, Mar 14 2006
Edited by Stefan Steinerberger, Aug 13 2007


STATUS

approved



