OFFSET
1,2
LINKS
Paolo Xausa, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = n^9 * Sum_{d|n, d<n, d odd} 1 / d^9.
G.f.: Sum_{k>=2} k^9 * x^k / (1 - x^(2*k)). - Ilya Gutkovskiy, May 19 2023
From Amiram Eldar, Oct 13 2023: (Start)
Sum_{k=1..n} a(k) = c * n^10 / 10, where c = 1023*zeta(10)/1024 = 1.0000170413... . (End)
EXAMPLE
a(10) = 10^9 * Sum_{d|10, d<10, d odd} 1 / d^9 = 10^9 * (1/1^9 + 1/5^9) = 1000000512.
MATHEMATICA
A352055[n_]:=DivisorSum[n, 1/#^9&, #<n&&OddQ[#]&]n^9; Array[A352055, 50] (* Paolo Xausa, Aug 10 2023 *)
a[n_] := DivisorSigma[-9, n/2^IntegerExponent[n, 2]] * n^9 - Mod[n, 2]; Array[a, 100] (* Amiram Eldar, Oct 13 2023 *)
PROG
(PARI) a(n) = n^9 * sigma(n >> valuation(n, 2), -9) - n % 2; \\ Amiram Eldar, Oct 13 2023
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Wesley Ivan Hurt, Mar 01 2022
STATUS
approved