login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A181970
Places of nonprimes in A050376.
2
3, 6, 9, 13, 20, 28, 37, 47, 63, 71, 83, 111, 127, 160, 177, 235, 280, 301, 348, 377, 430, 509, 542, 633, 700, 731, 838, 875, 915, 1030, 1194, 1284, 1327, 1415, 1458, 1559, 1752, 1915, 2015, 2181, 2231, 2531, 2590, 2773, 2960, 3089, 3154, 3289, 3485, 3562, 3919, 3997, 4142
OFFSET
1,1
COMMENTS
Or numbers not expressed in the form pi(p) + pi(sqrt(p)) + pi((sqrt(sqrt(p)))) +... with prime p.
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
FORMULA
a(n) ~ 0.5 n^2 log n. - Charles R Greathouse IV, Apr 11 2012
EXAMPLE
Show that 28 not expressed in form pi(p) + pi(sqrt(p)) + pi((sqrt(sqrt(p)))) +... with prime p. Indeed, for p=79, this sum is 22+4+1=27, while for p=83, it is 23+4+2=29.
PROG
(PARI) first_few(lim)=my(v=List(apply(n->n^2, primes(primepi(sqrtint(lim))))), u, t); forprime(p=2, (lim+.5)^(1/4), t=p^2; while((t=t^2)<=lim, listput(v, t))); listput(v, 1); v=vecsort(Vec(v)); u=vector(#v-1, i, sum(j=v[i]+1, v[i+1]-1, isprime(j))); u[1]++; for(i=2, #u, u[i]+=u[i-1]+1); u \\ Charles R Greathouse IV, Apr 10 2012
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Shevelev, Apr 06 2012
EXTENSIONS
a(30)-a(53) from Charles R Greathouse IV, Apr 10 2012
STATUS
approved