login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A304066 a(n) = Sum_{k=1..n} k*floor(n/prime(k)). 0
0, 1, 3, 4, 7, 10, 14, 15, 17, 21, 26, 29, 35, 40, 45, 46, 53, 56, 64, 68, 74, 80, 89, 92, 95, 102, 104, 109, 119, 125, 136, 137, 144, 152, 159, 162, 174, 183, 191, 195, 208, 215, 229, 235, 240, 250, 265, 268, 272, 276, 285, 292, 308, 311, 319, 324, 334, 345, 362, 368, 386, 398, 404, 405, 414 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,3

COMMENTS

Partial sums of A066328.

LINKS

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

Index entries for sequences computed from indices in prime factorization

FORMULA

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

a(p^k) = a(p^k-1) + pi(p), where p is a prime and pi() = A000720.

a(n) = A056239(A048803(n)).

MAPLE

seq(add(k*floor(n/ithprime(k)), k=1..n), n=1..65); # Paolo P. Lava, May 14 2018

MATHEMATICA

Table[Sum[k Floor[n/Prime[k]], {k, n}], {n, 65}]

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

a[n_] := Plus @@ (PrimePi[#[[1]]] & /@ FactorInteger[n]); a[1] = 0; Accumulate[Table[a[n], {n, 65}]]

CROSSREFS

Cf. A000040, A000720, A008472, A013939, A024916, A024924, A048803, A056239, A066328, A081401, A304038.

Sequence in context: A108579 A287410 A228643 * A296994 A236475 A287458

Adjacent sequences:  A304063 A304064 A304065 * A304067 A304068 A304069

KEYWORD

nonn

AUTHOR

Ilya Gutkovskiy, May 05 2018

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 20 23:22 EDT 2019. Contains 321354 sequences. (Running on oeis4.)