Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.
%I #15 Sep 08 2022 08:46:07
%S 1,0,10,-10,10,0,0,10,0,-10,10,0,10,-10,20,-10,0,10,-10,0,10,0,20,-10,
%T 0,0,0,0,10,0,0,0,10,-20,20,10,0,10,0,-20,0,0,20,-10,20,-10,0,10,-10,
%U 10,10,0,10,-10,0,0,0,0,0,0,10,0,20,-10,10,-10,0,10,10
%N Expansion of (f(-q^2, -q^3)^5 - 3 * q * f(-q, -q^4)^5) / f(-q)^3 in powers of q where f() is a Ramanujan theta function.
%C Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
%H G. C. Greubel, <a href="/A239051/b239051.txt">Table of n, a(n) for n = 0..1000</a>
%H Michael Somos, <a href="/A010815/a010815.txt">Introduction to Ramanujan theta functions</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/RamanujanThetaFunctions.html">Ramanujan Theta Functions</a>
%F Moebius transform is period 5 sequence [ 0, 10, -10, 0, 0, ...].
%F G.f.: 1 + 10 * ( Sum_{k>=0} x^(5*k + 2) / (1 - x^(5*k + 2)) - x^(5*k + 3) / (1 - x^(5*k + 3)) ).
%F a(n) = A227216(n) - 3 * A229802(n).
%F a(5*n) = a(n). a(5*n + 1) = 0.
%e G.f. = 1 + 10*q^2 - 10*q^3 + 10*q^4 + 10*q^7 - 10*q^9 + 10*q^10 + 10*q^12 + ...
%t a[ n_] := If[ n < 1, Boole[ n == 0], 10 Sum[ {0, 1, -1, 0, 0}[[ Mod[ d, 5, 1] ]] ], {d, Divisors @ n}]];
%o (PARI) {a(n) = if( n<1, n==0, 10 * sumdiv(n, d, (d%5==2) - (d%5==3)))};
%o (Sage) ModularForms( Gamma1(5), 1, prec=70).0;
%o (Magma) Basis( ModularForms( Gamma1(5), 1), 70) [1];
%Y Cf. A227216, A229802.
%K sign
%O 0,3
%A _Michael Somos_, Jun 13 2014