Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A122865
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A122865 Expansion of chi(q)* phi(q^3)* psi(-q^3) in powers of q where chi(), phi(), psi() are Ramanujan theta functions. +0
8
1, 1, 0, 2, 2, 1, 0, 0, 3, 0, 0, 2, 2, 2, 0, 0, 1, 2, 0, 2, 2, 1, 0, 0, 2, 0, 0, 2, 4, 0, 0, 0, 2, 3, 0, 2, 2, 0, 0, 0, 1, 0, 0, 4, 0, 2, 0, 0, 4, 2, 0, 0, 2, 2, 0, 0, 3, 0, 0, 2, 2, 0, 0, 0, 2, 1, 0, 2, 4, 2, 0, 0, 0, 0, 0, 2, 2, 2, 0, 0, 2, 2, 0, 4, 0, 1, 0, 0, 4, 0, 0, 2, 2, 0, 0, 0, 3, 2, 0, 2, 0, 0, 0, 0, 2 (list; graph; listen)
OFFSET

0,4

FORMULA

Expansion of q^(-1/3)*eta(q^2)^2*eta(q^6)^4/(eta(q)*eta(q^3)*eta(q^4)*eta(q^12)) in powers of q.

Euler transform of period 12 sequence [ 1, -1, 2, 0, 1, -4, 1, 0, 2, -1, 1, -2, ...]. - Michael Somos Apr 19 2007

PROGRAM

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

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

(PARI) {a(n)=if(n<0, 0, n=3*n+1; sumdiv(n, d, (d%4==1) -(d%4==3)))} /* Michael Somos Apr 19 2007 */

CROSSREFS

A002654(3n+1)=A035154(3n+1)=A113446(3n+1)=A122864(3n+1)=a(n).

Sequence in context: A029390 A108040 A137566 this_sequence A074080 A111407 A084440

Adjacent sequences: A122862 A122863 A122864 this_sequence A122866 A122867 A122868

KEYWORD

nonn

AUTHOR

Michael Somos, Sep 15 2006

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 November 25 08:46 EST 2009. Contains 167481 sequences.


AT&T Labs Research