OFFSET
1,2
COMMENTS
An amazing Ramanujan identity. Here sigma_m(n) denotes Sum_{d|n} d^m.
REFERENCES
Bruce Berndt, Ramanujan's Notebooks Part II, Springer-Verlag; page 301.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = (1/240)*(sigma_5(2n+1)-sigma(2n+1)) (see A081863(2)).
MAPLE
f:= n -> 1/240*(numtheory:-sigma[5](2*n+1)-numtheory:-sigma(2*n+1)):
map(f, [$1..100]); # Robert Israel, Aug 12 2018
MATHEMATICA
lst={}; Do[AppendTo[lst, DivisorSigma[5, 2 n + 1] - DivisorSigma[1, 2 n + 1]], {n, 40}]; lst / 240 (* Vincenzo Librandi, Aug 13 2018 *)
Table[Sum[DivisorSigma[1, 2k+1]DivisorSigma[3, n-k], {k, 0, n-1}], {n, 35}] (* Harvey P. Dale, Jul 25 2020 *)
PROG
(PARI) a(n) = sum(k=0, n-1, sigma(2*k+1)*sigma(n-k, 3)); \\ Michel Marcus, Dec 04 2013
(PARI) a(n) = (sigma(2*n+1, 5) - sigma(2*n+1))/240; \\ Michel Marcus, Dec 04 2013
(Magma) [(DivisorSigma(5, 2*n+1)-DivisorSigma(1, 2*n+1))/240: n in [1..40]]; // Vincenzo Librandi, Aug 13 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Apr 11 2003
EXTENSIONS
Three more terms from Michel Marcus, Dec 04 2013
STATUS
approved