Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A113416
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A113416 Expansion of eta(q^2)^7*eta(q^4)/(eta(q)*eta(q^8)^2) in powers of q. +0
3
1, 2, -2, -4, -2, -8, 4, 16, -2, 14, 8, -20, 4, -24, -16, 16, -2, 36, -14, -36, 8, -32, 20, 48, 4, 42, 24, -40, -16, -56, -16, 64, -2, 40, -36, -64, -14, -72, 36, 48, 8, 84, 32, -84, 20, -56, -48, 96, 4, 114, -42, -72, 24, -104, 40, 80, -16, 72, 56, -116, -16, -120, -64, 112, -2, 96, -40, -132, -36, -96, 64, 144, -14 (list; graph; listen)
OFFSET

0,2

FORMULA

Euler transform of period 8 sequence [2, -5, 2, -6, 2, -5, 2, -4, ...].

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

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

PROGRAM

(PARI) a(n)=if(n<=0, n==0, -2*sumdiv(n, d, d*(d%2)*(-1)^(n/d+(d+1)\4)))

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

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

CROSSREFS

Sequence in context: A067538 A096154 A084540 this_sequence A131999 A103178 A087909

Adjacent sequences: A113413 A113414 A113415 this_sequence A113417 A113418 A113419

KEYWORD

sign

AUTHOR

Michael Somos, Oct 29 2005

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 23 17:09 EST 2009. Contains 167438 sequences.


AT&T Labs Research