OFFSET
0,3
COMMENTS
There seem to be several typos in Jepsen's equations. The enumeration here is derived from the expression of p(n) as 1/8ths of Psi(e)+2*Psi(rho)+Psi(rho^2)+2*Psi(sigma)+2*Psi(rho*sigma) if n>=3.
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Charles H. Jepsen, Packing a box with bricks, Math. Mag. 64 (2) (1991) 92-97, Table 1.
Index entries for linear recurrences with constant coefficients, signature (6,-4,-26,33,8,-8,24,-31,-14,12,2,-1).
FORMULA
G.f.: 1 +x +2*x^2 -x^3*(-7 +16*x +57*x^2 -118*x^3 -38*x^4 +30*x^5 -53*x^6 +127*x^7 +42*x^8 -49*x^9 -7*x^10 +4*x^11) / ( (x-1)*(1+x) *(x^2+2*x-1) *(x^2+1) *(x^2-4*x+1) *(x^4-4*x^2+1) ).
MATHEMATICA
CoefficientList[Series[1+x+2*x^2-x^3*(-7+16*x+57*x^2-118*x^3-38*x^4+30*x^5-53*x^6+127*x^7+42*x^8-49*x^9-7*x^10+4*x^11)/((x-1)*(1+x)*(x^2+2*x-1)*(x^2+1)*(x^2-4*x+1)*(x^4-4*x^2+1)), {x, 0, 33}], x] (* Vincenzo Librandi, Sep 02 2025 *)
PROG
(Magma) m:=35; R<x>:=PowerSeriesRing(Integers(), m); Coefficients(R!(1 +x +2*x^2 -x^3*(-7 +16*x +57*x^2 -118*x^3 -38*x^4 +30*x^5 -53*x^6 +127*x^7 +42*x^8 -49*x^9 -7*x^10 +4*x^11) / ( (x-1)*(1+x) *(x^2+2*x-1) *(x^2+1) *(x^2-4*x+1) *(x^4-4*x^2+1)) )); // Vincenzo Librandi, Sep 02 2025
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
R. J. Mathar, Aug 31 2025
STATUS
approved
