OFFSET
0,2
LINKS
Antti Karttunen, Table of n, a(n) for n = 0..10000
M. Kaneko and D. Zagier, Supersingular j-invariants, hypergeometric series and Atkin's orthogonal polynomials, pp. 97-126 of D. A. Buell and J. T. Teitelbaum, eds., Computational Perspectives on Number Theory, Amer. Math. Soc., 1998
FORMULA
a(0) = 1; for n >= 1, a(n) = 4*A000005(n). [After the PARI-program of Michael Somos.] - Antti Karttunen, May 25 2017
EXAMPLE
1 + 4*q + 8*q^2 + 8*q^3 + 12*q^4 + 8*q^5 + 16*q^6 + 8*q^7 + 16*q^8 + ...
MAPLE
with(numtheory); E:=proc(k) series(1-(2*k/bernoulli(k))*add( sigma[k-1](n)*q^n, n=1..60), q, 61); end; E(1);
MATHEMATICA
terms = 61; CoefficientList[1+4*Sum[x^k/(1-x^k), {k, 1, terms}]+O[x]^terms, x] (* Jean-François Alcover, Feb 27 2018 *)
PROG
(PARI) {a(n) = if( n<1, n==0, 4 * numdiv(n))} /* Michael Somos, Jul 04 2011 */
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Feb 28 2009
STATUS
approved