Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: author:somos
Displaying 1-10 of 988 results found. page 1 2 3 4 5 6 7 8 9 10 11 ... 99
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A074664 Number of algebraically independent elements of degree n in the algebra of symmetric polynomials in noncommuting variables. +0
29
1, 1, 2, 6, 22, 92, 426, 2146, 11624, 67146, 411142, 2656052, 18035178, 128318314, 954086192, 7396278762, 59659032142, 499778527628, 4341025729290, 39035256389026, 362878164902216, 3482882959111530, 34472032118214598 (list; graph; listen)
OFFSET

1,3

COMMENT

Also the number of irreducible set partitions of size n (see A055105) {1}; {1,2}; {1,2,3}, {1,23}; ...; and also the number of set partitions of n which do not have a proper subset of parts with a union equal to a subset {1,2,...,j} with j<n (atomic set partitions, see A087903) {1}; {12}; {13,2}, {123}; ...

REFERENCES

N. Bergeron, C. Reutenauer, M. Rosas and M. Zabrocki, Invariants and Coinvariants of the Symmetric Group in Noncommuting Variables, arXiv:math.CO/0502082

D. E. Knuth, TAOCP, Vol. 4, Section 7.2.1.7, Problem 26.

M. C. Wolf, Symmetric Functions of Non-commutative Elements, Duke Math. J., 2 (1936), 626-637.

LINKS

T. D. Noe, Table of n, a(n) for n=1..100

M. Klazar, Bell numbers, their relatives and algebraic differential equations

FORMULA

G.f.: 1-1/B(x) where B(x) = g.f. for A000110 the Bell numbers.

a(n) = Sum_{k = 1, ..., n-1}A087903(n, k). a(n+1) = Sum{k = 0..n} A086329(n, k) . a(n+2) = Sum_{k = 0..n} A086211(n, k) . - DELEHAM Philippe (kolotoko(AT)wanadoo.fr), Jun 13 2004

G.f. x/(1-(x-x^2)/(1-x-(x-2x^2)/(1-2x-(x-3x^2)/...))) (a continued fraction). - Michael Somos Sep 22 2005

Hankel transform is A000142 . - Philippe DELEHAM (kolotoko(AT)wanadoo.fr), Jun 21 2007

EXAMPLE

m{1} = x1+x2+x3+..., so a(1) = 1

m{1,2} = x1 x2+x2 x1+x2 x3+x3 x2+x1 x3+..., m{12} = x1 x1+x2 x2+x3 x3+... where m{1} m{1} = m{1,2} + m{12}, so a(2)=2-1=1

m{1,2,3} = x1 x2 x3+x1 x2 x4+x1 x3 x4+..., m{12,3} = x1 x1 x2+x2 x2 x1+..., m{13,2} = x1 x2 x1+x2 x1 x2+..., m{1,23} = x1 x2 x2+x2 x1 x1+..., m{123}=x1 x1 x1+x2 x2 x2+... and there are 3 independent relations among these 5 elements m{12} m{1} = m{123} + m{12,3}, m{1} m{12} = m{123}+m{1,23}, m{1} m{1,1} = m{1,2,3}+m{12,3}+m{13,2} so a(3)=5-3=2

PROGRAM

(PARI) a(n)=if(n<0, 0, polcoeff(1-1/serlaplace(exp(exp(x+x*O(x^n))-1)), n))

CROSSREFS

Row sums of A055105, A055106, A055107. Cf. A098742, A003319.

Row sums of A087903, A055105, A055106, A055107

Sequence in context: A014330 A124294 A124295 this_sequence A091768 A150274 A109317

Adjacent sequences: A074661 A074662 A074663 this_sequence A074665 A074666 A074667

KEYWORD

nonn,easy,nice

AUTHOR

Michael Somos

EXTENSIONS

Edited by Mike Zabrocki (zabrocki(AT)mathstat.yorku.ca), Sep 03 2005

A088218 Total number of leaves in all rooted ordered trees with n edges. +0
28
1, 1, 3, 10, 35, 126, 462, 1716, 6435, 24310, 92378, 352716, 1352078, 5200300, 20058300, 77558760, 300540195, 1166803110, 4537567650, 17672631900, 68923264410, 269128937220, 1052049481860, 4116715363800, 16123801841550 (list; graph; listen)
OFFSET

0,3

COMMENT

Number of ordered partitions of n into n parts, allowing zeros (cf. A097070) is binomial(2*n-1,n) = a(n) = essentially A001700. - Vladeta Jovovic (vladeta(AT)eunet.rs), Sep 15 2004

a(n) = A110556(n)*(-1)^n, central terms in triangle A110555. - Reinhard Zumkeller (reinhard.zumkeller(AT)gmail.com), Jul 27 2005

Hankel transform is A000027; example: Det([1,1,3,10;1,3,10,35;3,10,35,126;10,35,126,462])=4 . - Philippe DELEHAM (kolotoko(AT)wanadoo.fr), Apr 13 2007

a(n) is the number of functions f:[n]->[n] such that for all x,y in [n] if x<y then f(x)<=f(y). So 2*a(n)-n=A045992(n) [From Geoffrey Critzer (critzer.geoffrey(AT)usd443.org), Apr 02 2009]

Hankel transform of the aeration of this sequence is A000027 doubled: 1,1,2,2,3,3,... [From Paul Barry (pbarry(AT)wit.ie), Sep 26 2009]

REFERENCES

Paul Barry, On Integer-Sequence-Based Constructions of Generalized Pascal Triangles, Journal of Integer Sequences, Vol. 9 (2006), Article 06.2.4.

FORMULA

a(n)=(0^n+C(2n, n))/2. - Paul Barry (pbarry(AT)wit.ie), May 21 2004

a(n) is the coefficient of x^n in 1/(1-x)^n and also the sum of the first n coefficients of 1/(1-x)^n. Given B(x) with the property that the coefficient of x^n in B(x)^n equals the sum of the first n coefficients of B(x)^n, then B(x)=B(0)/(1-x).

G.f.: 1/(2-C(x)) where C(x) is g.f. for Catalan numbers A000108.

G.f.: (1+1/sqrt(1-4x))/2. a(n)=binomial(2n-1,n).

a(n)=sum{k=0..n, binomial(2n, k)cos((n-k)*pi)}; a(n)=sum{k=0..n, binomial(n, (n-k)/2)(1+(-1)^(n-k))cos(k*pi/2)/2} (with interpolated zeros); a(n)=sum{k=0..floor(n/2), binomial(n, k)cos((n-2k)pi/2)} (with interpolated zeros); - Paul Barry (pbarry(AT)wit.ie), Nov 02 2004

a(n)=Sum_{k, 0<=k<=n}A094527(n,k)*(-1)^k . - Philippe DELEHAM (kolotoko(AT)wanadoo.fr), Mar 14 2007

EXAMPLE

The five rooted ordered trees with 3 edges have 10 leaves.

..x........................

..o..x.x..x......x.........

..o...o...o.x..x.o..x.x.x..

..r...r....r....r.....r....

MAPLE

seq(abs(binomial(-n, -2*n)), n=0..24); - Zerinvary Lajos (zerinvarylajos(AT)yahoo.com), Oct 03 2007

PROGRAM

(PARI) a(n)=sum(i=0, n, binomial(n+i-2, i))

(PARI) a(n)=if(n<0, 0, polcoeff((1+1/sqrt(1-4*x+x*O(x^n)))/2, n))

(PARI) a(n)=if(n<0, 0, polcoeff(1/(1-x+x*O(x^n))^n, n))

(PARI) a(n)=if(n<0, 0, binomial(2*n-1, n))

(PARI) {a(n)=if(n<1, n==0, polcoeff( subst((1-x)/(1-2*x), x, serreverse(x-x^2+x*O(x^n))), n))}

CROSSREFS

A001700(n)=a(n+1). a(n)=A024718(n)-A024718(n-1).

Sequence in context: A099908 A167403 A001700 this_sequence A110556 A072266 A085282

Adjacent sequences: A088215 A088216 A088217 this_sequence A088219 A088220 A088221

KEYWORD

nonn

AUTHOR

Michael Somos, Sep 24 2003

EXTENSIONS

Essentially the same as A001700, which has much more information.

A055615 a(n)=n*moebius(n) (cf. A008683). +0
18
1, -2, -3, 0, -5, 6, -7, 0, 0, 10, -11, 0, -13, 14, 15, 0, -17, 0, -19, 0, 21, 22, -23, 0, 0, 26, 0, 0, -29, -30, -31, 0, 33, 34, 35, 0, -37, 38, 39, 0, -41, -42, -43, 0, 0, 46, -47, 0, 0, 0, 51, 0, -53, 0, 55, 0, 57, 58, -59, 0, -61, 62, 0, 0, 65, -66, -67, 0, 69, -70, -71, 0 (list; graph; listen)
OFFSET

1,2

COMMENT

Dirichlet inverse of n.

Absolute values give n if n is square-free otherwise 0.

LINKS

T. D. Noe, Table of n, a(n) for n=1..1000

FORMULA

Dirichlet g.f.: 1/zeta(s-1).

Multiplicative with a(p^e) = -p*0^(e-1), e>0 and p prime. - Reinhard Zumkeller (reinhard.zumkeller(AT)gmail.com), Jul 17 2003

Conjectures: lim b->1+ Sum n=1..inf a(n)*b^(-n) = -12 and lim b->1- Sum n=1..inf a(n)*b^n = -12 (+ indicates that b decreases to 1, - indicates it increases to 1), both considering that zeta(-1) = -1/12 and calculations (more generally mu(n)*n^s is Abel summable to zeta(-s)). - Gerald McGarvey (Gerald.McGarvey(AT)comcast.net), Sep 26 2004

Dirichlet generating function for the absolute value: zeta(s-1)/zeta(2s-2). - Franklin T. Adams-Watters, Sep 11 2005.

PROGRAM

(PARI) a(n)=n*moebius(n)

(PARI) a(n)=if(n<1, 0, direuler(p=2, n, 1-p*X)[n])

CROSSREFS

Cf. A000027, A023900. Moebius transform of A023900.

Cf. A008683, A062004.

Sequence in context: A128214 A145105 A140700 this_sequence A049268 A004179 A122830

Adjacent sequences: A055612 A055613 A055614 this_sequence A055616 A055617 A055618

KEYWORD

sign,easy,nice,mult

AUTHOR

Michael Somos, Jun 04 2000

A073189 Integers 0..n three times then 0..n+1 three times etc. +0
16
0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 1, 2, 0, 1, 2, 0, 1, 2, 0, 1, 2, 3, 0, 1, 2, 3, 0, 1, 2, 3, 0, 1, 2, 3, 4, 0, 1, 2, 3, 4, 0, 1, 2, 3, 4, 0, 1, 2, 3, 4, 5, 0, 1, 2, 3, 4, 5, 0, 1, 2, 3, 4, 5, 0, 1, 2, 3, 4, 5, 6, 0, 1, 2, 3, 4, 5, 6, 0, 1, 2, 3, 4, 5, 6, 0, 1, 2, 3, 4, 5, 6, 7, 0, 1, 2, 3, 4, 5, 6, 7, 0, 1, 2, 3, 4 (list; graph; listen)
OFFSET

0,12

LINKS

M. Somos, Sequences used for indexing triangular or square arrays

EXAMPLE

0; 0; 0; 0,1; 0,1; 0,1; 0,1,2; 0,1,2; 0,1,2; 0,1,2,3; 0,1,2,3; ...

PROGRAM

(PARI) a(n)=local(m); m=floor(sqrt(6*n+6)-3/2)\3+1; (n-3*binomial(m, 2))%m

CROSSREFS

Sequence in context: A025858 A025684 A025678 this_sequence A025855 A097203 A025850

Adjacent sequences: A073186 A073187 A073188 this_sequence A073190 A073191 A073192

KEYWORD

nonn

AUTHOR

Michael Somos, Jul 19 2002

A084964 Follow n by n-2. Also solution of a(n+2)=a(n)+1, a(0)=2, a(1)=0. +0
16
2, 0, 3, 1, 4, 2, 5, 3, 6, 4, 7, 5, 8, 6, 9, 7, 10, 8, 11, 9, 12, 10, 13, 11, 14, 12, 15, 13, 16, 14, 17, 15, 18, 16, 19, 17, 20, 18, 21, 19, 22, 20, 23, 21, 24, 22, 25, 23, 26, 24, 27, 25, 28, 26, 29, 27, 30, 28, 31, 29, 32, 30, 33, 31, 34, 32, 35, 33, 36, 34, 37, 35, 38, 36, 39 (list; graph; listen)
OFFSET

0,1

LINKS

Index entries for two-way infinite sequences

FORMULA

G.f.: (2-2x+x^2)/((1-x)(1-x^2)). a(2n+1)=n. a(2n)=n+2. a(n+2)=a(n)+1. a(n)=-a(-3-n).

a(n) = floor(n/2) + 1 + (-1)^n. - Reinhard Zumkeller (reinhard.zumkeller(AT)gmail.com), Aug 27 2005

A112032(n)=2^a(n); A112033(n)=3*2^a(n); a(n)=A109613(n+2)-A052938(n). - Reinhard Zumkeller (reinhard.zumkeller(AT)gmail.com), Aug 27 2005

MATHEMATICA

lst={}; a=1; Do[a=n-a; AppendTo[lst, a], {n, 0, 6!}]; lst [From Vladimir Orlovsky (4vladimir(AT)gmail.com), Dec 14 2008]

PROGRAM

(PARI) a(n)=n\2-2*(n%2)+2

CROSSREFS

Cf. A030451.

Cf. A097065.

Cf. A152832 [From Vladimir Orlovsky (4vladimir(AT)gmail.com), Dec 14 2008]

Sequence in context: A025636 A025637 A097065 this_sequence A008720 A008734 A053445

Adjacent sequences: A084961 A084962 A084963 this_sequence A084965 A084966 A084967

KEYWORD

nonn

AUTHOR

Michael Somos, Jun 15 2003

A107920 Lucas and Lehmer numbers with parameters (1+-sqrt(-7))/2. +0
15
0, 1, 1, -1, -3, -1, 5, 7, -3, -17, -11, 23, 45, -1, -91, -89, 93, 271, 85, -457, -627, 287, 1541, 967, -2115, -4049, 181, 8279, 7917, -8641, -24475, -7193, 41757, 56143, -27371, -139657, -84915, 194399, 364229, -24569, -753027, -703889, 802165, 2209943, 605613, -3814273 (list; graph; listen)
OFFSET

0,5

COMMENT

This is an example of a sequence of Lehmer numbers. In this case, the two parameters, alpha and beta, are (1 +- i Sqrt(7))/2. Bilu, Hanrot, Voutier and Mignotte show that all terms of a Lehmer sequence a(n) have a primitive factor for n > 30. Note that for this sequence, a(30) = 24475 = 5*5*11*89 has no primitive factors. - T. D. Noe (noe(AT)sspectra.com), Oct 29 2003

Row sums of Riordan array (1/(1+2x^2),x/(1+2x^2)). - Paul Barry (pbarry(AT)wit.ie), Sep 10 2005

LINKS

Y. Bilu, G. Hanrot, P. M. Voutier and M. Mignotte, Existence of primitive divisors of Lucas and Lehmer numbers

Eric Weisstein's World of Mathematics, Lehmer Number

FORMULA

G.f.; x/(1-x+2x^2). a(n)=a(n-1)-2*a(n-2).

a(n+1)=sum{k=0..n, C((n+k)/2, k)*(-2)^((n-k)/2)*(1+(-1)^(n-k))/2}; a(n+1)=sum{k=0..floor(n/2), C(n-k, k)(-2)^k}; - Paul Barry (pbarry(AT)wit.ie), Sep 10 2005

a(n+1)=Sum_{k, 0<=k<=n} A109466(n,k)*2^(n-k). [From Philippe DELEHAM (kolotoko(AT)wanadoo.fr), Oct 26 2008]

MAPLE

a:= n-> (Matrix([[1, 1], [ -2, 0]])^n)[1, 2]: seq (a(n), n=0..45); [From Alois P. Heinz (heinz(AT)hs-heilbronn.de), Sep 03 2008]

PROGRAM

(PARI) a(n)=if(n<0, 0, imag(quadgen(-7)^n))

(Other) sage: [lucas_number1(n, 1, +2) for n in xrange(0, 46)] # [From Zerinvary Lajos (zerinvarylajos(AT)yahoo.com), Apr 22 2009]

CROSSREFS

A001607(n)=-(-1)^n*a(n).

Sequence in context: A001607 A167433 A077020 this_sequence A159285 A021080 A049764

Adjacent sequences: A107917 A107918 A107919 this_sequence A107921 A107922 A107923

KEYWORD

sign

AUTHOR

Michael Somos, May 28 2005

A078812 Triangle read by rows: T(n,k) = binomial(n+k-1,2*k-1). +0
14
1, 2, 1, 3, 4, 1, 4, 10, 6, 1, 5, 20, 21, 8, 1, 6, 35, 56, 36, 10, 1, 7, 56, 126, 120, 55, 12, 1, 8, 84, 252, 330, 220, 78, 14, 1, 9, 120, 462, 792, 715, 364, 105, 16, 1, 10, 165, 792, 1716, 2002, 1365, 560, 136, 18, 1, 11, 220, 1287, 3432, 5005, 4368, 2380, 816, 171, 20 (list; table; graph; listen)
OFFSET

0,2

COMMENT

Apart from signs, identical to A053122.

Coefficient array for Morgan-Voyce polynomial B(n,x); see A085478 for references. DELEHAM Philippe (kolotoko(AT)wanadoo.fr), Feb 16 2004

T(n,k)=number of compositions of n having k parts when there are q kinds of part q (q=1,2,...). Example: T(4,2)=10 because we have (1,3),(1,3'),(1,3"), (3,1),(3',1),(3",1),(2,2),(2,2'),(2',2) and (2',2'). - Emeric Deutsch (deutsch(AT)duke.poly.edu), Apr 09 2005

T(n, k) is also the number of idempotent order-preserving full transformations (of an n-chain) of height k (height(alpha) = |Im(alpha)|). [From A. Umar (aumarh(AT)squ.edu.om), Oct 02 2008]

REFERENCES

Laradji, A. and Umar, A. Combinatorial results for semigroups of order-preserving full transformations. Semigroup Forum 72 (2006), 51-62. [From A. Umar (aumarh(AT)squ.edu.om), Oct 02 2008]

LINKS

T. D. Noe, Rows n=0..50 of triangle, flattened

FORMULA

G.f.: x*y/(1-(2+y)*x+x^2). To get row n, expand this in powers of x then expand the coefficient of x^n in increasing powers of y.

If indexing begins at 0 we have: T(n, k) = (n+k+1)!/((n-k)!*(2k+1))!. T(n, k) = Sum_{j>=0} T(n-1-j, k-1)*(j+1) with T(n, 0) = n+1, T(n, k) = 0 if n<k. T(n, k) = T(n-1, k-1) + T(n-1, k) + Sum_{j>=0} (-1)^j*T(n-1, k+j)*A000108(j) with T(n, k) = 0 if k<0, T(0, 0)=1 and T(0, k) = 0 for k>0. G.f. for the column k : Sum_{n>=0} T(n, k)*x^n = (x^k)/(1-x)^(2k+2). Row sums : Sum_{k>=0} T(n, k) = A001906(n+1). - DELEHAM Philippe (kolotoko(AT)wanadoo.fr), Feb 16 2004

Diagonal sums are A000079(n)=sum{k=0..floor(n/2), binomial(n+k+1, n-k)}. - Paul Barry (pbarry(AT)wit.ie), Jun 21 2004

Riordan array (1/(1-x)^2, x/(1-x)^2). - Paul Barry (pbarry(AT)wit.ie), Oct 22 2006

EXAMPLE

Triangle begins:

.........................1

........................2,.1

......................3,.4,.1

....................4,.10,.6,.1

..................5,.20,.21,.8,.1

................6,.35,.56,.36,.10,.1

.............7,.56,.126,.120,.55,.12,.1

..........8,.84,.252,.330,.220,.78,.14,.1

MAPLE

for n from 1 to 11 do seq(binomial(n+k-1, 2*k-1), k=1..n) od; # yields sequence in triangular form (Deutsch)

PROGRAM

(PARI) T(n, k)=if(n<0, 0, binomial(n+k-1, 2*k-1))

(PARI) {T(n, k)=polcoeff( polcoeff( x*y/(1-(2+y)*x+x^2) +x*O(x^n), n), k)}

CROSSREFS

This triangle is formed from odd-numbered rows of triangle A011973 read in reverse order.

The column sequences are A000027, A000292, A000389, A000580, A000582, A001288 for k=1..6, resp. For k=7..24 they are A010966..(+2)..A011000 and for k=25..50 they are A017713..(+2)..A017763.

Cf. A053123, A049310. Row sums give A001906.

With signs: A053122.

Cf. A119900 [From Philippe DELEHAM (kolotoko(AT)wanadoo.fr), Dec 02 2008]

Sequence in context: A137614 A143326 A053122 this_sequence A104711 A133112 A159856

Adjacent sequences: A078809 A078810 A078811 this_sequence A078813 A078814 A078815

KEYWORD

easy,nice,nonn,tabl

AUTHOR

Michael Somos, Dec 05, 2002

EXTENSIONS

Edited by N. J. A. Sloane (njas(AT)research.att.com), Apr 28 2008

A080995 Characteristic function of generalized pentagonal numbers. +0
14
1, 1, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0 (list; graph; listen)
OFFSET

0,1

COMMENT

Repeatedly [1,[0,]^2k,1,[0,]^k], k>=0; characteristic function of generalized pentagonal numbers: a(A001318(n))=1, a(A118300(n))=0. - Reinhard Zumkeller (reinhard.zumkeller(AT)gmail.com), Apr 22 2006

REFERENCES

P. A. MacMahon, Combinatory Analysis, Cambridge Univ. Press, London and New York, Vol. 1, 1915 and Vol. 2, 1916; see vol. 2, p 81, Article 331.

LINKS

T. D. Noe, Table of n, a(n) for n=0..1001

Index entries for characteristic functions

Eric Weisstein's World of Mathematics, Jacobi Theta Functions

FORMULA

G.f.: Sum x^(n*(3n+1)/2), n=-inf..inf [the exponents are the pentagonal numbers, A000326].

a(n)=b(24n+1) where b(n) is multiplicative and b(2^e)=b(3^e)=0^e, b(p^e)=(1+(-1)^e)/2 if p>3. - Michael Somos Jun 06 2005

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

Expansion of phi(-q^3) / chi(-q) in powers of q where phi(), chi() are Ramanujan theta functions. - Michael Somos Sep 14 2007

Expansion of psi(q) - q * psi(q^9) in powers of q^3 where psi() is a Ramanujan theta function. - Michael Somos Sep 14 2007

Expansion of f(x, x^2) in powers of x where f() is Ramanujan's two-variable theta function.

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

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

G.f.: Product_{k>0} (1 - x^(3*k)) / (1 - x^k + x^(2*k)). - Michael Somos Jan 26 2008

EXAMPLE

q + q^25 + q^49 + q^121 + q^169 + q^289 + q^361 + q^529 + q^625 + ...

PROGRAM

(PARI) a(n)=if(n<0, 0, abs(polcoeff(eta(x+x*O(x^n)), n)))

(PARI) a(n)=issquare(24*n+1) /* Michael Somos Apr 13 2005 */

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

CROSSREFS

|A010815(n)| = a(n). A089806(2n) = a(n). A033683(24n+1) = a(n).

Sequence in context: A115513 A133080 A010815 this_sequence A121373 A133985 A143062

Adjacent sequences: A080992 A080993 A080994 this_sequence A080996 A080997 A080998

KEYWORD

nonn,easy

AUTHOR

Michael Somos, Feb 27, 2003

A092848 Expansion of reciprocal of Hauptmodul for Gamma_0(18). +0
14
1, -1, 0, 2, -2, -1, 4, -4, -1, 8, -8, -2, 14, -14, -4, 24, -23, -6, 40, -38, -10, 63, -60, -16, 98, -92, -24, 150, -140, -36, 224, -208, -54, 329, -304, -78, 478, -440, -112, 684, -627, -160, 968, -884, -224, 1358, -1236, -312, 1884, -1710, -432, 2592, -2346, -590, 3540, -3196, -801, 4796, -4320, -1082, 6454 (list; graph; listen)
OFFSET

0,4

REFERENCES

B. C. Berndt, Ramanujan's Notebooks Part III, Springer-Verlag, see p. 345 Entry 1(i).

W. Duke, Continued fractions and modular functions, Bull. Amer. Math. Soc. 42 (2005), 137-162. See page 155 Eq. (9.13)

FORMULA

G.f.: Product_{k>0} (1-x^(2k-1))/(1-x^(6k-3))^3.

Given g.f. A(x), then B(x)=x*A(x^3) satisfies 0=f(B(x), B(x^2)) where f(u, v)=u^2-v+2uv^2.

Given g.f. A(x), then B(x)=x*A(x^3) satisfies 0=f(B(x), B(x^3)) where f(u, v)=(v^3-v^2+v)-u^3(1+2v+4v^2).

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

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

Expansion of chi(-q)/chi(-q^3)^3 where chi() (g.f. A000700) is a Ramanujan theta function.

G.f.: 1/(1+ (x+x^2)/(1+ (x^2+x^4)/(1 +(x^3+x^6)/...))).

Expansion of q^(-1/3)c(q^2)/c(q) where c() is a cubic AGM analog function. - Michael Somos Oct 04 2006

EXAMPLE

q - q^4 + 2*q^10 - 2*q^13 - q^16 + 4*q^19 - 4*q^22 - q^25 + 8*q^28 + ...

PROGRAM

(PARI) {a(n)=local(A, m); if(n<0, 0, A=1+O(x); m=1; while(m<=n, m*=2; A=subst(A, x, x^2); A=sqrt(A+(x*A^2)^2)-x*A^2); polcoeff(A, n))}

(PARI) {a(n)=if(n<0, 0, polcoeff(prod(k=0, (n-1)\2, (1-x^(2*k+1))^if(k%3==1, -2, 1), 1+x*O(x^n)), n))}

CROSSREFS

A062242(2*n + 1) = a(n). A128111(n) = (-1)^n * a(n). Convolution inverse of A062242.

Sequence in context: A129862 A138189 A110090 this_sequence A128111 A107356 A124725

Adjacent sequences: A092845 A092846 A092847 this_sequence A092849 A092850 A092851

KEYWORD

sign

AUTHOR

Michael Somos, Mar 07 2004

A079006 Expansion of q^(-1/4) * (eta(q) * eta(q^4)^2 / eta(q^2)^3)^2 in powers of q. +0
13
1, -2, 5, -10, 18, -32, 55, -90, 144, -226, 346, -522, 777, -1138, 1648, -2362, 3348, -4704, 6554, -9056, 12425, -16932, 22922, -30848, 41282, -54946, 72768, -95914, 125842, -164402, 213901, -277204, 357904, -460448, 590330, -754368, 960948, -1220370 (list; graph; listen)
OFFSET

0,2

REFERENCES

A. Cayley, A memoir on the transformation of elliptic functions, Collected Mathematical Papers. Vols. 1-13, Cambridge Univ. Press, London, 1889-1897, Vol. 9, p. 128.

H. R. P. Ferguson, D. E. Nielsen and G. Cook, A partition formula for the integer coefficients of the theta function nome, Math. Comp., 29 (1975), 851-855.

N. J. Fine, Basic Hypergeometric Series and Applications, Amer. Math. Soc., 1988; Eq. (34.3).

FORMULA

a(n) = (2/n)*Sum_{k=1..n} (-1)^k*A046897(k)*a(n-k). - Vladeta Jovovic (vladeta(AT)eunet.rs), Dec 24 2002

Expansion of q^(-1/4)(1/2)k^(1/2) in powers of q.

Expansion of (1/q)(1/2)(1-sqrt(k'))/(1+sqrt(k')) in powers of q^4.

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

G.f. A(x) satisfies A(x)^2=A(x^2)/(1+4*x*A(x^2)^2). - Michael Somos Mar 19 2004

Given g.f. A(x), then B(x)=x*A(x^4) satisfies 0=f(B(x), B(x^2)) where f(u, v)=u^2(1+4v^2)-v . - Michael Somos Jul 09 2005

Given g.f. A(x), then B(x)=x*A(x^4) satisfies 0=f(B(x), B(x^2), B(x^3), B(x^6)) where f(u1, u2, u3, u6)=u1*u3*(u6+u2)^2-u2*u6 . - Michael Somos Jul 09 2005

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

Expansion of continued fraction 1/(1-x^2+(x^1+x^3)^2/(1-x^6+(x^2+x^6)^2/(1-x^10+(x^3+x^9)^2/...))). - Michael Somos Sep 01 2005

Given g.f. A(x), then B(x)=2*x*A(x^4) satisfies 0=f(B(x), B(x^3)) where f(u, v)=(1-u^4)(1-v^4)-(1-uv)^4 . - Michael Somos Jan 01 2006

Expansion of f(-q^4)^2 / f(q)^2 = psi(q)^2 / phi(q)^2 = psi(q^2)^2 / psi(q)^2 = psi(-q)^2 / phi(-q^2)^2 = 1 / (chi(q)^2 * chi(-q^2)^2) = 1 / (chi(q)^4 * chi(-q)^2) = chi(-q)^2 / chi(-q^2)^4 = psi(q^2) / phi(q) in powers of q where phi(), psi(), chi(), f() a re Ramanujan theta functions.

Expansion of psi(q^2) / phi(q) = psi(q)^2 / phi(q)^2 = psi(q^2)^2 / psi(q)^2 = psi(-q)^2 / phi(-q^2)^2 = chi(-q)^2 / chi(-q^2)^4 = 1 / (chi(q)^2 * chi(-q^2)^2) = 1 / (chi(q)^4 * chi(-q)^2) = f(-q^4)^2 / f(q)^2 in powers of q where phi(), psi(), chi(), f() are Ramanujan theta functions.

EXAMPLE

q - 2*q^5 + 5*q^9 - 10*q^13 + 18*q^17 - 32*q^21 + 55*q^25 - 90*q^29 + ...

PROGRAM

(PARI) {a(n)=local(N, A); if(n<0, 0, N=(sqrtint(16*n+1)+1)\2; A=contfracpnqn( matrix(2, N, i, j, if(i==1, if(j<2, 1+O(x^(N^2+N)), (x^(j-1)+x^(3*j-3))^2), 1-x^(4*j-2)))); polcoeff(A[2, 1]/A[1, 1], 4*n))} /* Michael Somos Sep 01 2005 */

(PARI) {a(n) = local(A, m); if( n<0, 0, A = 1 + O(x); m = 1; while( m<=n, m*=2; A = subst(A, x, x^2); A = sqrt(A / (1 + 4 * x*A^2))); polcoeff(A, n))}

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

CROSSREFS

Convolution inverse of A029839. Convolution square of A083365.

See A127391, A127392, A001936 for other versions of this sequence.

a(n)=(-1)^n A001936(n).

Cf. A002103, A046897.

Sequence in context: A034350 A006327 A103577 this_sequence A001936 A127297 A018739

Adjacent sequences: A079003 A079004 A079005 this_sequence A079007 A079008 A079009

KEYWORD

sign,easy

AUTHOR

Michael Somos, Dec 22 2002

page 1 2 3 4 5 6 7 8 9 10 11 ... 99

Search completed in 0.037 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 20 21:36 EST 2009. Contains 167244 sequences.


AT&T Labs Research