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!)
A064939 a(n) = Sum_{i=1..omega(n)} i*p_i, where {p_i}, i=1..omega(n) is the increasing sequence of prime divisors of n, where omega is the number of distinct prime factors of n (A001221). 2
0, 2, 3, 2, 5, 8, 7, 2, 3, 12, 11, 8, 13, 16, 13, 2, 17, 8, 19, 12, 17, 24, 23, 8, 5, 28, 3, 16, 29, 23, 31, 2, 25, 36, 19, 8, 37, 40, 29, 12, 41, 29, 43, 24, 13, 48, 47, 8, 7, 12, 37, 28, 53, 8, 27, 16, 41, 60, 59, 23, 61, 64, 17, 2, 31, 41, 67, 36, 49, 33, 71, 8, 73, 76, 13, 40 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

LINKS

Harry J. Smith, Table of n, a(n) for n = 1..1000

FORMULA

a(n) = Sum_{p|n, q|n, p and q prime, p <= q} q. - Wesley Ivan Hurt, Aug 20 2020

EXAMPLE

factorset(30) = {2,3,5}, thus a(30) = 1*2 + 2*3 + 3*5 = 23.

MAPLE

with(numtheory): seq(add(i*sort(convert(factorset(n), 'list'))[i], i=1..nops(factorset(n))), n=1..200);

MATHEMATICA

ispd[n_]:=Module[{f=FactorInteger[n][[All, 1]]}, Total[f Range[ Length[f]]]]; Join[{0}, Array[ispd, 80, 2]] (* Harvey P. Dale, Aug 06 2017 *)

PROG

(PARI) { for (n=1, 1000, f=factor(n)~; a=sum(i=1, length(f), i*f[1, i]); write("b064939.txt", n, " ", a) ) } \\ Harry J. Smith, Sep 30 2009

CROSSREFS

Cf. A001221 (omega).

Sequence in context: A085818 A270709 A323382 * A248012 A151549 A086989

Adjacent sequences:  A064936 A064937 A064938 * A064940 A064941 A064942

KEYWORD

nonn

AUTHOR

Vladeta Jovovic, Oct 27 2001

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 June 17 05:25 EDT 2021. Contains 345080 sequences. (Running on oeis4.)