login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A077454 a(n) = sigma_3(n^3)/sigma(n^3). 4
1, 39, 511, 2359, 12621, 19929, 101179, 149943, 368089, 492219, 1611831, 1205449, 4457701, 3945981, 6449331, 9588151, 22722609, 14355471, 44576623, 29772939, 51702469, 62861409, 141611691, 76620873, 196890121, 173850339, 268218727, 238681261, 574336533, 251523909 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

LINKS

Amiram Eldar, Table of n, a(n) for n = 1..10000

FORMULA

a(n) = A001158(n^3)/A000203(n^3).

Multiplicative with a(p^e) = (p^(6*e+2) + p^(3*e+1) + 1)/(p^2 + p + 1). - Amiram Eldar, Sep 09 2020

EXAMPLE

a(2) = sigma_3(2^3)/sigma(2^3) = 585/15 = 39.

MATHEMATICA

f[p_, e_] := (p^(6*e+2) + p^(3*e+1) + 1)/(p^2 + p + 1); a[n_] := Times @@ (f @@@ FactorInteger[n]); Array[a, 30] (* Amiram Eldar, Sep 09 2020 *)

PROG

(PARI) a(n)=sumdiv(n^3, d, d^3)/sigma(n^3)

(PARI) a(n) = my(f=factor(n^3)); sigma(f, 3)/sigma(f); \\ Michel Marcus, Sep 09 2020

CROSSREFS

Cf. A000203, A000578, A001158, A057660, A077455, A077456.

Sequence in context: A229517 A254871 A193072 * A142976 A200409 A341565

Adjacent sequences:  A077451 A077452 A077453 * A077455 A077456 A077457

KEYWORD

nonn,easy,mult

AUTHOR

Benoit Cloitre, Nov 30 2002

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified August 11 12:12 EDT 2022. Contains 356065 sequences. (Running on oeis4.)