login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A121454 Expansion of q * psi(-q) * psi(-q^7) in powers of q where psi(q) is a Ramanujan theta function. 2
1, -1, 0, -1, 0, 0, 1, -1, 1, 0, 2, 0, 0, -1, 0, -1, 0, -1, 0, 0, 0, -2, 2, 0, 1, 0, 0, -1, 2, 0, 0, -1, 0, 0, 0, -1, 2, 0, 0, 0, 0, 0, 2, -2, 0, -2, 0, 0, 1, -1, 0, 0, 2, 0, 0, -1, 0, -2, 0, 0, 0, 0, 1, -1, 0, 0, 2, 0, 0, 0, 2, -1, 0, -2, 0, 0, 2, 0, 2, 0, 1, 0, 0, 0, 0, -2, 0, -2, 0, 0, 0, -2, 0, 0, 0, 0, 0, -1, 2, -1, 0, 0, 0, 0, 0 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,11

COMMENTS

Ramanujan theta functions: f(q) := Product_{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) := Product_{k>=0} (1+q^(2k+1)) (A000700).

LINKS

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

Michael Somos, Introduction to Ramanujan theta functions

Eric Weisstein's World of Mathematics, Ramanujan Theta Functions

FORMULA

Expansion of eta(q)*eta(q^4)*eta(q^7)*eta(q^28)/(eta(q^2)*eta(q^14)) in powers of q.

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

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

Multiplicative with a(2^e) = -1 if e>0, a(7^e) = 1, a(p^e) = e+1 if p == 1, 2, 4 (mod 7), a(p^e) = (1+(-1)^e)/2 if p == 3, 5, 6 (mod 7).

G.f. A(x) satisfies 0=f(A(x), A(x^2), A(x^4)) where f(u, v, w)= u^4*w*v +2*u^3*w*v^2 +2*u^2*w^2*v^2 +4*u^3*w^3 +4*u^3*w^2*v +8*u*w^4*v +8*u*w^3*v^2 +8*u^2*w^3*v -u^2*v^4 -2*u^2*w*v^3 -4*w^2*v^4 -4*u*w^2*v^3.

G.f.: Sum_{k>0} (-1)^k *x^k*(1-x^(2*k))*(1-x^(4*k))*(1-x^(6*k))/(1-x^(14*k)) = x * Product_{k>0} (1-x^k)*(1+x^(2*k))*(1-x^(7*k))*(1+ x^(14*k)).

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

a(n) = (-1)^(n+1)*A035162(n).

MATHEMATICA

eta[q_]:= q^(1/24)*QPochhammer[q]; Rest[CoefficientList[Series[eta[q] eta[q^4] eta[q^7] eta[q^28]/(eta[q^2] eta[q^14]), {q, 0, 100}], q]] (* G. C. Greubel, Apr 19 2018 *)

PROG

(PARI) {a(n)=if(n<1, 0, -(-1)^n*sumdiv(n, d, kronecker(-28, d)))}

(PARI) {a(n)=local(A); if(n<1, 0, n--; A=x*O(x^n); polcoeff( eta(x+A)*eta(x^4+A)*eta(x^7+A)*eta(x^28+A)/eta(x^2+A)/eta(x^14+A), n))}

CROSSREFS

Cf. A035162.

Sequence in context: A087476 A307505 A035162 * A025462 A024879 A024316

Adjacent sequences:  A121451 A121452 A121453 * A121455 A121456 A121457

KEYWORD

sign,mult

AUTHOR

Michael Somos, Jul 30 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 October 21 14:42 EDT 2021. Contains 348155 sequences. (Running on oeis4.)