Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A138741
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A138741 Expansion of q^(-1/2) * eta(q)^3*eta(q^4)*eta(q^12)/(eta(q^2)^2*eta(q^3)) in powers of q. +0
2
1, 3, 2, 0, 1, 0, 2, 6, 2, 0, 0, 0, 3, 3, 2, 0, 0, 0, 2, 6, 2, 0, 2, 0, 1, 6, 2, 0, 0, 0, 2, 0, 4, 0, 0, 0, 2, 9, 0, 0, 1, 0, 4, 6, 2, 0, 0, 0, 2, 0, 2, 0, 0, 0, 2, 6, 2, 0, 2, 0, 1, 6, 4, 0, 0, 0, 0, 6, 2, 0, 0, 0, 4, 3, 2, 0, 2, 0, 2, 6, 0, 0, 0, 0, 3, 0, 2, 0, 0, 0, 2, 6, 4, 0, 0, 0, 2, 12, 2, 0, 0, 0, 4, 0, 0 (list; graph; listen)
OFFSET

0,2

FORMULA

Expansion of q^(-1/2) * (theta_2(q)^2 + 3 * theta_2(q^3)^2) / 4 in powers of q.

Expansion of phi(q) * psi(q) * psi(q^3) / phi(q^3) in powers of q where phi(), psi() are Ramanujan theta functions.

Euler transform of period 12 sequence [ 3, -4, 2, -2, 3, -2, 3, -2, 2, -4, 3, -2, ...].

Moebius transform is period 24 sequence [ 1, -1, 2, 0, 1, -2, -1, 0, -2, -1, -1, 0, 1, 1, 2, 0, 1, 2, -1, 0, -2, 1, -1, 0, ...].

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

G.f. is a period 1 Fourier series which satisfies f(-1 / (48 t)) = 6 (t/i) g(t) where q = exp(2 pi i t) and g() is g.f. for A132003.

a(6*n + 3) = a(6*n + 5) = 0.

EXAMPLE

q + 3*q^3 + 2*q^5 + q^9 + 2*q^13 + 6*q^15 + 2*q^17 + 3*q^25 + 3*q^27 + ...

PROGRAM

(PARI) {a(n) = if( n<0 | n%6==3 | n%6==5, 0, if( n%2, 3, 1) * sumdiv( (n \ if( n%2, 3, 1)) * 2 + 1, d, kronecker( -4, d) ) )}

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

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

CROSSREFS

A116604(n) = (-1)^n * a(n). A008441(n) = a(2*n). A002175(n) = a(6*n). 3 * A008441(n) = a(6*n + 1). 2 * A121444(n) = a(6*n + 2).

Sequence in context: A133209 A131290 A116604 this_sequence A079618 A117406 A008783

Adjacent sequences: A138738 A138739 A138740 this_sequence A138742 A138743 A138744

KEYWORD

nonn

AUTHOR

Michael Somos, Mar 27 2008

page 1

Search completed in 0.002 seconds

Lookup | Welcome | Find friends | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
More pages | Superseeker | Maintained by N. J. A. Sloane (njas@research.att.com)

Last modified July 24 12:00 EDT 2008. Contains 142294 sequences.


AT&T Labs Research