cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-7 of 7 results.

A069910 Expansion of Product_{i in A069908} 1/(1 - x^i).

Original entry on oeis.org

1, 0, 1, 1, 2, 2, 3, 3, 5, 5, 7, 8, 11, 12, 16, 18, 23, 26, 33, 37, 46, 52, 63, 72, 87, 98, 117, 133, 157, 178, 209, 236, 276, 312, 361, 408, 471, 530, 609, 686, 784, 881, 1004, 1126, 1279, 1433, 1621, 1814, 2048, 2286, 2574, 2871, 3223, 3590, 4022, 4472, 5000
Offset: 0

Views

Author

N. J. A. Sloane, May 05 2002

Keywords

Comments

Number 39 of the 130 identities listed in Slater 1952.
Number of partitions of 2*n into distinct odd parts. - Vladeta Jovovic, May 08 2003

Examples

			G.f. = 1 + x^2 + x^3 + 2*x^4 + 2*x^5 + 3*x^6 + 3*x^7 + 5*x^8 + 5*x^9 + ...
G.f. = q^-1 + q^95 + q^143 + 2*q^191 + 2*q^239 + 3*q^287 + 3*q^335 + ...
		

References

  • M. D. Hirschhorn, The Power of q, Springer, 2017. Chapter 19, Exercises p. 173.

Crossrefs

Programs

  • Maple
    a:= proc(n) option remember; `if`(n=0, 1,
          add(add(d*[0$2, 1$4, 0$5, 1$4, 0][irem(d, 16)+1],
          d=numtheory[divisors](j))*a(n-j), j=1..n)/n)
        end:
    seq(a(n), n=0..80);  # Alois P. Heinz, Apr 01 2014
  • Mathematica
    max = 56; p = Product[1/(1-x^i), {i, Select[Range[max], MemberQ[{2, 3, 4, 5, 11, 12, 13, 14}, Mod[#, 16]]&]}]; s = Series[p, {x, 0, max}]; a[n_] := Coefficient[s, x, n]; Table[a[n], {n, 0, max}] (* Jean-François Alcover, Apr 09 2014 *)
    nmax=60; CoefficientList[Series[Product[(1-x^(8*k-1))*(1-x^(8*k-7))*(1-x^(8*k))*(1-x^(16*k-6))*(1-x^(16*k-10))/(1-x^k), {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Oct 04 2015 *)
    a[ n_] := SeriesCoefficient[ Product[ (1 - x^k)^-{ 0, 0, 1, 1, 1, 1, 0, 0, 0, 0, 0, 1, 1, 1, 1, 0 }[[ Mod[k, 16] + 1]], {k, n}], {x, 0, n}]; (* Michael Somos, Apr 14 2016 *)
  • PARI
    {a(n) = my(A); if( n<0,0, n=2*n; A = x * O(x^n); polcoeff( eta(-x + A) / eta(x^2 + A), n))}; /* Michael Somos, Apr 11 2004 */
    
  • PARI
    N=66;  q='q+O('q^N);  S=1+sqrtint(N);
    gf=sum(n=0, S, q^(2*n^2) / prod(k=1, 2*n, 1-q^k ) );
    Vec(gf)  \\ Joerg Arndt, Apr 01 2014
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( prod(k=1, n, (1 - x^k + x * O(x^n))^-[ 0, 0, 1, 1, 1, 1, 0, 0, 0, 0, 0, 1, 1, 1, 1, 0][k%16 + 1]), n))}; /* Michael Somos, Apr 14 2016 */

Formula

Euler transform of period 16 sequence [0, 1, 1, 1, 1, 0, 0, 0, 0, 0, 1, 1, 1, 1, 0, 0, ...]. - Michael Somos, Apr 11 2004
G.f.: Sum_{n>=0} q^(2*n^2) / Product_{k=1..2*n} (1 - q^k). - Joerg Arndt, Apr 01 2014
a(n) ~ exp(sqrt(n/3)*Pi) / (2^(5/2) * 3^(1/4) * n^(3/4)). - Vaclav Kotesovec, Oct 04 2015
Expansion of f(x^3, x^5) / f(-x^2) in powers of x where f(, ) is Ramanujan's general theta function. - Michael Somos, Apr 14 2016
a(n) = A000700(2*n).
a(n) = A027356(4n+1,2n+1). - Alois P. Heinz, Oct 28 2019
From Peter Bala, Feb 08 2021: (Start)
G.f.: A(x) = Product_{n >= 1} (1 + x^(4*n))^2*(1 + x^(4*n-2))*(1 + x^(8*n-3))*(1 + x^(8*n-5)).
The 2 X 2 matrix Product_{k >= 0} [1, x^(2*k+1); x^(2*k+1), 1] = [A(x^2), x*B(x^2); x*B(x)^2, A(x^2)], where B(x) is the g.f. of A069911.
A(x^2) + x*B(x^2) = A^2(-x) + x*B^2(-x) = Product_{k >= 0} 1 + x^(2*k+1), the g.f. of A000700.
A^2(x) + x*B^2(x) is the g.f. of A226622.
(A^2(x) + x*B^2(x))/(A^2(x) - x*B^2(x)) is the g.f. of A208850.
A^4(sqrt(x)) - x*B^4(sqrt(x)) is the g.f. of A029552.
A(x)*B(x) is the g.f. of A226635; A(-x)/B(-x) is the g.f. of A111374; B(-x)/A(-x) is the g.f. of A092869. (End)

A083906 Table read by rows: T(n, k) is the number of length n binary words with exactly k inversions.

Original entry on oeis.org

1, 2, 3, 1, 4, 2, 2, 5, 3, 4, 3, 1, 6, 4, 6, 6, 6, 2, 2, 7, 5, 8, 9, 11, 9, 7, 4, 3, 1, 8, 6, 10, 12, 16, 16, 18, 12, 12, 8, 6, 2, 2, 9, 7, 12, 15, 21, 23, 29, 27, 26, 23, 21, 15, 13, 7, 4, 3, 1, 10, 8, 14, 18, 26, 30, 40, 42, 48, 44, 46, 40, 40, 30, 26, 18, 14, 8, 6, 2, 2
Offset: 0

Views

Author

Alford Arnold, Jun 19 2003

Keywords

Comments

There are A033638(n) values in the n-th row, compliant with the order of the polynomial.
In the example for n=6 detailed below, the orders of [6, k]_q are 1, 6, 9, 10, 9, 6, 1 for k = 0..6,
the maximum order 10 defining the row length.
Note that 1 6 9 10 9 6 1 and related distributions are antidiagonals of A077028.
A083480 is a variation illustrating a relationship with numeric partitions, A000041.
The rows are formed by the nonzero entries of the columns of A049597.
If n is even the n-th row converges to n+1, n-1, n-4, ..., 19, 13, 7, 4, 3, 1 which is A029552 reversed, and if n is odd the sequence is twice A098613. - Michael Somos, Jun 25 2017

Examples

			When viewed as an array with A033638(r) entries per row, the table begins:
. 1 ............... : 1
. 2 ............... : 2
. 3 1 ............. : 3 + q = (1) + (1+q) + (1)
. 4 2 2 ........... : 4 + 2q + 2q^2 = 1 + (1+q+q^2) + (1+q+q^2) + 1
. 5 3 4 3 1 ....... : 5 + 3q + 4q^2 + 3q^3 + q^4
. 6 4 6 6 6 2 2
. 7 5 8 9 11 9 7 4 3 1
. 8 6 10 12 16 16 18 12 12 8 6 2 2
. 9 7 12 15 21 23 29 27 26 23 21 15 13 7 4 3 1
...
The second but last row is from the sum over 7 q-polynomials coefficients:
. 1 ....... : 1 = [6,0]_q
. 1 1 1 1 1 1 ....... : 1+q+q^2+q^3+q^4+q^5 = [6,1]_q
. 1 1 2 2 3 2 2 1 1 ....... : 1+q+2q^2+2q^3+3q^4+2q^5+2q^6+q^7+q^8 = [6,2]_q
. 1 1 2 3 3 3 3 2 1 1 ....... : 1+q+2q^2+3q^3+3q^4+3q^5+3q^6+2q^7+q^8+q^9 = [6,3]_q
. 1 1 2 2 3 2 2 1 1 ....... : 1+q+2q^2+2q^3+3q^4+2q^5+2q^6+q^7+q^8 = [6,4]_q
. 1 1 1 1 1 1 ....... : 1+q+q^2+q^3+q^4+q^5 = [6,5]_q
. 1 ....... : 1 = [6,6]_q
		

References

  • George E. Andrews, 'Theory of Partitions', 1976, page 242.

Crossrefs

Programs

  • Magma
    R:=PowerSeriesRing(Rationals(), 100);
    qBinom:= func< n,k,x | n eq 0 or k eq 0 select 1 else (&*[(1-x^(n-j))/(1-x^(j+1)): j in [0..k-1]]) >;
    A083906:= func< n,k | Coefficient(R!((&+[qBinom(n,k,x): k in [0..n]]) ), k) >;
    [A083906(n,k): k in [0..Floor(n^2/4)], n in [0..12]]; // G. C. Greubel, Feb 13 2024
    
  • Maple
    QBinomial := proc(n,m,q) local i ; factor( mul((1-q^(n-i))/(1-q^(i+1)),i=0..m-1) ) ; expand(%) ; end:
    A083906 := proc(n,k) add( QBinomial(n,m,q),m=0..n ) ; coeftayl(%,q=0,k) ; end:
    for n from 0 to 10 do for k from 0 to A033638(n)-1 do printf("%d,",A083906(n,k)) ; od: od: # R. J. Mathar, May 28 2009
    T := proc(n, k) if n < 0 or k < 0 or k > floor(n^2/4) then return 0 fi;
    if n < 2 then return n + 1 fi; 2*T(n-1, k) - T(n-2, k) + T(n-2, k - n + 1) end:
    seq(print(seq(T(n, k), k = 0..floor((n/2)^2))), n = 0..8);  # Peter Luschny, Feb 16 2024
  • Mathematica
    Table[CoefficientList[Total[Table[FunctionExpand[QBinomial[n, k, q]], {k, 0, n}]],q], {n, 0, 10}] // Grid (* Geoffrey Critzer, May 14 2017 *)
  • PARI
    {T(n, k) = polcoeff(sum(m=0, n, prod(k=0, m-1, (x^n - x^k) / (x^m - x^k))), k)}; /* Michael Somos, Jun 25 2017 */
    
  • SageMath
    def T(n,k): # T = A083906
        if k<0 or k> (n^2//4): return 0
        elif n<2 : return n+1
        else: return 2*T(n-1, k) - T(n-2, k) + T(n-2, k-n+1)
    flatten([[T(n,k) for k in range(int(n^2//4)+1)] for n in range(13)]) # G. C. Greubel, Feb 13 2024

Formula

T(n, k) is the coefficient [q^k] of the Sum_{m=0..n} [n, m]_q over q-Binomial coefficients.
Row sums: Sum_{k=0..floor(n^2/4)} T(n,k) = 2^n.
For n >= k, T(n+1,k) = T(n, k) + A000041(k). - Geoffrey Critzer, Feb 12 2021
Sum_{k=0..floor(n^2/4)} (-1)^k*T(n, k) = A060546(n). - G. C. Greubel, Feb 13 2024
From Mikhail Kurkov, Feb 14 2024: (Start)
T(n, k) = 2*T(n-1, k) - T(n-2, k) + T(n-2, k - n + 1) for n >= 2 and 0 <= k <= floor(n^2/4).
Sum_{i=0..n} T(n-i, i) = A000041(n+1). Note that upper limit of the summation can be reduced to A083479(n) = (n+2) - ceiling(sqrt(4*n)).
Both results were proved (see MathOverflow link for details). (End)
From G. C. Greubel, Feb 17 2024: (Start)
T(n, floor(n^2/4)) = A000034(n).
Sum_{k=0..floor(n^2/4)} (-1)^k*T(n, k) = A016116(n+1).
Sum_{k=0..(n + 2) - ceiling(sqrt(4*n))} (-1)^k*T(n - k, k) = (-1)^n*A000025(n+1) = -A260460(n+1). (End)

Extensions

Edited by R. J. Mathar, May 28 2009
New name using a comment from Geoffrey Critzer by Peter Luschny, Feb 17 2024

A098613 Expansion of psi(x^2) / f(-x) in powers of x where psi(), f() are Ramanujan theta functions.

Original entry on oeis.org

1, 1, 3, 4, 7, 10, 17, 23, 35, 48, 69, 93, 131, 173, 236, 310, 413, 536, 704, 903, 1170, 1489, 1904, 2403, 3044, 3811, 4784, 5951, 7409, 9157, 11325, 13912, 17095, 20891, 25519, 31029, 37708, 45632, 55184, 66495, 80050, 96064, 115173, 137680, 164425, 195860
Offset: 0

Views

Author

Michael Somos, Sep 17 2004

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
This sequence convolved with A000009 gives A001936. - Gary W. Adamson, Mar 24 2011
a(n) is the number of partitions of n in which each odd part can occur any number of times but each even part is of two kinds and each kind can occur at most once. - Michael Somos, Dec 01 2019

Examples

			G.f. = 1 + x + 3*x^2 + 4*x^3 + 7*x^4 + 10*x^5 + 17*x^6 + 23*x^7 + ...
G.f. = q^5 + q^29 + 3*q^53 + 4*q^77 + 7*q^101 + 10*q^125 + 17*q^149 + ...
		

Crossrefs

Programs

  • Mathematica
    a[ n_] := SeriesCoefficient[ EllipticTheta[ 2, 0, x] / (2 x^(1/4) QPochhammer[ x]), {x, 0, n}]; (* Michael Somos, Oct 29 2013 *)
    a[ n_] := SeriesCoefficient[ 1 / (QPochhammer[ x, x^2] QPochhammer[ x^2, x^4]^2), {x, 0, n}]; (* Michael Somos, Oct 29 2013 *)
    nmax = 40; CoefficientList[ Series[Product[(1 + x^k) * (1 + x^(2*k))^2, {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Sep 07 2015 *)
    a[ n_] := SeriesCoefficient[ QPochhammer[ -x, x^2] QPochhammer[ -x^2, x^2]^3, {x, 0, n}]; (* Michael Somos, Sep 07 2015 *)
    a[ n_] := SeriesCoefficient[ QPochhammer[ x, -x]^2 QPochhammer[ -x, x]^3, {x, 0, n}]; (* Michael Somos, Sep 07 2015 *)
    a[ n_] := SeriesCoefficient[ QPochhammer[ -x, x] QPochhammer[ -x^2, x^2]^2, {x, 0, n}]; (* Michael Somos, Sep 07 2015 *)
  • PARI
    {a(n) = if( n<0, 0, polcoeff( sum(k=0, (sqrtint(4*n+1)-1)\2, x^(k^2+k)) / eta(x + x * O(x^n)), n))};
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^4 + A)^2 / (eta(x + A) * eta(x^2 + A)), n))};
    
  • PARI
    {a(n) = if( n<0,0, polcoeff( sum(k=0, 2*n+1, prod(i=1, k, (1 - x^(2*n+2 - i)) / (1 - x^i))) / 2, n^2))};

Formula

Expansion of chi(x) / chi(-x^2)^3 = 1 / (chi(-x)* chi(-x^2)^2) = 1 / (chi(x)^2 * chi(-x)^3) in powers of x where chi() is a Ramanujan theta function. - Michael Somos, Sep 07 2015
Expansion of q^(-5/24) * eta(q^4)^2 / (eta(q) * eta(q^2)) in powers of q.
Euler transform of period 4 sequence [1, 2, 1, 0, ...].
G.f. A(x) is the limit of x^(n^2+n) * P_{2*n+1}(1/x)/2 where P_n(q) = Sum_{k=0..n} C(n, k; q) and C(n, k; q) is the q-binomial coefficients. See A083906 for P_n(q).
G.f.: (Sum_{k>0} x^(k^2-k)) / (Product_{k>0} (1 - x^k)).
a(n) ~ exp(Pi*sqrt(2*n/3)) / (2^(11/4) * 3^(1/4) * n^(3/4)). - Vaclav Kotesovec, Sep 07 2015
G.f. is a period 1 Fourier series which satisfies f(-1 / (576 t)) = 8^(-1/2) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A143161. - Michael Somos, Sep 07 2015
G.f.: Product_{k>=1} (1 + x^(2*k))^2 / (1 - x^(2*k-1)). - Michael Somos, Dec 01 2019

A143161 Expansion of chi(-x)^2 * chi(-x^2) in powers of x where chi() is a Ramanujan theta function.

Original entry on oeis.org

1, -2, 0, 0, 3, -2, 0, 0, 4, -6, 0, 0, 7, -8, 0, 0, 13, -14, 0, 0, 19, -20, 0, 0, 29, -34, 0, 0, 43, -46, 0, 0, 62, -70, 0, 0, 90, -96, 0, 0, 126, -138, 0, 0, 174, -186, 0, 0, 239, -262, 0, 0, 325, -346, 0, 0, 435, -472, 0, 0, 580, -620, 0, 0, 769, -826, 0, 0
Offset: 0

Views

Author

Michael Somos, Jul 27 2008

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).

Examples

			G.f. = 1 - 2*x + 3*x^4 - 2*x^5 + 4*x^8 - 6*x^9 + 7*x^12 - 8*x^13 + 13*x^16 + ...
G.f. = 1/q - 2*q^5 + 3*q^23 - 2*q^29 + 4*q^47 - 6*q^53 + 7*q^71 - 8*q^77 + ...
		

Crossrefs

Programs

  • Mathematica
    a[ n_] := SeriesCoefficient[ QPochhammer[ x, x^2]^2 QPochhammer[ x^2, x^4], {x, 0, n}]; (* Michael Somos, Sep 07 2015 *)
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x + A)^2 / (eta(x^2 + A) * eta(x^4 + A)), n))};

Formula

Expansion of q^(1/6) * eta(q)^2 / (eta(q^2) * eta(q^4)) in powers of q.
Euler transform of period 4 sequence [ -2, -1, -2, 0, ...].
G.f. is a period 1 Fourier series which satisfies f(-1 / (576 t)) = 8^(1/2) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A098613. - Michael Somos, Sep 07 2015
a(4*n + 2) = a(4*n + 3) = 0.
G.f.: (Product_{k>0} (1 + x^k)^2 * (1 + x^(2*k)))^-1.
a(4*n) = A029552(n). a(4*n + 1) = -2 * A098613(n).

A225853 Expansion of phi(x) / f(-x^4) in powers of x where phi(), f() are Ramanujan theta functions.

Original entry on oeis.org

1, 2, 0, 0, 3, 2, 0, 0, 4, 6, 0, 0, 7, 8, 0, 0, 13, 14, 0, 0, 19, 20, 0, 0, 29, 34, 0, 0, 43, 46, 0, 0, 62, 70, 0, 0, 90, 96, 0, 0, 126, 138, 0, 0, 174, 186, 0, 0, 239, 262, 0, 0, 325, 346, 0, 0, 435, 472, 0, 0, 580, 620, 0, 0, 769, 826, 0, 0, 1007, 1072, 0
Offset: 0

Views

Author

Michael Somos, May 17 2013

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).

Examples

			1 + 2*x + 3*x^4 + 2*x^5 + 4*x^8 + 6*x^9 + 7*x^12 + 8*x^13 + 13*x^16 + ...
1/q + 2*q^5 + 3*q^23 + 2*q^29 + 4*q^47 + 6*q^53 + 7*q^71 + 8*q^77 + 13*q^95 + ...
		

Crossrefs

Programs

  • Mathematica
    a[n_]:= SeriesCoefficient[EllipticTheta[3,0,q]/QPochhammer[q^4],{q,0,n}];
    a[n_]:= SeriesCoefficient[QPochhammer[q^2,q^4]^3/QPochhammer[q,q^2]^2, {q,0,n}];
  • PARI
    {a(n) = local(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A)^5 / (eta(x + A)^2 * eta(x^4 + A)^3), n))}

Formula

Expansion of chi(x)^2 * chi(-x^2) = chi(x)^3 * chi(-x) = chi(-x^2)^3 / chi(-x)^2 in powers of x where chi() is a Ramanujan theta function.
Expansion of q^(1/4) * eta(q^2)^5 / (eta(q)^2 * eta(q^4)^3) in powers of q.
Euler transform of period 4 sequence [ 2, -3, 2, 0, ...].
G.f. is a period 1 Fourier series which satisfies f(-1 / (576 t)) = 2^(1/2) g(t) where q = exp(2 Pi i t) and g() is the g.f. of A029552.
G.f.: Product_{k>0} (1 - x^(4*k-2))^3 / (1 - x^(2*k-1))^2 = (Sum_{k in Z} x^k^2) / (Product_{k>0} (1 - x^(4*k))).
a(n) = (-1)^n * A143161(n). a(4*n + 2) = a(4*n + 3) = 0.

A385672 Irregular triangle read by rows: T(n, k) is the number of n-step walks on the square lattice having algebraic area k; n >= 0, 0 <= k <= floor(n^2/4).

Original entry on oeis.org

1, 4, 12, 2, 40, 8, 4, 124, 42, 16, 6, 2, 416, 160, 92, 28, 16, 4, 4, 1348, 678, 362, 174, 88, 34, 22, 8, 6, 2, 4624, 2548, 1624, 756, 460, 200, 156, 56, 40, 20, 12, 4, 4, 15632, 10062, 6336, 3586, 2110, 1106, 742, 388, 278, 152, 82, 46, 34, 14, 8, 6, 2
Offset: 0

Views

Author

Andrei Zabolotskii, Aug 04 2025

Keywords

Comments

Rows can be extended to negative k with T(n, -k) = T(n, k). Sums of such extended rows give 4^n.
The algebraic area is Integral y dx over the walk, which equals (Sum_{steps right} y) - (Sum_{steps left} y).

Examples

			The triangle begins:
     1
     4
    12,   2
    40,   8,   4
   124,  42,  16,   6,  2
   416, 160,  92,  28, 16,  4,  4
  1348, 678, 362, 174, 88, 34, 22, 8, 6, 2
   ...
T(3, 1) = 8: RUR (right, up, right), LUR, RDL, LDL, URU, URD, DLU, DLD.
		

Crossrefs

Row lengths are A033638 = A002620 + 1.
A352838 is an analog that gives the number of closed walks.

Programs

  • Python
    d = [{((0, 0), 0): 1}]
    for _ in range(10):
        nd = {}
        for key, nw in d[-1].items():
            pos, ar = key
            x, y = pos
            for key in [
                ((x+1, y), ar + y),
                ((x-1, y), ar - y),
                ((x, y+1), ar),
                ((x, y-1), ar)
                ]:
                if key in nd:
                    nd[key] += nw
                else:
                    nd[key] = nw
        d.append(nd)
    t = []
    for nd in d:
        a = [0] * (max(ar for _, ar in nd) + 1)
        for key, nw in nd.items():
            _, ar = key
            if ar >= 0:
                a[ar] += nw
        t.append(a)
    print(t)

Formula

It appears that T(2*n, n^2 - k) = 2 * A029552(k) for k < n and T(2*n+1, n^2+n - k) = 4 * A098613(k) for k < n.

A059982 Symmetric array of numeric partitions related to 1 4 9 16 ... and 1 3 4 7 13 ..., read by rows.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 2, 3, 2, 1, 3, 5, 3, 1, 1, 5, 7, 5, 1, 2, 7, 11, 7, 2, 3, 11, 15, 11, 3, 5, 15, 22, 15, 5, 1, 7, 22, 30, 22, 7, 1, 1, 11, 30, 42, 30, 11, 1, 2, 15, 42, 56, 42, 15, 2, 3, 22, 56, 77, 56, 22, 3, 5, 30, 77, 101, 77, 30, 5, 7, 42, 101, 135, 101, 42, 7, 11, 56, 135, 176
Offset: 0

Views

Author

Alford Arnold, Mar 06 2001

Keywords

Comments

Begin with row zero and generate a column of values using the sequence of numeric partitions (A000041). At rows 1 4 9 16 25 ... A000290, generate two new columns one each to the left and right of existing columns. Note that the row sums appear to be A029552.

Examples

			The array begins:
........................1
................1.......1.......1
................1.......2.......1
................2.......3.......2
........1.......3.......5.......3.......1
........1.......5.......7.......5.......1
........2.......7.......11......7.......2
........3.......11......15......11......3
........5.......15......22......15......5
1.......7.......22......30......22......7.......1
1.......11......30......42......30......11......1
2.......15......42......56......42......15......2
3.......22......56......77......56......22......3
5.......30......77......101.....77......30......5
		

References

  • Kass, Moody, Patera and Slansky (1990), Affine Lie Algebras, Weight Multiplicities and Branching Rules. University of California Press. Vol. I, page 108.

Crossrefs

Showing 1-7 of 7 results.