Delannoy number

From HandWiki
Short description: Number of paths between grid corners, allowing diagonal steps
Delannoy number
Named afterHenri–Auguste Delannoy
No. of known termsinfinity
Formula[math]\displaystyle{ D(m,n) = \sum_{k=0}^{\min(m,n)} \binom{m}{k} \binom{n}{k} 2^k }[/math]
OEIS index

In mathematics, a Delannoy number [math]\displaystyle{ D }[/math] describes the number of paths from the southwest corner (0, 0) of a rectangular grid to the northeast corner (m, n), using only single steps north, northeast, or east. The Delannoy numbers are named after French army officer and amateur mathematician Henri Delannoy.[1]

The Delannoy number [math]\displaystyle{ D(m,n) }[/math] also counts the number of global alignments of two sequences of lengths [math]\displaystyle{ m }[/math] and [math]\displaystyle{ n }[/math],[2] the number of points in an m-dimensional integer lattice or cross polytope which are at most n steps from the origin,[3] and, in cellular automata, the number of cells in an m-dimensional von Neumann neighborhood of radius n[4] while the number of cells on a surface of an m-dimensional von Neumann neighborhood of radius n is given with (sequence A266213 in the OEIS).

Example

The Delannoy number D(3,3) equals 63. The following figure illustrates the 63 Delannoy paths from (0, 0) to (3, 3):

Delannoy3x3.svg

The subset of paths that do not rise above the SW–NE diagonal are counted by a related family of numbers, the Schröder numbers.

Delannoy array

The Delannoy array is an infinite matrix of the Delannoy numbers:[5]

 m
n 
0 1 2 3 4 5 6 7 8
0 1 1 1 1 1 1 1 1 1
1 1 3 5 7 9 11 13 15 17
2 1 5 13 25 41 61 85 113 145
3 1 7 25 63 129 231 377 575 833
4 1 9 41 129 321 681 1289 2241 3649
5 1 11 61 231 681 1683 3653 7183 13073
6 1 13 85 377 1289 3653 8989 19825 40081
7 1 15 113 575 2241 7183 19825 48639 108545
8 1 17 145 833 3649 13073 40081 108545 265729
9 1 19 181 1159 5641 22363 75517 224143 598417

In this array, the numbers in the first row are all one, the numbers in the second row are the odd numbers, the numbers in the third row are the centered square numbers, and the numbers in the fourth row are the centered octahedral numbers. Alternatively, the same numbers can be arranged in a triangular array resembling Pascal's triangle, also called the tribonacci triangle,[6] in which each number is the sum of the three numbers above it:

            1
          1   1
        1   3   1
      1   5   5   1
    1   7  13   7   1
  1   9  25  25   9   1
1  11  41  63  41  11   1

Central Delannoy numbers

The central Delannoy numbers D(n) = D(n,n) are the numbers for a square n × n grid. The first few central Delannoy numbers (starting with n=0) are:

1, 3, 13, 63, 321, 1683, 8989, 48639, 265729, ... (sequence A001850 in the OEIS).

Computation

Delannoy numbers

For [math]\displaystyle{ k }[/math] diagonal (i.e. northeast) steps, there must be [math]\displaystyle{ m-k }[/math] steps in the [math]\displaystyle{ x }[/math] direction and [math]\displaystyle{ n-k }[/math] steps in the [math]\displaystyle{ y }[/math] direction in order to reach the point [math]\displaystyle{ (m, n) }[/math]; as these steps can be performed in any order, the number of such paths is given by the multinomial coefficient [math]\displaystyle{ \binom{m+n-k}{k , m-k , n-k} = \binom{m+n-k}{m} \binom{m}{k} }[/math]. Hence, one gets the closed-form expression

[math]\displaystyle{ D(m,n) = \sum_{k=0}^{\min(m,n)} \binom{m+n-k}{m} \binom{m}{k} . }[/math]

An alternative expression is given by

[math]\displaystyle{ D(m,n) = \sum_{k=0}^{\min(m,n)} \binom{m}{k} \binom{n}{k} 2^k }[/math]

or by the infinite series

[math]\displaystyle{ D(m,n) = \sum_{k=0}^\infty \frac{1}{2^{k+1}} \binom{k}{n} \binom{k}{m}. }[/math]

And also

[math]\displaystyle{ D(m,n) = \sum_{k=0}^{n} A(m,k), }[/math]

where [math]\displaystyle{ A(m,k) }[/math] is given with (sequence A266213 in the OEIS).

The basic recurrence relation for the Delannoy numbers is easily seen to be

[math]\displaystyle{ D(m,n)=\begin{cases}1 &\text{if }m=0\text{ or }n=0\\D(m-1,n) + D(m-1,n-1) + D(m,n-1)&\text{otherwise}\end{cases} }[/math]

This recurrence relation also leads directly to the generating function

[math]\displaystyle{ \sum_{m,n = 0}^\infty D(m, n) x^m y^n = (1 - x - y - xy)^{-1} . }[/math]

Central Delannoy numbers

Substituting [math]\displaystyle{ m = n }[/math] in the first closed form expression above, replacing [math]\displaystyle{ k \leftrightarrow n-k }[/math], and a little algebra, gives

[math]\displaystyle{ D(n) = \sum_{k=0}^n \binom{n}{k} \binom{n+k}{k} , }[/math]

while the second expression above yields

[math]\displaystyle{ D(n) = \sum_{k=0}^n \binom{n}{k}^2 2^k . }[/math]

The central Delannoy numbers satisfy also a three-term recurrence relationship among themselves,[7]

[math]\displaystyle{ n D(n) = 3(2n-1)D(n-1) - (n-1)D(n-2) , }[/math]

and have a generating function

[math]\displaystyle{ \sum_{n = 0}^\infty D(n) x^n = (1-6x+x^2)^{-1/2} . }[/math]

The leading asymptotic behavior of the central Delannoy numbers is given by

[math]\displaystyle{ D(n) = \frac{c \, \alpha^n}{\sqrt{n}} \, (1 + O(n^{-1})) }[/math]

where [math]\displaystyle{ \alpha = 3 + 2 \sqrt{2} \approx 5.828 }[/math] and [math]\displaystyle{ c = (4 \pi (3 \sqrt{2} - 4))^{-1/2} \approx 0.5727 }[/math].

See also

References

  1. Banderier, Cyril; Schwer, Sylviane (2005), "Why Delannoy numbers?", Journal of Statistical Planning and Inference 135 (1): 40–54, doi:10.1016/j.jspi.2005.02.004 
  2. Covington, Michael A. (2004), "The number of distinct alignments of two strings", Journal of Quantitative Linguistics 11 (3): 173–182, doi:10.1080/0929617042000314921 
  3. Luther, Sebastian; Mertens, Stephan (2011), "Counting lattice animals in high dimensions", Journal of Statistical Mechanics: Theory and Experiment 2011 (9): P09026, doi:10.1088/1742-5468/2011/09/P09026, Bibcode2011JSMTE..09..026L, http://stacks.iop.org/1742-5468/2011/i=09/a=P09026 
  4. Breukelaar, R.; Bäck, Th. (2005), "Using a Genetic Algorithm to Evolve Behavior in Multi Dimensional Cellular Automata: Emergence of Behavior", Proceedings of the 7th Annual Conference on Genetic and Evolutionary Computation (GECCO '05), New York, NY, USA: ACM, pp. 107–114, doi:10.1145/1068009.1068024, ISBN 1-59593-010-8 
  5. Sulanke, Robert A. (2003), "Objects counted by the central Delannoy numbers", Journal of Integer Sequences 6 (1): Article 03.1.5, Bibcode2003JIntS...6...15S, http://www.emis.de/journals/JIS/VOL6/Sulanke/delannoy.pdf 
  6. Sloane, N. J. A., ed. "Sequence A008288 (Square array of Delannoy numbers D(i,j) (i >= 0, j >= 0) read by antidiagonals)". OEIS Foundation. https://oeis.org/A008288. 
  7. Peart, Paul; Woan, Wen-Jin (2002). "A bijective proof of the Delannoy recurrence". Congressus Numerantium 158: 29–33. ISSN 0384-9864. 

External links