This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A088840 #21 Oct 03 2023 13:11:49 %S A088840 1,1,1,7,1,1,1,5,1,1,1,7,1,1,1,31,1,1,1,7,1,1,1,5,1,1,1,7,1,1,1,21,1, %T A088840 1,1,7,1,1,1,5,1,1,1,7,1,1,1,31,1,1,1,7,1,1,1,5,1,1,1,7,1,1,1,127,1,1, %U A088840 1,7,1,1,1,5,1,1,1,7,1,1,1,31,1,1,1,7,1,1,1,5,1,1,1,7,1,1,1,21,1,1,1,7,1,1 %N A088840 Denominator of sigma(4n)/sigma(n). %H A088840 Antti Karttunen, <a href="/A088840/b088840.txt">Table of n, a(n) for n = 1..16384</a> %H A088840 <a href="/index/Si#SIGMAN">Index entries for sequences related to sigma(n)</a>. %F A088840 From _Amiram Eldar_, Oct 03 2023: (Start) %F A088840 Multiplicative with a(2^e) = (((-1)^e+2)*(2^(e+1)-1))/3 = A213243(e+1), and a(p^e) = 1 for an odd prime p. %F A088840 a(n) = A213243(A007814(n+1)). %F A088840 Dirichlet g.f.: ((8^s + 4^s + 2^(s+1))/(8^s + 4^s - 2^(s+2) - 4)) * zeta(s). %F A088840 Sum_{k=1..n} a(k) = (2*n/(3*log(2))) * (log(n) + gamma - 1 + 7*log(2)/12), where gamma is Euler's constant (A001620). (End) %t A088840 Table[Denominator[DivisorSigma[1, 4*n]/DivisorSigma[1, n]], {n, 1, 128}] %t A088840 a[n_] := Module[{e = IntegerExponent[n, 2]}, (((-1)^e+2)*(2^(e+1)-1))/3]; Array[a, 100] (* _Amiram Eldar_, Oct 03 2023 *) %o A088840 (PARI) A088840(n) = denominator(sigma(4*n)/sigma(n)); \\ _Antti Karttunen_, Nov 18 2017 %o A088840 (PARI) a(n) = {my(e = valuation(n, 2)); (((-1)^e+2) * (2^(e+1)-1))/3;} \\ _Amiram Eldar_, Oct 03 2023 %Y A088840 See A088839 for numerator. %Y A088840 Cf. A088837, A088838, A088841, A038712, A080278, A000203, A001620, A007814, A193553, A213243. %K A088840 easy,nonn,mult,frac %O A088840 1,4 %A A088840 _Labos Elemer_, Nov 04 2003 %E A088840 Typo in definition corrected by _Antti Karttunen_, Nov 18 2017