OFFSET
1,13
COMMENTS
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..10000
Michael Somos, Introduction to Ramanujan theta functions, 2019.
Eric Weisstein's World of Mathematics, Ramanujan Theta Functions.
FORMULA
From Michael Somos, May 05 2015: (Start)
Expansion of (phi(q) * phi(q^22) + phi(q^2) * phi(q^11)) / 2 - 1 in powers of q where phi() is a Ramanujan theta function.
a(n) is multiplicative with a(p^e) = 1 if p = 2 or 11, a(p^e) = e+1 if Kronecker(-22, p) = +1, a(p^e) = (1 + (-1)^e)/2 if Kronecker(-22, p) = -1.
G.f.: Sum_{k>0} x^k / (1 - x^k) * Kronecker(-22, k).
a(n) = A255647(n) unless n = 0. (End)
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Pi/sqrt(22) = 0.669789... . - Amiram Eldar, Nov 17 2023
EXAMPLE
G.f. = q + q^2 + q^4 + q^8 + q^9 + q^11 + 2*q^13 + q^16 + q^18 + 2*q^19 + ...
MATHEMATICA
a[ n_] := If[ n < 1, 0, DivisorSum[ n, KroneckerSymbol[ -22, #] &]]; (* Michael Somos, May 05 2015 *)
a[ n_] := If[ n < 1, 0, Sum[ KroneckerSymbol[ -22, d], { d, Divisors[ n]}]]; (* Michael Somos, May 05 2015 *)
a[ n_] := SeriesCoefficient[ (EllipticTheta[ 3, 0, q] EllipticTheta[ 3, 0, q^22] + EllipticTheta[ 3, 0, q^2] EllipticTheta[ 3, 0, q^11]) / 2 - 1, {q, 0, n}]; (* Michael Somos, May 05 2015 *)
PROG
(PARI) {a(n) = if( n<1, 0, sumdiv(n, d, kronecker( -22, d)))}; /* Michael Somos, May 05 2015 */
(PARI) {a(n) = if( n<1, 0, direuler(p=2, n, 1 / ((1 - X) * (1 - kronecker( -22, p) * X)))[n])}; /* Michael Somos, May 05 2015 */
(PARI) {a(n) = my(A, p, e); if( n<1, 0, A = factor(n); prod( k=1, matsize(A)[1], [p, e] = A[k, ]; if( p==2 || p==11, 1, kronecker( -22, p) == 1, e+1, 1-e%2)))}; /* Michael Somos, May 05 2015 */
CROSSREFS
KEYWORD
nonn,easy,mult
AUTHOR
STATUS
approved