cp's OEIS Frontend

A326368 Number of tilings of an equilateral triangle of side length n with unit triangles (of side length 1) and exactly three unit "lozenges" or "diamonds" (also of side length 1).

%I A326368 #30 Sep 18 2019 04:58:47
%S A326368 0,0,18,434,2814,11127,33365,83568,184254,369254,686952,1203930,
%T A326368 2009018,3217749,4977219,7471352,10926570,15617868,21875294,30090834,
%U A326368 40725702,54318035,71490993,92961264,119547974,152182002,191915700,239933018,297560034,366275889
%N A326368 Number of tilings of an equilateral triangle of side length n with unit triangles (of side length 1) and exactly three unit "lozenges" or "diamonds" (also of side length 1).
%H A326368 Colin Barker, <a href="/A326368/b326368.txt">Table of n, a(n) for n = 1..1000</a>
%H A326368 Richard J. Mathar, <a href="https://arxiv.org/abs/1909.06336">Lozenge tilings of the equilateral triangle</a>, arXiv:1909.06336 [math.CO], 2019.
%H A326368 <a href="/index/Rec#order_07">Index entries for linear recurrences with constant coefficients</a>, signature (7,-21,35,-35,21,-7,1).
%F A326368 a(n) = (1/16)*(n-2)*(9*n^5 - 9*n^4 - 81*n^3 + 81*n^2 + 160*n - 192) for n >= 2 (proved by _Greg Dresden_ and E. Sijaric).
%F A326368 From _Colin Barker_, Jul 02 2019: (Start)
%F A326368 G.f.: x^3*(18 + 308*x + 154*x^2 - 87*x^3 + 10*x^4 + 2*x^5) / (1 - x)^7.
%F A326368 a(n) = 7*a(n-1) - 21*a(n-2) + 35*a(n-3) - 35*a(n-4) + 21*a(n-5) - 7*a(n-6) + a(n-7) for n>8.
%F A326368 (End)
%e A326368 We can represent a unit triangle this way:
%e A326368        o
%e A326368       / \
%e A326368      o - o
%e A326368 and a unit "lozenge" or "diamond" has these three orientations:
%e A326368      o
%e A326368     / \          o - o            o - o
%e A326368    o   o  and   /   /   and also   \   \
%e A326368     \ /        o - o                o - o
%e A326368      o
%e A326368 and for n=3, here is one of the 18 different tiling of the triangle of side length 3 with exactly three lozenges:
%e A326368           o
%e A326368          / \
%e A326368         o   o
%e A326368        / \ / \
%e A326368       o - o   o
%e A326368      /   / \ / \
%e A326368     o - o - o - o
%t A326368 Rest@ CoefficientList[Series[x^3*(18 + 308 x + 154 x^2 - 87 x^3 + 10 x^4 + 2 x^5)/(1 - x)^7, {x, 0, 30}], x] (* _Michael De Vlieger_, Jul 07 2019 *)
%o A326368 (PARI) concat([0,0], Vec(x^3*(18 + 308*x + 154*x^2 - 87*x^3 + 10*x^4 + 2*x^5) / (1 - x)^7 + O(x^40))) \\ _Colin Barker_, Jul 02 2019
%Y A326368 Cf. A273464, A326367, A326369.
%K A326368 nonn,easy
%O A326368 1,3
%A A326368 _Greg Dresden_, Jul 01 2019