login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A121363 Expansion of q^(-1/4)(eta(q)*eta(q^6)*eta(q^9)/eta(q^3))^2/(eta(q^2)eta(q^18)) in powers of q. 5
1, -2, 0, 2, -2, 0, 3, -2, 0, 2, -2, 0, 1, -2, 0, 2, -4, 0, 2, 0, 0, 4, -2, 0, 2, -2, 0, 2, -2, 0, 1, -4, 0, 0, -2, 0, 4, -2, 0, 2, 0, 0, 3, -2, 0, 2, -4, 0, 2, -2, 0, 4, 0, 0, 0, -4, 0, 2, -2, 0, 2, -2, 0, 0, -2, 0, 4, -2, 0, 2, -2, 0, 3, -2, 0, 0, -4, 0, 2, -2, 0, 6, 0, 0, 2, 0, 0, 2, -2, 0, 1, -4, 0, 2, -4, 0, 0, -2, 0, 2, -2, 0, 2, 0, 0 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,2

COMMENTS

Ramanujan theta functions: f(q) := Prod_{k>=1} (1-(-q)^k) (see A121373), phi(q) := theta_3(q) := Sum_{k=-oo..oo} q^(k^2) (A000122), psi(q) := Sum_{k=0..oo} q^(k*(k+1)/2) (A010054), chi(q) := Prod_{k>=0} (1+q^(2k+1)) (A000700).

LINKS

G. C. Greubel, Table of n, a(n) for n = 0..1000

M. Somos, Introduction to Ramanujan theta functions

Eric Weisstein's World of Mathematics, Ramanujan Theta Functions

FORMULA

Euler transform of period 18 sequence [ -2, -1, 0, -1, -2, -1, -2, -1, -2, -1, -2, -1, -2, -1, 0, -1, -2, -2, ...].

a(n) = b(4n+1) where b(n) is multiplicative and b(2^e) = b(3^e) = 0^e, b(p^e) = (1+(-1)^e)/2 if p == 3 (mod 4), b(p^e) = e+1 if p == 1 (mod 12), b(p^e) = (e+1)(-1)^e if p == 5 (mod 12).

G.f.: Product_{k>0} (1+x^(3n))^2(1-x^n)(1-x^(9n))/((1+x^n)(1+x^9n)).

a(3n+2) = 0.

Expansion of phi(-q)*phi(-q^9)/chi(-q^3)^2 in powers of q where phi(),chi() are Ramanujan theta functions.

MATHEMATICA

QP = QPochhammer; s = (QP[q]*QP[q^6]*(QP[q^9]/QP[q^3]))^2/QP[q^2]/QP[q^18]+ O[q]^105; CoefficientList[s, q] (* Jean-Fran├žois Alcover, Nov 30 2015, adapted from PARI *)

PROG

(PARI) {a(n)=if(n<0, 0, n=4*n+1; dirmul(vector(n, k, kronecker(12, k)), vector(n, k, kronecker(-12, k)))[n])}

(PARI) {a(n)=local(A, p, e); if(n<0, 0, n=4*n+1; A=factor(n); prod(k=1, matsize(A)[1], if(p=A[k, 1], e=A[k, 2]; if(p<5, 0, if(p%4==3, (1+(-1)^e)/2, (e+1)*if(p%3==2, (-1)^e, 1)))))) }

(PARI) {a(n)=local(A); if(n<0, 0, A=x*O(x^n); polcoeff( (eta(x+A)*eta(x^6+A)*eta(x^9+A)/ eta(x^3+A))^2/ eta(x^2+A)/ eta(x^18+A), n))}

CROSSREFS

A002175(n) = a(3n). A121444(n) = -a(3n+1)/2.

Sequence in context: A308654 A143613 A208955 * A214258 A138021 A166065

Adjacent sequences:  A121360 A121361 A121362 * A121364 A121365 A121366

KEYWORD

sign

AUTHOR

Michael Somos, Jul 22 2006

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified June 25 17:47 EDT 2019. Contains 324353 sequences. (Running on oeis4.)