OFFSET
0,2
COMMENTS
This is the case k = 2 of the parametric family A_{3n+1}(k) = (k^n-1)^2*k^(n+1)+2*k^n*(k^n - 1)*(k^(n+1)-1), which counts axis-aligned 2 X 2 squares in the 3D (orthogonal) projection of a d-dimensional hypercube with d = 3n+1. Thus a(n) = (2^n-1)^2*2^(n+1) +2*2^n*(2^n-1)*(2^(n+1)-1) is the number of 2 X 2 tiles in the 3D projection of the 4, 7, 10, ...-dimensional hypercubes (side 2).
LINKS
Salvatore Ferraro, Counting axis-aligned 2X2 tiles in the 3D projection of hypercubes (short note), Zenodo, 2025
Index entries for linear recurrences with constant coefficients, signature (14,-56,64).
FORMULA
a(n) = 3*2^(3*n+1) - 5*2^(2*n+1) + 2^(n+2).
From Elmo R. Oliveira, Mar 20 2026: (Start)
G.f.: -16*x*(x + 1)/((8*x - 1)*(4*x - 1)*(2*x - 1)).
E.g.f.: 2*exp(2*x)*(2 - 5*exp(2*x) + 3*exp(6*x)).
a(n) = 14*a(n-1) - 56*a(n-2) + 64*a(n-3) for n > 2. (End)
MAPLE
a := n -> (2^n-1)^2*2^(n+1) +2*2^n*(2^n-1)*(2^(n+1)-1);
seq(a(n), n=0..39);
MATHEMATICA
Table[(2^n - 1)^2*2^(n + 1) + 2*2^n*(2^n - 1)*(2^(n + 1) - 1), {n, 0, 39}]
PROG
(Python)
def a(n:int)->int:
return (2**n-1)**2*2**(n+1) + 2*2**n*(2**n-1)*(2**(n+1) - 1)
print([a(n) for n in range(40)])
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Salvatore Ferraro, Sep 27 2025
STATUS
approved
