login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A333700 a(n) = Sum_{k=1..n} pi(k) * pi(n-k). 0
0, 0, 0, 1, 4, 8, 14, 22, 32, 45, 58, 73, 90, 110, 132, 158, 184, 214, 246, 282, 320, 363, 406, 455, 506, 562, 618, 678, 738, 804, 872, 944, 1018, 1099, 1180, 1269, 1358, 1450, 1544, 1644, 1744, 1852, 1962, 2078, 2196, 2321, 2446, 2581, 2718, 2863 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,5

COMMENTS

Convolution of A000720 with itself.

LINKS

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

Eric Weisstein's World of Mathematics, Prime Counting Function

FORMULA

G.f.: (1/(1 - x)^2) * (Sum_{k>=1} x^prime(k))^2.

a(n) = Sum_{k=1..n} A046992(k) * A010051(n-k).

a(n) = Sum_{k=1..n} k * A073610(n-k+1).

MATHEMATICA

Table[Sum[PrimePi[k] PrimePi[n - k], {k, n}], {n, 50}]

nmax = 50; CoefficientList[Series[(1/(1 - x)^2) Sum[x^Prime[k], {k, 1, nmax}]^2, {x, 0, nmax}], x] // Rest

PROG

(PARI) a(n) = sum(k=1, n, primepi(k)*primepi(n-k)); \\ Michel Marcus, Apr 03 2020

CROSSREFS

Cf. A000720, A002815, A010051, A034387, A046992, A073610.

Sequence in context: A054347 A194149 A003682 * A011897 A110895 A049628

Adjacent sequences:  A333697 A333698 A333699 * A333701 A333702 A333703

KEYWORD

nonn

AUTHOR

Ilya Gutkovskiy, Apr 02 2020

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 March 2 06:13 EST 2021. Contains 341742 sequences. (Running on oeis4.)