Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.
%I #23 Sep 07 2021 04:28:15
%S 1,2,4,3,7,17,4,10,26,68,5,13,35,95,259,6,16,44,122,340,950,7,19,53,
%T 149,421,1193,3387,8,22,62,176,502,1436,4116,11814,9,25,71,203,583,
%U 1679,4845,14001,40503,10,28,80,230,664,1922,5574,16188,47064,136946,11,31,89,257,745,2165,6303,18375,53625,156629,457795
%N Triangle I(m,n) read by rows: number of perfect lattice paths on the m*n board.
%H D. Yaqubi, M. Farrokhi D. G., and H. Ghasemian Zoeram, <a href="https://arxiv.org/abs/1612.08697">Lattice paths inside a table, I</a>, arXiv:1612.08697 [math.CO], 2016-2017, array I(m,n).
%F I(m,n) = (n+2)*3^(n-2) + (m-n)*Sum_{i=0..n-1} A005773(i)*A005773(n-i) + 2*Sum_{k=0..n-3) (n-k-2)*3^(n-k-3)*A001006(k). [Yaqubi Corr. 2.10]
%F I(m,n) = A188866(m-1,n) for m > 1. - _Pontus von Brömssen_, Sep 06 2021
%e Triangle begins:
%e 1;
%e 2, 4;
%e 3, 7, 17;
%e 4, 10, 26, 68;
%e 5, 13, 35, 95, 259;
%e 6, 16, 44, 122, 340, 950;
%e 7, 19, 53, 149, 421, 1193, 3387;
%e 8, 22, 62, 176, 502, 1436, 4116, 11814;
%e 9, 25, 71, 203, 583, 1679, 4845, 14001, 40503;
%e 10, 28, 80, 230, 664, 1922, 5574, 16188, 47064, 136946;
%p Inm := proc(n,m)
%p if m >= n then
%p (n+2)*3^(n-2)+(m-n)*add(A005773(i)*A005773(n-i),i=0..n-1)
%p +2*add((n-k-2)*3^(n-k-3)*A001006(k),k=0..n-3) ;
%p else
%p 0 ;
%p end if;
%p end proc:
%p for m from 1 to 13 do
%p for n from 1 to m do
%p printf("%a,",Inm(n,m)) ;
%p end do:
%p printf("\n") ;
%p end do:
%p # Second program:
%p A296449row := proc(n) local gf, ser;
%p gf := n -> 1 + (x*(n - (3*n + 2)*x) + (2*x^2)*(1 +
%p ChebyshevU(n - 1, (1 - x)/(2*x))) / ChebyshevU(n, (1 - x)/(2*x)))/(1 - 3*x)^2;
%p ser := n -> series(expand(gf(n)), x, n + 1);
%p seq(coeff(ser(n), x, k), k = 1..n) end:
%p for n from 0 to 11 do A296449row(n) od; # _Peter Luschny_, Sep 07 2021
%t (* b = A005773 *) b[0] = 1; b[n_] := Sum[k/n*Sum[Binomial[n, j] * Binomial[j, 2*j - n - k], {j, 0, n}], {k, 1, n}];
%t (* c = A001006 *) c[0] = 1; c[n_] := c[n] = c[n-1] + Sum[c[k] * c[n-2-k], {k, 0, n-2}];
%t Inm[n_, m_] := If[m >= n, (n + 2)*3^(n - 2) + (m - n)*Sum[b[i]*b[n - i], {i, 0, n - 1}] + 2*Sum[(n - k - 2)*3^(n - k - 3)*c[k], {k, 0, n-3}], 0];
%t Table[Inm[n, m], {m, 1, 13}, {n, 1, m}] // Flatten (* _Jean-François Alcover_, Jan 23 2018, adapted from first Maple program. *)
%Y Cf. A081113 (diagonal), A000079 (2nd row), A001333 (3rd row), A126358, A057960, A126360, A002714, A126362, A188866.
%K nonn,tabl,easy
%O 1,2
%A _R. J. Mathar_, Dec 13 2017