|
|
A167546
|
|
The ED1 array read by antidiagonals
|
|
17
|
|
|
1, 1, 1, 2, 4, 1, 6, 12, 7, 1, 24, 48, 32, 10, 1, 120, 240, 160, 62, 13, 1, 720, 1440, 960, 384, 102, 16, 1, 5040, 10080, 6720, 2688, 762, 152, 19, 1, 40320, 80640, 53760, 21504, 6144, 1336, 212, 22, 1
(list;
table;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,4
|
|
COMMENTS
|
The coefficients in the upper right triangle of the ED1 array (m > n) were found with the a(n,m) formula while the coefficients in the lower left triangle of the ED1 array (m <= n) were found with the recurrence relation, see below. We use for the array rows the letter n (>= 1) and for the array columns the letter m (>= 1).
Our procedure for finding the coefficients in the lower left triangle can be compared with the procedure that De Smit and Lenstra used to fill in the hole in the middle of 'The Print Gallery' by M. C. Escher, see the links. In this lithograph Escher made use of the so-called Droste effect, hence we propose to call this square array of numbers the ED1 array.
|
|
LINKS
|
Johannes W. Meijer, The four Escher-Droste arrays, jpg image, Mar 08 2013.
|
|
FORMULA
|
a(n,m) = (2*(m-1)!/(m-n-1)!)*Integral_{y>=0} sinh(y*(2*n-1))/cosh(y)^(2*m-1) for m > n.
The (n-1)-differences of the n-th array row lead to the recurrence relation
Sum_{k=0..n-1} (-1)^k*binomial(n-1,k)*a(n,m-k) = (2*n-1)*(n-1)!
which in its turn leads to, see also A167557,
a(n,m) = 4^(m-1)*(m-1)!*(n+m-2)!/(2*m-2)! for m <= n.
|
|
EXAMPLE
|
The ED1 array begins with:
1, 1, 1, 1, 1, 1, 1, 1, 1, 1
1, 4, 7, 10, 13, 16, 19, 22, 25, 28
2, 12, 32, 62, 102, 152, 212, 282, 362, 452
6, 48, 160, 384, 762, 1336, 2148, 3240, 4654, 6432
24, 240, 960, 2688, 6144, 12264, 22200, 37320, 59208, 89664
120, 1440, 6720, 21504, 55296, 122880, 245640, 452880, 783144, 1285536
|
|
MAPLE
|
nmax:=10; mmax:=10; for n from 1 to nmax do for m from 1 to n do a(n, m) := 4^(m-1)*(m-1)!*(n-1+m-1)!/(2*m-2)! od; for m from n+1 to mmax do a(n, m):= (2*n-1)*(n-1)! + sum((-1)^(k-1)*binomial(n-1, k)*a(n, m-k), k=1..n-1) od; od: for n from 1 to nmax do for m from 1 to n do d(n, m):=a(n-m+1, m) od: od: T:=1: for n from 1 to nmax do for m from 1 to n do a(T):= d(n, m): T:=T+1: od: od: seq(a(n), n=1..T-1);
|
|
MATHEMATICA
|
nmax = 10; mmax = 10; For[n = 1, n <= nmax, n++, For[m = 1, m <= n, m++, a[n, m] = 4^(m - 1)*(m - 1)!*((n - 1 + m - 1)!/(2*m - 2)!)]; For[m = n + 1, m <= mmax, m++, a[n, m] = (2*n - 1)*(n - 1)! + Sum[(-1)^(k - 1)*Binomial[n - 1, k]*a[n, m - k], {k, 1, n - 1}]]; ]; For[n = 1, n <= nmax, n++, For[m = 1, m <= n, m++, d[n, m] = a[n - m + 1, m]]; ]; t = 1; For[n = 1, n <= nmax, n++, For[m = 1, m <= n, m++, a[t] = d[n, m]; t = t + 1]]; Table[a[n], {n, 1, t - 1}] (* Jean-François Alcover, Dec 20 2011, translated from Maple *)
|
|
CROSSREFS
|
A000142 equals the first column of the array.
A167550 equals the a(n, n+1) diagonal of the array.
A047053 equals the a(n, n) diagonal of the array.
A167558 equals the a(n+1, n) diagonal of the array.
A167551 equals the row sums of the ED1 array read by antidiagonals.
A167552 is a triangle related to the a(n) formulas of rows of the ED1 array.
A167556 is a triangle related to the GF(z) formulas of the rows of the ED1 array.
A167557 is the lower left triangle of the ED1 array.
Cf. A068424 (the (m-1)!/(m-n-1)! factor), A007680 (the (2*n-1)*(n-1)! factor).
|
|
KEYWORD
|
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|