A008731 -id:A008731 - cp's OEIS frontend

A008763 Expansion of g.f.: x^4/((1-x)(1-x^2)^2(1-x^3)).

0, 0, 0, 0, 1, 1, 3, 4, 7, 9, 14, 17, 24, 29, 38, 45, 57, 66, 81, 93, 111, 126, 148, 166, 192, 214, 244, 270, 305, 335, 375, 410, 455, 495, 546, 591, 648, 699, 762, 819, 889, 952, 1029, 1099, 1183, 1260, 1352, 1436, 1536, 1628, 1736, 1836, 1953, 2061, 2187, 2304, 2439

Offset: 0

N. J. A. Sloane, Simon Plouffe, Stephen P. Humphries

Number of 2 X 2 square partitions of n.
1/((1-x^2)*(1-x^4)^2*(1-x^6)) is the Molien series for 4-dimensional representation of a certain group of order 192 [Nebe, Rains, Sloane, Chap. 7].
Number of ways of writing n as n = p+q+r+s so that p >= q, p >= r, q >= s, r >= s with p, q, r, s >= 1. That is, we can partition n as
pq
rs
with p >= q, p >= r, q >= s, r >= s.
The coefficient of s(2n) in s(n,n) * s(n,n) * s(n,n) * s(n,n) is a(n+4), where s(n) is the Schur function corresponding to the trivial representation, s(n,n) is a Schur function corresponding to the two row partition and * represents the inner or Kronecker product of symmetric functions. - Mike Zabrocki, Dec 22 2005
Let F() be the Fibonacci sequence A000045. Let f([x, y, z, w]) = F(x) * F(y) * F(z) * F(w). Let N([x, y, z, w]) = x^2 + y^2 + z^2 + w^2. Let Q(k) = set of all ordered quadruples of integers [x, y, z, w] such that 1 <= x <= y <= z <= w and N([x, y, z, w]) = k. Let P(n) = set of all unordered triples {q1, q2, q3} of elements of some Q(k) such that max(w1, w2, w3) = n and f(q1) + f(q2) = f(q3). Then a(n-1) is the number of elements of P(n). - Michael Somos, Jan 21 2015
Number of partitions of 2n+2 into 4 parts with alternating parity from smallest to largest (or vice versa). - Wesley Ivan Hurt, Jan 19 2021

			a(7) = 4:
41 32 31 22
11 11 21 21
G.f. = x^4 + x^5 + 3*x^6 + 4*x^7 + 7*x^8 + 9*x^9 + 14*x^10 + 17*x^11 + ...
a(5-1) = 1 because P(5) has only one triple {[1,1,1,5], [2,2,2,4], [1,3,3,3]} of elements from Q(28) where f([1,1,1,5]) = 5, f([2,2,2,4]) = 3, f([1,3,3,3]) = 8, and 5 + 3 = 8. - _Michael Somos_, Jan 21 2015
a(6-1) = 1 because P(6) has only one triple {[1,1,2,6], [2,2,3,5], [1,3,4,4]} of elements from Q(42) where f([1,1,2,6]) = 8, f([2,2,3,5]) = 10, f([1,3,4,4]) = 18 and 8 + 10 = 18. - _Michael Somos_, Jan 21 2015
a(7-1) = 3 because P(7) has three triples. The triple {[1,1,1,7], [2,4,4,4], [3,3,3,5]} from Q(52) where f([1,1,1,7]) = 13, f([2,4,4,4]) = 27, f([3,3,3,5]) = 40 and 13 + 27 = 40. The triple {[1,2,2,7], [2,3,3,6], [1,4,4,5]} from Q(58) where f([1,2,2,7]) = 13, f([2,3,3,6]) = 32, f([1,4,4,5]) = 45 and 13 + 32 = 45. The triple {[1,1,3,7], [2,2,4,6], [1,3,5,5]} from Q(60) where f([1,1,3,7]) = 26, f([2,2,4,6]) = 24, f([1,3,5,5]) = 50 and 26 + 24 = 50. - _Michael Somos_, Jan 21 2015

G. E. Andrews, MacMahon's Partition Analysis II: Fundamental Theorems, Annals Combinatorics, 4 (2000), 327-338.
G. E. Andrews, P. Paule and A. Riese, MacMahon's Partition Analysis VIII: Plane partition diamonds, Advances Applied Math., 27 (2001), 231-242 (Cor. 2.1, n=1).
S. P. Humphries, Braid groups, infinite Lie algebras of Cartan type and rings of invariants, Topology and its Applications, 95 (3) (1999) pp. 173-205.

T. D. Noe, Table of n, a(n) for n = 0..1000
Nesrine Benyahia-Tani, Zahra Yahi, and Sadek Bouroubi, Ordered and non-ordered non-congruent convex quadrilaterals inscribed in a regular n-gon, Rostocker Math. Kolloq. 68, 71-79 (2013), Theorem 5.
W. Duke, On codes and Siegel modular forms, Int. Math. Res. Notes 1993, No. 5, Theorem 2. [MR1219862 (94d:11029)]
Michele Graffeo, Sergej Monavari, Riccardo Moschetti, and Andrea T. Ricolfi, Enumeration of partitions via socle reduction, arXiv:2501.10267 [math.CO], 2025. See p. 40.
S. P. Humphries, Home page
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 450
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 232
G. Nebe, E. M. Rains and N. J. A. Sloane, Self-Dual Codes and Invariant Theory, Springer, Berlin, 2006.
Michael Somos, In the Elliptic Realm
Index entries for linear recurrences with constant coefficients, signature (1,2,-1,-2,-1,2,1,-1).

See A266769 for a version without the four leading zeros.
Cf. A001993, A070557, A070558, A070559, A089299, A001970, A089292, A026810, A001400.
First differences of A097701.
Cf. A082424, A082437.

a:=[0,0,0,0,1,1,3,4];; for n in [9..60] do a[n]:=a[n-1]+2*a[n-2]-a[n-3]-2*a[n-4]-a[n-5]+2*a[n-6]+a[n-7]-a[n-8]; od; a; # G. C. Greubel, Sep 10 2019

K:=Rationals(); M:=MatrixAlgebra(K,4); q1:=DiagonalMatrix(M,[1,-1,1,-1]); p1:=DiagonalMatrix(M,[1,1,-1,-1]); q2:=DiagonalMatrix(M,[1,1,1,-1]); h:=M![1,1,1,1, 1,1,-1,-1, 1,-1,1,-1, 1,-1,-1,1]/2; H:=MatrixGroup<4,K|q1,q2,h,p1>; MolienSeries(H);

R:=PowerSeriesRing(Integers(), 60); [0,0,0,0] cat Coefficients(R!( x^4/((1-x)*(1-x^2)^2*(1-x^3)) )); // G. C. Greubel, Sep 10 2019

a:= n-> (Matrix(8, (i,j)-> if (i=j-1) then 1 elif j=1 then [1,2,-1,-2,-1,2,1,-1][i] else 0 fi)^n)[1,5]: seq(a(n), n=0..60); # Alois P. Heinz, Jul 31 2008

CoefficientList[Series[x^4/((1-x)*(1-x^2)^2*(1-x^3)), {x,0,60}], x] (* Jean-François Alcover, Mar 30 2011 *)
LinearRecurrence[{1,2,-1,-2,-1,2,1,-1},{0,0,0,0,1,1,3,4},60] (* Harvey P. Dale, Mar 04 2012 *)
a[ n_]:= Quotient[9(n+1)(-1)^n +2n^3 -9n +65, 144]; (* Michael Somos, Jan 21 2015 *)
a[ n_]:= Sign[n] SeriesCoefficient[ x^4/((1-x)(1-x^2)^2(1-x^3)), {x, 0, Abs@n}]; (* Michael Somos, Jan 21 2015 *)

{a(n) = (9*(n+1)*(-1)^n + 2*n^3 - 9*n + 65) \ 144}; /* Michael Somos, Jan 21 2015 */

a(n)=([0,1,0,0,0,0,0,0; 0,0,1,0,0,0,0,0; 0,0,0,1,0,0,0,0; 0,0,0,0,1,0,0,0; 0,0,0,0,0,1,0,0; 0,0,0,0,0,0,1,0; 0,0,0,0,0,0,0,1; -1,1,2,-1,-2,-1,2,1]^n*[0;0;0;0;1;1;3;4])[1,1] \\ Charles R Greathouse IV, Feb 06 2017

def AA008763_list(prec):
    P. = PowerSeriesRing(ZZ, prec)
    return P(x^4/((1-x)*(1-x^2)^2*(1-x^3))).list()
AA008763_list(60) # G. C. Greubel, Sep 10 2019

Let f4(n) = number of partitions n = p+q+r+s into exactly 4 parts, with p >= q >= r >= s >= 1 (see A026810, A001400) and let g4(n) be the number with q > r (so that g4(n) = f4(n-2)). Then a(n) = f4(n) + g4(n).
a(n) = (1/144)*( 2*n^3 + 9*n*((-1)^n - 1) - 16*((n is 2 mod 3) - (n is 1 mod 3)) ).
a(n) = (1/72)*(n+3)*(n+2)*(n+1)-(1/12)*(n+2)*(n+1)+(5/144)*(n+1)+(1/16)*(n+1)*(-1)^n+(1/16)*(-1)^(n+1)+(7/144)+(2*sqrt(3)/27)*sin(2*Pi*n/3). - Richard Choulet, Nov 27 2008
a(n) = a(n-1) + 2*a(n-2) - a(n-3) - 2*a(n-4) - a(n-5) + 2*a(n-6) + a(n-7) - a(n-8), n>7. - Harvey P. Dale, Mar 04 2012
a(n) = floor((9*(n+1)*(-1)^n + 2*n^3 - 9*n + 65)/144). - Tani Akinari, Nov 06 2012
a(n+1) - a(n) = A008731(n-3). - R. J. Mathar, Aug 06 2013
a(n) = -a(-n) for all n in Z. - Michael Somos, Jan 21 2015
Euler transform of length 3 sequence [1, 2, 1]. - Michael Somos, Jun 26 2017

Entry revised Dec 25 2003

A164680 Expansion of x/((1-x)^3(1-x^2)^3(1-x^3)).

Original entry on oeis.org

1, 3, 9, 20, 42, 78, 139, 231, 372, 573, 861, 1254, 1791, 2499, 3432, 4629, 6162, 8085, 10492, 13455, 17094, 21503, 26832, 33201, 40795, 49764, 60333, 72687, 87096, 103785, 123075, 145236, 170646, 199626, 232617, 269997, 312277, 359898, 413448, 473438

Offset: 1

Alford Arnold, Aug 21 2009

Convolution of A006918 with A001399, or of A002625 with A059841 (A000035 if offsets are respected),
or of A038163 with A022003 or of A057524 with A027656 or of A014125 with the aerated version of A000217,
or of A002624 with A103221, or of A002623 with A008731, or of other combinations of splitting the signature -/3,3,1 into two components.
If we apply the enumeration of Molien series as described in A139672,
this is row 45=9*5 of a table of values related to Molien series, i.e., the
product of the sequence on row 9 (A006918) with the sequence on row 5 (A001399).
This is associated with the root system E6, and can be described using the additive function on the affine E6 diagram:
      1
      |
      2
      |
1--2--3--2--1

			To calculate a(3), we consider the first three terms of A001399 = (1 1 2...)
and the first three terms of A006918 = (1 2 5 ...), to get the convolved a(3) = 1*5+1*2+2*1 = 9.

G. C. Greubel, Table of n, a(n) for n = 1..1000
Index entries for Molien series
Index entries for linear recurrences with constant coefficients, signature (3,0,-7,3,6,0,-6,-3,7,0,-3,1).

Cf. A139672 (row 21).
For G2, the corresponding sequence is A001399.
For F4, the corresponding sequence is A115264.
For E7, the corresponding sequence is A210068.
For E8, the corresponding sequence is A045513.
See A210634 for a closely related sequence.

R:=PowerSeriesRing(Integers(), 40); Coefficients(R!( 1/((1-x)^3*(1-x^2)^3*(1-x^3)) )); // G. C. Greubel, Jan 13 2020

seq(coeff(series(x/((1-x)^3*(1-x^2)^3*(1-x^3)), x, n+1), x, n), n = 1..40); # G. C. Greubel, Jan 13 2020

Rest@CoefficientList[Series[x/((1-x)^3*(1-x^2)^3*(1-x^3)), {x,0,40}], x] (* G. C. Greubel, Jan 13 2020 *)
LinearRecurrence[{3,0,-7,3,6,0,-6,-3,7,0,-3,1},{1,3,9,20,42,78,139,231,372,573,861,1254},40] (* Harvey P. Dale, Aug 03 2025 *)

Vec(1/(1-x)^3/(1-x^2)^3/(1-x^3)+O(x^40)) \\ Charles R Greathouse IV, Sep 23 2012

x=PowerSeriesRing(QQ, 'x', 40).gen()
1/((1-x)^3*(1-x^2)^3*(1-x^3))

a(n) = round( -(-1)^n*(n+3)*(n+7)/256 +(6*n^6 +180*n^5 +2070*n^4 +11400*n^3 +30429*n^2 +34290*n +9785)/103680 ) - R. J. Mathar, Mar 19 2012

Edited and extended by R. J. Mathar, Aug 22 2009

Corrected link to index entries - R. J. Mathar, Aug 26 2009

A114209 Number of permutations of [n] having exactly two fixed points and avoiding the patterns 123 and 231.

Original entry on oeis.org

0, 1, 0, 2, 1, 3, 2, 5, 3, 7, 5, 9, 7, 12, 9, 15, 12, 18, 15, 22, 18, 26, 22, 30, 26, 35, 30, 40, 35, 45, 40, 51, 45, 57, 51, 63, 57, 70, 63, 77, 70, 84, 77, 92, 84, 100, 92, 108, 100, 117, 108, 126, 117, 135, 126, 145, 135, 155, 145, 165, 155, 176, 165, 187, 176, 198, 187

Offset: 1

Emeric Deutsch, Nov 17 2005

nonn

			a(2)=1 because we have 12; a(3)=0 because no permutation of [3] can have exactly two fixed points; a(4)=2 because we have 1432 and 3214.

T. Mansour and A. Robertson, Refined restricted permutations avoiding subsets of patterns of length three, Annals of Combinatorics, 6, 2002, 407-418.
Index entries for linear recurrences with constant coefficients, signature (0,2,1,-1,-2,0,1).

Cf. A008731, A114208, A114210.

a:=proc(n) if n mod 6 = 0 then n*(n+6)/24 elif n mod 6 = 1 or n mod 6 = 5 then (n^2-1)/24 elif n mod 6 = 2 or n mod 6 = 4 then (n+2)*(n+4)/24 else (n^2-9)/24 fi end: seq(a(n),n=1..70);

a(n) = n(n+6)/24 if n mod 6 = 0; (n^2-1)/24 if n mod 6 = 1 or 5; (n+2)(n+4)/24 if n mod 6 = 2 or 4; (n^2-9)/24 if n mod 6 = 3.

a(n) = A008731(n-2). O.g.f.: x^2/((1-x)^3(1+x)^2(1+x+x^2)). [R. J. Mathar, Aug 11 2008]

cp's OEIS Frontend

A008763 Expansion of g.f.: x^4/((1-x)*(1-x^2)^2*(1-x^3)).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

GAP

Magma

Magma

Maple

Mathematica

PARI

PARI

Sage

Formula

Extensions

A164680 Expansion of x/((1-x)^3*(1-x^2)^3*(1-x^3)).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Sage

Formula

Extensions

A114209 Number of permutations of [n] having exactly two fixed points and avoiding the patterns 123 and 231.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Formula

A008763 Expansion of g.f.: x^4/((1-x)(1-x^2)^2(1-x^3)).

A164680 Expansion of x/((1-x)^3(1-x^2)^3(1-x^3)).