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!)
A373105 a(n) = sigma_10(n^2)/sigma_5(n^2). 4
1, 993, 58807, 1016801, 9762501, 58395351, 282458443, 1041204193, 3472494301, 9694163493, 25937263551, 59795016407, 137858120557, 280481233899, 574103396307, 1066193093601, 2015992480593, 3448186840893, 6131063781703, 9926520779301 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
LINKS
FORMULA
a(n) = Sum_{1 <= x_1, x_2, x_3, x_4, x_5 <= n} ( n/gcd(x_1, x_2, x_3, x_4, x^5, n) )^5.
a(n) = Sum_{d|n} mu(n/d) * (n/d)^5 * sigma_10(d).
From Amiram Eldar, May 25 2024: (Start)
Multiplicative with a(p^e) = (p^(10*e+5) + 1)/(p^5 + 1).
Dirichlet g.f.: zeta(s)*zeta(s-10)/zeta(s-5).
Sum_{k=1..n} a(k) ~ c * n^11 / 11, where c = zeta(11)/zeta(6) = 0.9834383562... . (End)
MATHEMATICA
f[p_, e_] := (p^(10*e+5) + 1)/(p^5 + 1); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 20] (* Amiram Eldar, May 25 2024 *)
PROG
(PARI) a(n) = sigma(n^2, 10)/sigma(n^2, 5);
(PARI) a(n) = sumdiv(n, d, moebius(n/d)*(n/d)^5*sigma(d, 10));
CROSSREFS
Sequence in context: A181707 A035857 A028513 * A267119 A210315 A182694
KEYWORD
nonn,mult
AUTHOR
Seiichi Manyama, May 25 2024
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 September 18 07:41 EDT 2024. Contains 375996 sequences. (Running on oeis4.)