OFFSET
8,3
LINKS
G. C. Greubel, Table of n, a(n) for n = 8..60
Saibal Mitra and Bernard Nienhuis, Osculating Random Walks on Cylinders, in Discrete Random Walks, DRW'03, Cyril Banderier and Christian Krattenthaler (eds.), Discrete Mathematics and Theoretical Computer Science Proceedings AC, pp. 259-264.
FORMULA
Even n: Q(n, m) = C_{n/2-m}(n) + Sum_{r=1..(n-2*m)/4} (-1)^r * ((m+2*r)/(m+r)) * binomial(m+r, r) * C_{n/2 - m - 2*r}(n).
Odd n: Q(n, m) = Sum_{r=0..(n-2*m-1)/4)} (-1)^r * binomial(m+r,r) * ( C_{(n-1)/2 - m - 2*r}(n) - C_{(n-1)/2 - m - 2*r - 1}(n) ), where the c_{k}(n) are the absolute values of the coefficients of the characteristic polynomial of the n X n Pascal matrix P_{i, j} = binomial(i+j-2, i-1). The sequence is given by Q(n, 4).
MATHEMATICA
M[n_, k_]:= Table[Binomial[i+j-2, i-1], {i, n}, {j, k}];
c[k_, n_]:= Coefficient[CharacteristicPolynomial[M[n, n], x], x, k]//Abs;
Q[n_?EvenQ, m_]:= Sum[(-1)^r*((m+2*r)/(m+r))*Binomial[m +r, r]*c[n/2 -m-2*r, n], {r, 0, (n-2*m)/4}];
Q[n_?OddQ, m_]:= Sum[(-1)^r*Binomial[m+r, r]*(c[(n-1)/2 -m-2*r, n] - c[(n-1)/2 -m-2*r-1, n]), {r, 0, (n-2*m-1)/4}];
Table[Q[n, 4], {n, 8, 26}] (* G. C. Greubel, Nov 15 2019 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Saibal Mitra (smitra(AT)zonnet.nl), Mar 20 2004
STATUS
approved