Arxiv:Math/0609596V1 [Math.CO] 21 Sep 2006 29Suharc -Al Ngs@Ukzn.Ac.Za E-Mail: Africa South 3209 Number Ento 1

Home , Binomial coefficient

arXiv:math/0609596v1 [math.CO] 21 Sep 2006 29SuhArc -al [email protected] e-mail: Africa South 3209 number ento 1. Definition natural following the across came we [7] In [4]). (Cf interpretations In fsltost h altpolm n eie that verifies One problem. ballot the ( to solutions to of origin the from paths of ihu rsigtedaoa.I sgvnby given is It diagonal. the crossing without ya By endi [5].) in defined ubro ah from paths of number n sdntdby denoted is and a oss fhrzna rvria oiieui movements. unit positive vertical or horizontal of consist Fraohrdrvto,oecncekthat check can one derivation, another (For i (1 + drs:Sho fMteaia cecs nvriyo Kw of University Sciences, Mathematical of School Address: 1 ete oeo opoignwiette n formulas. and identities new In proving to on move then We h lsia aaa ubri h ubro ah rmteoii t origin the from paths of number the is number Catalan classical The o opeees egv in give we completeness, For § upre yNra-Agat2067063. grant Norway-SA by Supported edrv oebnma dniisb onigCtlnnmesadt and numbers Catalan counting by identities binomial some derive we 4 § path ecnie oevrain of variations some consider we 3 , utvraegnrtn ucin ipeepeso for expression simple derived. A is determinant r close function. a New generating and considered. relations multivariate is recursive numbers identities, binomial Catalan clude the of generalization A )or 0) C OEIETTE N OMLSINVOLVING FORMULAS AND IDENTITIES SOME ema nt sequence finite a mean we , n,m e words Key a a nacethsoy n[]Brrn studied Bertrand [1] In history. ancient an has i +1 Let EEAIE AAA NUMBERS CATALAN GENERALIZED C = n,m 0 a ≤ aaa ubr,pts eeaigfunctions generating paths, numbers, Catalan : (0 i . (0 + m , − ≤ 2 , m 1) ,m n, n. ) . 1. to hti osy akfrom walk a say, to is That h eeaie aaa ubri enda the as defined is number Catalan generalized The § ro of proof a 2 ihu rsigtedaoa,a h number the as diagonal, the crossing without ) Introduction ( n i-hNg Siu-Ah − { C a n,m ,n m, 1 i } 0 ≤ n opt lsdfrua o them. for formulas closed compute and n − i ≤ 1 + C 1 n m C n,m n,m where ) ihu rsigteline the crossing without 2 = n C = n n,m d Zl-aa,Pietermaritzburg, aZulu-Natal, (1 a n n a aywell-known many has and i − 2 + = ∈ m 2 m m oml o the for formula Z a 1 + a )(1+ 0 2 n 1 + h Catalan the + a to n either and 1 m,n n n,m slsin- esults − a m generalization: − n n , ) m 2 + hs steps whose h number the o the for n n , m ( o a . y i ≥ +1 ,n n, The = heir d m. qk = x ) 2 variants. In 5 we obtain recursive relations for the generalized Catalan numbers, § from which identities for the central binomial coefficients and an identity for certain sum of products of the generalized Catalan numbers are given. In the main section 6 we give a close formula for the multivariate generating function for the § generalized Catalan numbers. In 7 more formulas are produced, including certain § product of the generating functions. We end in the final section 8 with an elegant § formula for the Catalan determinant.

2. Generalized Catalan numbers

For completeness, we give a direct derivation of a close formula for n,m although C as remarked earlier it is already included in references such as [1] and [5].

Theorem 2. Given integers 0 m n, we have ≤ ≤ 2m +1 2n (1) n,m = . C n + m +1 n + m Proof. Fix 0 m n. Welet Λ bethesetofpathsfrom (0, 2m) to (n m,n m) ≤ ≤ − − − and Γ be the subset of “bad” paths, i.e. those that cross the line y = x. So

n,m = Λ Γ . C | | − | |

Given a path P Γ, let k be smallest such that P crosses y = x at (k, k), i.e. both ∈ (k, k) and (k, k +1) belong to P. Let Q be the path obtained from P by reflecting across the line y = x + 1 only the portion of P from (0, 2m) to (k, k + 1). The − result is a path from ( 2m 1, 1) to (n m,n m). See Figure 1 above. This − − − − is simply the André’s reflection method and one can easily see that it establishes a bijection between Γ and the set of paths from ( 2m 1, 1) to (n m,n m). − − − − Therefore 2(n m) ( 2m 1+1) 2n Γ = − − − − = . | | (n m) 1 n m 1 − − − − On the other hand, 2(n m) ( 2m) 2n 2n Λ = − − − = = . | | n m n m n + m − − Hence 2n 2n n,m = C n + m − n m 1 − − 3 Hn-m, n-mL

Hk, k+1L

Q P

H-2m-1, 1L O

y= x+1

H0, -2mL

Figure 1. Portion of P reﬂected across y = x + 1 to form Q. and the conclusion follows by noticing that 2n n m 2n n m 2n = − = − . n m 1 n + m +1 n m n + m +1 n + m − − −

Bertrand’s ballot problem: Suppose the ﬁnal result of an election is such that party A gets 2m more votes than its opposition B, say A gets n + m and B gets n m votes. What is the probability that A is at least 2m votes ahead of B in all − partial counts throughout the whole process?

Corollary 3. The probability of the event in the generalized ballot problem is 2m +1 . n + m +1

Proof. By shifting all paths 2m units up, it is clear that n,m is the number of C paths from (0, 0) to (n m,n + m) without crossing y = x +2m. Then identify − each vote for party A with (0, 1) and that for party B with (1, 0). The conclusion follows by noticing that each possible ballot outcome in the problem corresponds to a path from (0, 0) to (n m,n + m), hence there are − 4

2n of them in total, and those ballot outcome with A leading B by at n + m least 2m votes corresponds to paths from (0, 0) to (n m,n + m) without crossing − y = x +2m.

3. Some variations

Superﬁcially, the deﬁnition of n,m seems to be too restricted, since it considers C only the even number 2m and allowing touching y = x. Here in this section we consider some variations and they can be computed in similar ways.

Deﬁnition 4. Let 0 m n. ≤ ≤ ¯n,m := the number of paths from (0, 2m 1) to (n m,n m) without crossing C − − − − y = x.

n,m := the number of paths from (0, 2m) to (n m,n m) strictly below D − − − y = x until reaching (n m,n m). − − ¯n,m := the number of paths from (0, 2m 1) to (n m,n m) strictly below D − − − − y = x until reaching (n m,n m). − −

Theorem 5. Suppose 0 m n. Then ≤ ≤ 2m +2 2n +1 (1) ¯n,m = n ,m = ; C D +1 +1 n + m +2 n m − 2m 2n 1 (2) n,m = − ; D n + m n m − 2m +1 2n (3) ¯n,m = n,m = . D C n + m +1 n m − Proof. We first prove (2). Notice that a path P from the set defining n,m must D pass (n m,n m 1) before reaching (n m,n m). Similar to the proof of − − − − − Theorem 2, P can be reflected across the line y = x to form a path Q from ( 2m, 0) − to (n m,n m 1), as shown in Figure 2. − − − 5

Hn-m, n-mL Hn-m, n-m-1L

Q O P H-2m, 0 L

y= x-1 y= x H0, -2mL

Figure 2. Reﬂection across y = x.

We let Λ be the set of paths from (0, 2m) to (n m,n m 1) and Γ be − − − − the subset of “bad” paths, i.e. those that cross the line y = x 1. Clearly, 2n 1 − Λ = − . By the reﬂection trick used in the proof of Theorem 2, we can | | n m identify paths− in Γ with paths from ( 2m, 0) to (n m,n m 1). Therefore 2n 1 − − − − Γ = − . | | n + m 2n 1 2n 1 2m 2n 1 So n,m = Λ Γ = − − = − . D | | − | | n m − n + m n + m n m − −

To prove (1), by shifting 1 unit down, we can identify paths deﬁning ¯n,m with C paths from (0, 2m 2) to (n m,n m) strictly below y = x until reaching − − − − (n m,n m). Therefore − − ¯n,m = n ,m C D +1 +1 and the conclusion follows.

Now for (3), we notice that by shifting 1 unit up, paths deﬁning ¯n,m are D identiﬁed with paths from (0, 2m) to (n m,n m) without crossing y = x. − − − Therefore ¯n,m = n,m. D C 6

4. Identities

In this section we derive some interesting identities which are corollaries of Theorem 2 or Theorem 5.

Corollary 6. n− 1 1 2k +1 2(n k) 2 2n (2) − = , 2k +1 k m n k 2m +1 n m 1 k m X= − − − − where 0 m n 1. ≤ ≤ −

Proof. A “bad” path P from (0, 2m) to (n m,n m), i.e. one that crosses − − − y = x, is a path from (0, 2m) to some (k, k), immediately before crossing y = x, − then proceed to (k, k +1) and takes any path to (n m,n m) as shown in Figure 3 − − below.

Hn-m, n-mL

Hk, k+1L Hk, kL

O P

y= x

H0, -2mL

Figure 3. Decomposition of a “bad” path P. 7

The number of paths from (0, 2m) to (k, k) = ((m + k) m, (m + k) m) − − − without crossing y = x is given by m+k,m, while the number of paths from (k, k+1) 2(nC m k) 1 to (n m,n m) is given by − − − . Hence the number of the “bad” − − n m k paths are − −

n−m− 1 2(n m k) 1 m k,m − − − C + n m k k X=0 − − n−1 2(n k) 1 = k,m − − C n k k X=0 − n−1 2m +1 2k 2(n k) 1 = − − . k + m +1 k m n k k X=0 − − 2n Since this number is also n,m the conclusion follows from the last n m −C part of the proof of Theorem 2. −

Corollary 7. n 1 2k 2(n k) 2m +1 2n (3) − = k(n k + 1) k m n k m(n + m + 1) n m k m X= − − − − where 1 m n 1. ≤ ≤ − Proof. Note that a path from (0, 2m) to (n m,n m) without crossing the line − − − y = x can be decomposed into a path from (0, 2m) to the ﬁrst touch at y = x, − some (k, k), 0 k n m, followed by a path from (k, k), to (n m,n m) ≤ ≤ − − − without crossing the line y = x, as shown in Figure 4 below. 8

Hn-m, n-mL

Hk, kL

y= x

H0, -2mL

Figure 4. The decomposition at (k, k).

The number of the former paths is k m,m and the number for the latter is D + n−m−k, , therefore we have C 0

n−m (4) n,m = k m,m n−m−k, . C D + C 0 k X=0 Hence, by Theorem 2 and Theorem 5,

n−m 2m +1 2n 2m 2(k + m) 1 1 2(n m k) = − − − n + m +1 n m k +2m k n m k +1 n m k − k=0 − − − − Xn 2m 2k 1 1 2(n k) = − − k + m k m n k +1 n k k=m − − − Xn m 2k 2(n k) = − , k(n k + 1) k m n k k m X= − − − and equation (3) follows. 9

Remark 8. It is easy to check that equation (2) can be transformed into the following form: n 2n 2k +1 2k 2(n k) (5) − − (k + m + 1)(n k + 1) k m n k k m X= − − − (2n + 1)(2n + 2) 2n = . (2m + 1)(n + m + 1)(n + m + 2) n m − Although equations (3) and (5) are similar, it is not clear how they are related.

Corollary 9. n−m k 2h + k 1 2n 2h k 2n (6) − − − = h + k h n m h n m h X=0 − − − where 0 m n 1 and 1 k 2m. ≤ ≤ − ≤ ≤ 2n Proof. The number of paths from (0, 2m) to (n m,n m) is , the − − − n m right side of the above equation (6). − For each 1 k 2m we can decompose these paths into the path from (0, 2m) ≤ ≤ − ﬁrst touching the line y = x 2m + k at (h,h 2m + k), for some 0 h n m, − − ≤ ≤ − followed by a path from (h,h 2m + k) to (n m,n m), as shown in Figure 5 − − − below. Hn-m, n-mL

y= x-2m+k

Hh, h-2m+kL O y= x

H0, -2mL

Figure 5. The decomposition at (h,h 2m + k). − 10

2n 2h k The number of paths from (h,h 2m + k) to (n m,n m) is − − . − − − n m h By shifting up 2m k units, we can see that a path from (0, 2m) ﬁrst− touching− − − the line y = x 2m + k at (h,h 2m + k) corresponds to a path from (0, k) ﬁrst − − − touching the line y = x at (h,h). Therefore, if k is even, the number is given by k 2h + k 1 h+ k , k = − ; D 2 2 h + k h if k is odd, the number is given by k 2h + k 1 ¯ k− k− h+ 1 , 1 = − ; D 2 2 h + k h both producing the same number. k 2h + k 1 2n 2h k Summing the product − − − over 0 h n m, h + k h n m h ≤ ≤ − equation (6) follows. − −

Remark 10. It appears that equation (6) actually holds for all k> 0. The current proof doesn’t seem to work for k> 2m.

5. Recursive relations and identities

In this section we will derive two recursive relations for the generalized Catalan numbers. As corollaries, we obtain some identities for the central binomial coeﬃ- cients. The ﬁrst recursive relation for the generalized Catalan numbers we consider is the following.

Theorem 11. For 0 m < n, we have the following recursive relation ≤ n−m−1 n(2m + 1) (7) n,m = k m,m n−m−k− , . C (n m)(m + 1) C + C 1 0 k − X=0 Proof. A path deﬁning n,m is either a path touching the line y = x only once at C the terminal point (n m,n m) (the number of such is given by n,m ) or a path − − D touching some (k, k), where 0 k n m 1 for the last time before hitting the ≤ ≤ − − terminal point (n m,n m), as shown in Figure 6. − − 11 Hn-m, n-mL Hn-m, n-m-1L

Hk, kL Hk+1, kL

O P

y= x

H0, -2mL Figure 6. A path P touching y = x at (k, k) for the last time before touching the terminal point (n m,n m). − −

Write (k, k) as ((k + m) m, (k + m) m), we have − − n−m−1 n,m = n,m + k m,m n−m−k− , . C D C + C 1 0 k X=0 But 2m 2n 1 m(n + m + 1) n,m = − = n,m, D n + m n m n(2m + 1) C − hence the result follows.

As a consequence, we have an identity for the nth central binomial coeﬃcient.

Corollary 12. Let 0 m < n, then ≤ m n−m−1 2n n + m +1 n + k (8) = k m,m n−m−k− , . n m +1 n k C + C 1 0 k ! k Y=0 − X=0 Proof. By equation (7) in Theorem 11 we have: n−m− 2m +1 (2n)! n(2m + 1) 1 = k m,m n−m−k− , , n + m +1 (n + m)! (n m)! (n m)(m + 1) C + C 1 0 k − − X=0 i.e. n−m− (2n)! n + m +1 n 1 = k m,m n−m−k− , , (n + m)! (n m)! m +1 n m C + C 1 0 k − − X=0 12 hence the result follows if m = 0; for m> 0, we get n−m−1 (2n)! n + m +1 n(n + 1) (n + m) = ··· k m,m n−m−k− , , n!2 m +1 n(n 1) (n m + 1)(n m) C + C 1 0 k − ··· − − X=0 and equation (8) is proved.

The following is an identity to be used in 7. §

Corollary 13. Let 1 m < n, then ≤ n−m− m n 1 m +1 (9) k m,m n−m−k− , = (n k) (n + k) C + C 1 0 n! − k k ! k m ! X=0 Y=1 =Y+2 m n n k n + k = − . k k k ! k m ! Y=1 =Y+2 Proof. This is just another form of equation (8) in Corollary 12 by noticing that n 2n n + k = . n k k Y=1

Now we consider another recursive relation for the generalized Catalan numbers.

Theorem 14. Let 0 m n. Then ≤ ≤ n m k + m +1 (10) n,m = k,m n−k, . C 2m +1 k C C 0 k m X= Proof. From Theorem 5 in [7] we can obtain m(n + m + 1) n,m = n,m. D n(2m + 1) C The result then follows from the identity (4) of Corollary 7 in [7].

From the above we have another expression for the nth central binomial coeﬃ- cient. The proof is similar to that of Corollary 12.

Corollary 15. Let 1 m n, then ≤ ≤ m+1 n 2n m(n m) n + k k + m +1 (11) = − k,m n−k, . n (2m + 1)2 n k +1 k C C 0 k ! k m Y=1 − X= 13

6. Generating functions

In this section we derive a formula for the generating function of generalized Catalan numbers. A comprehensive teatment of generating functions can be found in [3].

For ﬁxed m 0 we let γm(x) denote the generating function of n m,m, i.e. ≥ C + the power series ∞ n (12) γm(x) = n m,m x . C + n=0 X It is well-known for the classical Catalan numbers n, that C 0 1 √1 4x γ (x) = − − . 0 2x (See for example [4].) We denote the corresponding multivariate generating function by: ∞ n m (13) Γ(x, y) = n m,m x y . C + n,m=0 X In [2] a formula for the generating function for an,m is given, i.e. 1 yγ (xy) a(x, y) := a xnym = − 0 . n,m 1 x y n≥m≥ X 0 − − Using the relation n,m = an m,n−m, it can be veriﬁed that C + 1 x x Γ(x, y) = a √y, + a √y, . 2 √y − −√y However we have another derivation here:

Theorem 16. γ (x) x(1 √1 4x ) x, y 0 . (14) Γ( ) = 2 i.e. − − 1 y (γ0(x)) 2x(x + y) y(1 √1 4x ) − − − −

We ﬁrst need the following lemmas:

Lemma 17. Let n 0 then ≥

(15) n m ,m = a m,m b, c, , C + +1 +1 C + C 0 C 0 a b c n + X+ = where we sum over a, b, c 0. ≥ 14

Proof. We count paths deﬁning n m ,m , i.e. paths from (0, 2(m + 1)) to C + +1 +1 − (n,n) = (n + m +1 (m + 1),n + m +1 (m + 1)) without crossing y = x, by − − using their ﬁrst touch (a,a 1) and (a + b,a + b) at the lines y = x 1 and y = x − − respectively as in Figure 7 below.

H n, n L

H a+b, a+b L

H a, a-1 L

y= x O

y= x-1

H 0, -2m-2 L

Figure 7. A path defining n+m+1,m+1 whose first touch at y = x 1 is (a,a 1) and whoseC first touch at y = x is (a + b,a + b) for− some a,b − 0 and a + b n. ≥ ≤

The number of paths from (0, 2m 2) that ﬁrst touch y = x 1 at (a,a 1) − − − − is the same as those reaching (a,a 2) without crossing y = x 2, i.e. the number − − is a m,m. Similarly, the number of paths from (a,a 1) ﬁrst reaching y = x at C + − (a + b,a + b) is b, and the number of paths from (a + b,a + b) to (n,n) without C 0 crossing y = x is c, , where c = n a b, hence the result is proved. C 0 − −

Now we have a recursive relation for the generating function γm(x): 15

Lemma 18. For m 0, we have ≥ 2 (16) γm+1(x) = γm(x) (γ0(x)) .

Proof. We use Lemma 17: ∞ n γm (x) = n m ,m x +1 C + +1 +1 n=0 X∞ n = a m,m b, c, x C + C 0 C 0 n=0 a b c n X +X+ = ∞ ∞ 2 n n = n m,m x n, x C + C 0 n=0 ! n=0 ! X 2 X = γm(x) (γ0(x)) .

Now we prove Theorem 16 from Lemma 18. Note that: ∞ 2 2 m+1 y Γ(x, y)(γ0(x)) = γm(x) (γ0(x)) y m=0 X∞ m+1 = γm+1(x) y m=0 X = Γ(x, y) γ (x), − 0 from which the result follows.

7. More formulas

In this section we prove a few more fomulas relating to the generating functions.

Theorem 19. Let m 1. The generating function γm(x) is equal to ≥

∞ m n+m+1 2x n + m k +1 n + m + k +1 1+ − xn . 1 √1 4x k k n=1 k ! k m ! ! − − X Y=1 =Y+2 16

Proof. Replacing n m 1 by n in Corollary 13 equation (9), we obtain for n 1 − − ≥ that n m n+m+1 n + m k +1 n + m + k +1 k m,m n−k, = − . C + C 0 k k k k ! k m ! X=0 Y=1 =Y+2 Then the theorem follows from ∞ n γm(x) γ (x) = k m,m h, x 0 C + C 0 n=0 k h n X +X= ∞ n n = k m,m n−k, x C + C 0 n=0 k ! X X=0 ∞ m n+m+1 n + m k +1 n + m + k +1 = 1+ − xn. k k n=1 k ! k m ! X Y=1 =Y+2

We now the product γm(x) γ0(x) and the following power series: ∞ (n + 1)(m + 1) n θ(x) = n m ,m x . (n + m + 1)(2m + 1) C + +1 n=0 X

Theorem 20. Let 0 m n. Then ≤ ≤ 1 √1 4x θ(x) = γ (x) γ (x) = γ (x) − − . m 0 m 2x

Proof. From Theorem 11 equation (7) we have: (n + m + 1)(2m + 1) n m ,m = k m,m h, , C + +1 (n + 1)(m + 1) C + C 0 k h n +X= therefore ∞ n 1 √1 4x θ(x) = k m,m h, x = γ(x) − − . C + C 0 2x n=0 k h n ! X +X=

Corollary 21. For 1 m n, we have ≤ ≤ m n+m (n + m)(2m + 1) n + m k n + m + k (17) n m,m = − . C + n(m + 1) k k k ! k m ! Y=1 =Y+2

Proof. Compare the power series expansion of γm(x) γ0(x) given in the last line of the proof of Theorem 19 and the one given by θ(x). Then equation (17) is obtained by replacing n +1 for n. 17

8. A beautiful determinant

In this ﬁnal section, we prove the following:

Theorem 22. For integer N 0, ≥

(2i + 1) (2(i + j))! N(N+1) det = 2 2 . j! (2i + j + 1)! 0≤i,j≤N

We suggest the name Catalan determinant for the above, for the entries of the matrix are the generalized Catalan numbers i j, i. C + First we need a slight modiﬁcation of [8] Theorem 1. (See also [6].) An almost identical proof is given here for completeness.

2 Lemma 23. Let f(x) = 1+ a1x + a2x + ... be a formal power series and let j 2i+1 j j ci,j = [x ]f (x), where [x ] denotes the coeﬃcient of x in the series. Then

N(N+1) 2 det [ci,j ]0≤i,j≤N = (2a1) .

Proof. Let

i+j i C := [ci,j ] and B := [bi,j ] where bi,j = ( 1) . 0≤i,j≤N 0≤i,j≤N − j The (i, k)-entry of BC is given by N N i k j i 2j+1 bi,j cj,k = ( 1) [x ] ( 1) f (x) − − j j=0 j=0 X X i i i = (( 1)i [xk] ( 1)j f 2j+1(x) (since = 0 for j>i) − − j j j=0 X i = ( 1)i [xk] 1 f 2(x) f(x) − − i = ( 1)i [xk] ( 2a x + ... ) (1 + a x + ... ) , − − 1 1 i which equals (2a1) if k = i and equals 0 if k < i, hence BC is an upper i triangular matrix with diagonal entries (2a1) . In particular N N N i ( +1) det BC = (2a1) = (2a1) 2 . i=0 Y But B is lower triangular with diagonal entries 1, so det B = 1 and the result follows. 18

To prove Theorem 22, we let, as in 6, γi(x) be the generating function of i j, i, § C + i.e. the power series ∞ j γi(x) = i j, i x . C + n=0 X Then 1 √1 4x γ (x) = − − = 1+ x +2x2 +5x3 + 14x4 + ... 0 2x and Lemma 18 gives that

2i+1 γi(x) = γ (x) for i 0. 0 ≥ Therefore j j 2i+1 i j, i = [x ]γi(x) = [x ]γ (x) C + 0 and Theorem 22 is now proved by applying Lemma 23 to f(x) = γ0(x).

References

1. J. Bertrand, Solution d’un probleme, Comptes Rendus Acad´emie des Sciences, Paris, 105 (1887), 369. 2. L. Carlitz & J. Riordan, Two element lattice permutation numbers and their q-generalization, Duke Mathematical journal, 31 (1964), 317-388. 3. P. Flajolet & R. Sedgewick, Analytic Combinatorics. 4. P. Hilton & J. Pederson, Catalan Numbers, Their Generalization, and Their Uses, Mathemat- ical Intelligencer, 13(2) (1991), 64-75. 5. P. Hilton & J. Pederson, A generalization of the Catalan numbers, Congressus Numerantium, 165(2003), 5-18. 6. C. Krattenthaler, Advanced determinant calculus, The Andrews Festschrift (Maratea, 1998). Sm. Lothar. Combin. 42 (1999), Art. B42q, 67 pp. (electronic). 7. S.-A. Ng, Margin Trading through Hyper Timeline, to appear in International Journal of Theoretical and Applied Finance. 8. H. Wilf, A combinatorial determinant, J. Combin. Theory Ser. A, in press. ( http://arXiv.org/abs/math/9809120 ) (1998).