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

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A261884 Expansion of (a(q) - a(q^2) - 2*a(q^3) + 2*a(q^6)) / 6 in powers of q where a() is a cubic AGM function. 1
1, -1, -1, 1, 0, 1, 2, -1, -1, 0, 0, -1, 2, -2, 0, 1, 0, 1, 2, 0, -2, 0, 0, 1, 1, -2, -1, 2, 0, 0, 2, -1, 0, 0, 0, -1, 2, -2, -2, 0, 0, 2, 2, 0, 0, 0, 0, -1, 3, -1, 0, 2, 0, 1, 0, -2, -2, 0, 0, 0, 2, -2, -2, 1, 0, 0, 2, 0, 0, 0, 0, 1, 2, -2, -1, 2, 0, 2, 2, 0 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,7

COMMENTS

Cubic AGM theta functions: a(q) (see A004016), b(q) (A005928), c(q) (A005882).

LINKS

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

FORMULA

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

a(n) is multiplicative with a(2^e) = (-1)^e, a(3^e) = -1 if e>0, a(p^e) = e+1 if p == 1 (mod 6), a(p^e) = (1 + (-1)^e)/2 if p == 5 (mod 6).

G.f.: Sum_{k>0} F(x^(6*k - 5)) - F(x^(6*k - 3)) + F(x^(6*k - 1)) where F(x) := x / (1 + x + x^2).

a(n) = A093829(n) unless n == 0 (mod 3). a(2*n) = - a(n). a(3*n + 1) = A033687(n).

a(6*n + 1) = A097195(n). a(6*n + 3) = - A033762(n). a(6*n + 5) = 0.

EXAMPLE

G.f. = x - x^2 - x^3 + x^4 + x^6 + 2*x^7 - x^8 - x^9 - x^12 + 2*x^13 + ...

MATHEMATICA

A004016[q_] := (QPochhammer[q]^3 + 9*q*QPochhammer[q^9]^3)/ QPochhammer[q^3]; A261884[n_] := SeriesCoefficient[(A004016[q] - A004016[q^2] - 2*A004016[q^3] + 2*A004016[q^6])/6, {q, 0, n}]; Table[A261884[n], {n, 1, 50}] (* G. C. Greubel, Sep 24 2017 *)

PROG

(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, (-1)^e, p==3, -1, p%6==1, e+1, 1-e%2)))};

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

CROSSREFS

Cf. A033687, A033762, A093829, A097195.

Sequence in context: A107110 A061197 A196346 * A035178 A093829 A113447

Adjacent sequences:  A261881 A261882 A261883 * A261885 A261886 A261887

KEYWORD

sign,mult

AUTHOR

Michael Somos, Sep 04 2015

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 22 00:51 EDT 2019. Contains 328315 sequences. (Running on oeis4.)