A050463 a(n) = Sum_{d|n, n/d=1 mod 4} d^4.

1, 16, 81, 256, 626, 1296, 2401, 4096, 6562, 10016, 14641, 20736, 28562, 38416, 50706, 65536, 83522, 104992, 130321, 160256, 194482, 234256, 279841, 331776, 391251, 456992, 531522, 614656, 707282, 811296, 923521, 1048576, 1185922

Offset

1,2

Links

Seiichi Manyama, Table of n, a(n) for n = 1..10000

Formula

From Amiram Eldar, Nov 05 2023: (Start)

a(n) = A285989(n) - A050467(n).

a(n) = A050468(n) + A050467(n).

a(n) = (A050468(n) + A285989(n))/2.

Sum_{k=1..n} a(k) ~ c * n^5 / 5, where c = 5*Pi^5/3072 + 31*zeta(5)/64 = 1.000340795436113... . (End)

Mathematica

a[n_] := DivisorSum[n, #^4 &, Mod[n/#, 4] == 1 &]; Array[a, 50] (* Amiram Eldar, Jul 08 2023 *)

Prog

(PARI) a(n) = sumdiv(n, d, (n/d % 4 == 1) * d^4); \\ Amiram Eldar, Nov 05 2023

Crossrefs

Cf. A050448, A050467, A050468, A013663, A285989.

Cf. A050460, A050461, A050462.

Sequence in context: A050751 A014188 A352050 * A075578 A113316 A113317

Adjacent sequences: A050460 A050461 A050462 * A050464 A050465 A050466

Keyword

nonn,easy

Author

N. J. A. Sloane, Dec 23 1999

Extensions

Offset changed from 0 to 1 by Seiichi Manyama, Jul 08 2023

Status

approved