Exercises from Section 1.2.3

Tord M. Johnson

February 12, 2014

1. [10] The text says that a1 + a2 + ⋯ + a0 = 0. What then, is a2 + ⋯ + a0?

The identities

∑ 1≤k≤0ak = 0 = a1 + ∑ 2≤k≤0ak

imply that

∑ 2≤k≤0ak = a2 + ⋯ + a0 = −a1.

2. [01] What does the notation ∑ ⁡ 1≤j≤naj mean, if n = 3.14?

The notation ∑ ⁡ 1≤j≤naj if n = 3.14 represents the sum of all aj such that 1 ≤ j ≤ n = 3.14, or equivalently, j ∈{1,2,3}. That is,

∑ 1≤j≤naj = a1 + a2 + a3.

▸ 3. [13] Without using the ∑ ⁡ -notation, write out the equivalent of

∑ 0≤n≤5 1 2n + 1

and also the equivalent of

∑ 0≤n2≤5 1 2n2 + 1.

Explain why the two results are different, in spite of rule (b).

For the first notation,

∑ 0≤n≤5 1 2n + 1 = 1 1 + 1 3 + 1 5 + 1 7 + 1 9 + 1 11,

while for the second notation,

∑ 0≤n2≤5 1 2n2 + 1 = 1 9 + 1 3 + 1 1 + 1 3 + 1 9.

The two results are different, in spite of rule (b), as P(n) = n2 is not a permutation, since P is neither one-to-one ( − 12 = 12 serves as a counterexample) nor onto (there is no integer n such that n2 = 3).

4. [10] Without using the ∑ ⁡ -notation, write out the equivalent of each side of Eq. (10) as a sum of sums for the case n = 3.

For n = 3, the equivalent of the left side of Eq. (10) is

∑ i=13 ∑ j=1ia ij = (a11) + (a21 + a22) + (a31 + a32 + a33)

and of the right side is

∑ j=13 ∑ i=j3 = (a 11 + a21 + a31) + (a22 + a32) + (a33).

▸ 5. [HM20] Prove that rule (a) is valid for arbitrary infinite series, provided the series converge.

Proposition. The distributive law, for the products of sums (∑ ⁡ R(i)ai)(∑ ⁡ S(j)bj) = ∑ ⁡ R(i)(∑ ⁡ S(j)aibj) holds for arbitrary infinite series, provided the sums converge.

Proof. Assume that we have two infinite series ai, bj whose sums converge; that is, that

∑ R(i)ai =(lim ⁡ n→∞∑ R(i) 0≤i<n ai)+(lim ⁡ n→∞∑ R(i) −n≤i<0 ai) = α

and

∑ S(j)bj =(lim ⁡ n→∞∑ S(j) 0≤j<n bj)+(lim ⁡ n→∞∑ S(j) −n≤j<0 bj) = β

for arbitary limits α, β. We must show that

(∑ R(i)ai)(∑ S(j)bj) = ∑ R(i)(∑ S(j)aibj).

But

(∑ R(i)ai)(∑ S(j)bj) = lim ⁡ n→∞((∑ R(i) 0≤i<n ai + ∑ R(i) −n≤i<0 ai)(∑ S(j) 0≤j<n bj + ∑ S(j) −n≤j<0 bj)) = lim ⁡ n→∞((∑ R(i) 0≤i<n ai)(∑ S(j) 0≤j<n bj + ∑ S(j) −n≤j<0 bj) +(∑ R(i) −n≤i<0 ai)(∑ S(j) 0≤j<n bj + ∑ S(j) −n≤j<0 bj)) = lim ⁡ n→∞(∑ R(i) 0≤i<n ai(∑ S(j) 0≤j<n bj + ∑ S(j) −n≤j<0 bj) + ∑ R(i) −n≤i<0 ai(∑ S(j) 0≤j<n bj + ∑ S(j) −n≤j<0 bj)) = lim ⁡ n→∞(∑ R(i) 0≤i<n (∑ S(j) 0≤j<n aibj + ∑ S(j) −n≤j<0 aibj) + ∑ R(i) −n≤i<0 (∑ S(j) 0≤j<n aibj + ∑ S(j) −n≤j<0 aibj)) = lim ⁡ n→∞(∑ R(i) 0≤i<n (∑ S(j)aibj) + ∑ R(i) −n≤i<0 (∑ S(j)aibj)) = ∑ R(i)(∑ S(j)aibj)

as we needed to show. □

6. [HM20] Prove that rule (d) is valid for an arbitrary infinite series, provided that any three of the four sums exist.

Proposition. Manipulating the domain, ∑ ⁡ R(j)aj + ∑ ⁡ S(j)aj = ∑ ⁡ R(j)∨S(j)aj + ∑ ⁡ R(j)∧S(j)aj is valid for arbitrary infinite series, provided that any three of the four sums exist.

Proof. Assume that R and S are arbitrary infinite series. We must show that if any three of the sums ∑ ⁡ R(j)aj, ∑ ⁡ S(j)aj, ∑ ⁡ R(j)∨S(j)aj, ∑ ⁡ R(j)∧S(j)aj exist, then we may manipulate the domain so that ∑ ⁡ R(j)aj + ∑ ⁡ S(j)aj = ∑ ⁡ R(j)∨S(j)aj + ∑ ⁡ R(j)∧S(j)aj. It is sufficient to show that the convergence of three sums is sufficient for the convergence of the fourth.

Case 1. [∑ ⁡ R(j)aj, ∑ ⁡ S(j)aj, ∑ ⁡ R(j)∨S(j)aj converge.] We must show that ∑ ⁡ R(j)∧S(j)aj converges. But if the sums ∑ ⁡ R(j)aj, ∑ ⁡ S(j)aj exist, then clearly their conjunction ∑ ⁡ R(j)∧S(j)aj exists.

Case 2. [∑ ⁡ R(j)aj, ∑ ⁡ S(j)aj, ∑ ⁡ R(j)∧S(j)aj converge.] We must sow that ∑ ⁡ R(j)∨S(j)aj converges. But if the sums ∑ ⁡ R(j)aj, ∑ ⁡ S(j)aj exist, then clearly their disjunction ∑ ⁡ R(j)∨S(j)aj exists.

Case 3. [∑ ⁡ R(j)aj, ∑ ⁡ R(j)∨S(j)aj, ∑ ⁡ R(j)∧S(j)aj converge.] We must show that ∑ ⁡ S(j)aj converges. But if the conjunction ∑ ⁡ R(j)∧S(j)aj exists, then clearly the single sum ∑ ⁡ S(j)aj exists.

Case 4. [∑ ⁡ S(j)aj, ∑ ⁡ R(j)∨S(j)aj, ∑ ⁡ R(j)∧S(j)aj converge.] We must show that ∑ ⁡ R(j)aj converges. But if the conjunction ∑ ⁡ R(j)∧S(j)aj exists, then clearly the single sum ∑ ⁡ R(j)aj exists.

Therefore, in all cases, we have shown that the convergence of three sums is sufficient for the convergence of the fourth. □

7. [HM23] Given that c is an integer, show that ∑ ⁡ R(j)aj = ∑ ⁡ R(c−j)ac−j, even if both series are infinite.

Proposition. ∑ ⁡ R(j)aj = ∑ ⁡ R(c−j)ac−j, even if both series are infinite.

Proof. Assume that R is an arbitrary relation. We must show that ∑ ⁡ R(j)aj = ∑ ⁡ R(c−j)ac−j, even if both series are infinite. Since R′(j) = c − j is a permutation of the integers, proving the finite case, we must show the infinite case. But

∑ R(j)aj =(lim ⁡ n→∞∑ R(j) 0≤j<n aj) +(lim ⁡ n→∞∑ R(j) −n≤j<0 aj) =(lim ⁡ n→∞∑ R(c−j) 0≤c−j<n ac−j) +(lim ⁡ n→∞∑ R(c−j) −n≤c−j<0 ac−j) = ∑ R(c−j)ac−j

as we needed to show. □

8. [HM25] Find an example of infinite series in which Eq. (7) is false.

An example of an infinite series in which ∑ ⁡ R(i) ∑ ⁡ S(j)aij≠∑ ⁡ S(j) ∑ ⁡ R(j)aij is given by

aij = -1if i − 1 = j 1 if i + 1 = j 0 otherwise

for S(i) ≡ i ≥ 0, R(j) ≡ j ≥ 0. Then

∑ R(i) ∑ S(j)aij = lim ⁡ n→∞∑ 0≤i<n ∑ 0≤j<naij = lim ⁡ n→∞(∑ 0≤i<1 ∑ 0≤j<naij + ∑ 1≤i<n ∑ 0≤j<naij) = lim ⁡ n→∞(a01 + ∑ 1≤i<n ∑ 0≤j<naij) = lim ⁡ n→∞( − 1 + ∑ 1≤i<n ∑ 0≤j<naij) = −1

and

∑ S(j) ∑ R(i)aij = lim ⁡ n→∞∑ 0≤j<n ∑ 0≤i<naij = lim ⁡ n→∞(∑ 0≤j<1 ∑ 0≤i<naij + ∑ 1≤j<n ∑ 0≤i<naij) = lim ⁡ n→∞(a10 + ∑ 1≤j<n ∑ 0≤i<naij) = lim ⁡ n→∞(1 + ∑ 1≤j<n ∑ 0≤i<naij) = 1.

▸ 9. [05] Is the derivation of Eq. (14) valid even if n = −1?

No, the derivation for Eq. (14) is not valid if n = −1, since the application of rule (d) assumes n ≥ 0. That is, if n = −1

∑ 0≤j≤−1axj≠a + ∑ 1≤j≤−1axj

(depsite the fact that ∑ ⁡ 0≤j≤−1axj = 0 = a(1−x−1+1 1−x )).

10. [05] Is the derivation of Eq. (14) valid even if n = −2?

No, the derivation for Eq. (14) is not valid if n = −2, since the application of rule (d) assumes n ≥ 0. That is, if n = −2

∑ 0≤j≤−2axj≠a + ∑ 1≤j≤−2axj.

11. [03] What should the right-hand side of Eq. (14) be if x = 1?

In the case x = 1

∑ 0≤j≤naxj = ∑ 0≤j≤na(1)j = ∑ 0≤j≤na = (n + 1)a.

12. [10] What is 1 + 1 7 + 1 49 + 1 343 + ⋯ + (1 7)n?

The sum 1 + 1 7 + 1 49 + 1 343 + ⋯ + (1 7)n is simply a geometric progression whose closed form is given by Eq. (14) with a = 1 and x = 1 7. This yields

∑ 0≤j≤n(1) 1 7j = (1) 1 −1 7 n+1 1 −1 7 = 1 −1 7 n+1 6∕7 = 7 6 1 − 1∕7n+1 .

13. [10] Using Eq. (15) and assuming that m ≤ n, evaluate ∑ ⁡ j=mnj.

Since

∑ m≤j≤nj = ∑ 0≤j≤nj −∑ 0≤j≤m−1j

we can use Eq. (15) with a = 0 and b = 1 to evaluate each term on the right-hand side as

∑ m≤j≤nj = ∑ 0≤j≤nj −∑ 0≤j≤m−1j = 1 2n(n + 1) −1 2(m − 1)m = 1 2(n(n + 1) − m(m − 1)).

14. [11] Using the result of the previous exercise, evaluate ∑ ⁡ j=mn ∑ ⁡ k=rsjk.

Since ∑ ⁡ m≤j≤nj = 1 2(n(n + 1) − m(m − 1)), we can evaluate the sum as

∑ m≤j≤n ∑ r≤k≤sjk = ∑ m≤j≤nj ∑ r≤k≤sk = 1 2(n(n + 1) − m(m − 1)) 1 2(s(s + 1) − r(r − 1)) = 1 4 (n(n + 1) − m(m − 1))(s(s + 1) − r(r − 1)).

▸ 15. [M22] Compute the sum 1 × 2 + 2 × 22 + 3 × 23 + ⋯ + n × 2n for small values of n. Do you see the pattern developing in these numbers? If not, discover it by manipulations similar to those leading up to Eq. (14).

We can manipulate the sum as

∑ 0≤k≤nk2k = 0 + ∑ 1≤k≤nk2k = ∑ 1≤k≤nk2k = 2∑ 1≤k≤nk2k−1 = 2∑ 0≤k≤n−1(k + 1)2k = 2 ∑ 0≤k≤n−1k2k + ∑ 0≤k≤n−12k = 2 ∑ 0≤k≤nk2k − n2n + ∑ 0≤k≤n−12k = 2 ∑ 0≤k≤nk2k − n2n − (1 − 2n) = 2∑ 0≤k≤nk2k − n2n+1 − (2 − 2n+1).

Comparing the first with the last yields

∑ 0≤k≤nk2k = n2n+1 + 2 − 2n+1 = 2(n2n − 2n + 1).

16. [M22] Prove that

∑ j=0njxj = nxn+2 − (n + 1)xn+1 + x (x − 1)2 ,

if x≠1, without using mathematical induction.

Proposition. ∑ ⁡ 0≤j≤njxj = nxn+2−(n+1)xn+1+x (x−1)2 if x≠1.

Proof. Assume that x≠1, n ≥ 0. Then

∑ 0≤j≤njxj = 0 + ∑ 1≤j≤njxj = ∑ 1≤j≤njxj = x∑ 1≤j≤njxj−1 = x∑ 0≤j≤n−1(j + 1)xj = x ∑ 0≤j≤n−1jxj + ∑ 0≤j≤n−1xj = x ∑ 0≤j≤njxj − nxn + ∑ 0≤j≤n−1xj = x ∑ 0≤j≤njxj − nxn + 1 − xn 1 − x = x∑ 0≤j≤njxj − nxn+1 + x1 − xn 1 − x .

Comparing the first relation with the last, we have

(1 − x)∑ 0≤j≤njxj = −nxn+1 + x1 − xn 1 − x ;

hence we obtain the formula

∑ 0≤j≤njxj = − (1 − x)nxn+1 (1 − x)2 + x 1 − xn (1 − x)2 = − nxn+1 + nxn+2 + x − xn+1 (1 − x)2 = nxn+2 − (n + 1)xn+1 + x (1 − x)2

as we needed to show. □

▸ 17. [M00] Let S be a set of integers. What is ∑ ⁡ j∈S1?

∑ ⁡ j∈S1 is the cardinality—the number of elements—of S; that is,

∑ j∈S1 = |S|.

18. [M20] Show how to interchange the order of summation as in Eq. (9) given that R(i) is the relation “n is a multiple of i” and S(i,j) is the relation “1 ≤ j < i.”

In Eq. (9), we have

∑ R(i) ∑ S(i,j)aij = ∑ S′(j) ∑ R′(i,j)aij

for S′(j) ≡ (∃ ⁡i)(R(i) ∧ S(i,j)) and R′(i,j) ≡ R(i) ∧ S(i,j).

In the case that R(i) ≡ (∃ ⁡k)(ki = n) and S(i,j) ≡ 1 ≤ j < i, we find

S′(j) ≡ (∃ ⁡i)(R(i) ∧ S(i,j)) ≡ (∃ ⁡i)((∃ ⁡k)(ki = n) ∧ 1 ≤ j < i) ≡ 1 ≤ j < n for i ≤ n, k ≥ 1

and

R′(i,j) ≡ R(i) ∧ S(i,j) ≡ (∃ ⁡k)(ki = n) ∧ 1 ≤ j < i ≡ (∃ ⁡k)(ki = n) ∧ j < i;

that is, S′(j) the relation “1 ≤ j < n” and R′(i,j) the relation “n is a multiple of i and i > j.”

19. [20] What is ∑ ⁡ j=mn(aj − aj−1)?

We have that

∑ m≤j≤n(aj − aj−1) = ∑ m≤j≤naj −∑ m≤j≤naj−1 = ∑ m≤j≤naj −∑ m−1≤j≤n−1aj = ∑ m≤j≤naj − am−1 −∑ m≤j≤naj + an = an − am−1

assuming m ≤ n.

▸ 20. [25] Dr. I. J. Matrix has observed a remarkable sequence of formulas:

9 × 1 + 2 = 11,9 × 12 + 3 = 111,9 × 123 + 4 = 1111,9 × 1234 + 5 = 11111.
a.
Write the good doctor’s great discovery in terms of the ∑ ⁡ -notation.
b.
Your answer to part (a) undoubtedly involves the number 10 as a base of the decimal system; generalize this formula so that you get a formula that will perhaps work in any base b.
c.
Prove your formula from part (b) by using formulas derived in the text or in exercise 16 above.

The remarkable sequence of formulas observed by Dr. I. J. Matrix are analyzed below.

a.
The good doctor’s great discovery in terms of the ∑ ⁡ -notation is
(10 − 1)∑ 0≤k≤n(n − k)10k + n + 1 = ∑ 0≤k≤n10k.
b.
The above formula may be generalized for perhaps any base b as
(b − 1)∑ 0≤k≤n(n − k)bk + n + 1 = ∑ 0≤k≤nbk.
c.
We may prove the above formula.

Proposition. (b − 1)∑ ⁡ 0≤k≤n(n − k)bk + n + 1 = ∑ ⁡ 0≤k≤nbk.

Proof. Assume an arbitrary base b ≥ 2, and n ≥ 0. We must show that

(b − 1)∑ 0≤k≤n(n − k)bk + n + 1 = ∑ 0≤k≤nbk.

But

(b − 1)∑ 0≤k≤n(n − k)bk + n + 1 = (b − 1) n∑ 0≤k≤nbk −∑ 0≤k≤nkbk + n + 1 = (b − 1) n 1 − bn+1 1 − b −∑ 0≤k≤nkbk + n + 1 by Eq. (14) = (b − 1) n 1 − bn+1 1 − b −nbn+2 − (n + 1)bn+1 + b (b − 1)2 + n + 1by exercise 16 = n(bn+1 − 1) + b(bn(−nb + n + 1) − 1) b − 1 + n + 1 = bn+1 + n − b(n + 1) b − 1 + n + 1 = 1 − bn+1 1 − b = ∑ 0≤k≤nbk by Eq. (14)

as we needed to show. □

▸ 21. [M25] Derive rule (d) from (8) and (17).

We may derive rule (d) for manipulating the domain

∑ R(j)aj + ∑ S(j)aj = ∑ R(j)∨S(j)aj + ∑ R(j)∧S(j)aj

from Eq. (8)

∑ R(i)(bi + ci) = ∑ R(i)bi + ∑ R(i)ci

and Eq. (17)

∑ R(j)aj = ∑ jaj[R(j)].

Given that [p] + [q] = [p ∨ q] + [p ∧ q], as evidenced by the following table

[p][q][p] + [q][p ∨ q][p ∧ q][p ∨ q] + [p ∧ q]






000000
0 1 1 1 0 1
101101
1 1 2 1 1 2

we have

∑ R(j)aj + ∑ S(j)aj = ∑ jaj[R(j)] + ∑ jaj[S(j)] by Eq. (17) = ∑ j(aj[R(j)] + aj[S(j)]) by Eq. (8) = ∑ jaj([R(j)] + [S(j)]) = ∑ jaj([R(j) ∨ S(j)] + [R(j) ∧ S(j)]) = ∑ jaj[R(j) ∨ S(j)] + ∑ jaj[R(j) ∧ S(j)] by Eq. (8) = ∑ R(j)∨S(j)aj + ∑ R(j)∧S(j)aj. by Eq. (17)

▸ 22. [20] State the appropriate analogs of Eqs. (5), (7), (8), and (11) for products instead of sums.

We have the following analogs for products: change of variable:

∏ R(i)ai = ∏ R(j) = ∏ R(p(j))ap(j);

interchanging order of production:

∏ R(i) ∏ S(j)aij = ∏ S(j) ∏ R(i)aij;

a special case of the above:

∏ R(i)(bici) = ∏ R(i)bi ∏ R(i)ci ;

and manipulating the domain:

∏ R(j)aj ∏ S(j)aj = ∏ R(j)∨S(j)aj ∏ R(j)∧S(j)aj .

23. [10] Explain why it is a good idea to define ∑ ⁡ R(j)aj and ∏ ⁡ R(j)aj as zero and one, respectively, when no integers satisify R(j).

It is a good idea to define ∑ ⁡ j∈∅aj = 0 and ∏ ⁡ j∈∅aj = 1 as they are the identity elements for the operations of addition and multiplication, respectively. This way, ∑ ⁡ j∈∅aj + ∑ ⁡ R(j)aj = ∑ ⁡ R(j)aj and ∏ ⁡ j∈∅aj ∏ ⁡ R(j)aj = ∏ ⁡ R(j)aj.

24. [20] Suppose that R(j) is true for only finitely many j. By induction on the number of integers satisfying R(j), prove that log ⁡ b ∏ ⁡ R(j)aj = ∑ ⁡ R(j)(log ⁡ baj), assuming that all aj > 0.

Proposition. log ⁡ b ∏ ⁡ R(j)aj = ∑ ⁡ R(j)(log ⁡ baj) for all aj > 0.

Proof. Suppose aj > 0 for 0 ≤ j ≤ n. We must show that log ⁡ b ∏ ⁡ 0≤j≤naj = ∑ ⁡ 0≤j≤n(log ⁡ baj).

If n = 0, clearly log ⁡ b ∏ ⁡ 0≤j≤0aj = log ⁡ ba0 = ∑ ⁡ 0≤j≤0(log ⁡ baj). Then, assuming

log ⁡ b ∏ 0≤j≤kaj = ∑ 0≤j≤k(log ⁡ baj),

we must show that

log ⁡ b ∏ 0≤j≤k+1aj = ∑ 0≤j≤k+1(log ⁡ baj).

But

log ⁡ b ∏ 0≤j≤k+1aj = log ⁡ b ak+1 ∏ 0≤j≤kaj = log ⁡ b ∏ 0≤j≤kaj + log ⁡ bak+1 = ∑ 0≤j≤k(log ⁡ baj) + log ⁡ bak+1 = ∑ 0≤j≤k+1(log ⁡ baj)

as we needed to show. □

▸ 25. [15] Consider the following derivation; is anything amiss?

(∑ i=1na i)(∑ j=1n 1 aj) = ∑ 1≤i≤n ∑ 1≤j≤n ai aj = ∑ 1≤i≤n ∑ 1≤i≤nai ai = ∑ i=1n1 = n.

Yes. For one,

∑ 1≤i≤n ∑ 1≤j≤n ai aj≠∑ 1≤i≤n ∑ 1≤i≤nai ai;

and for another,

∑ 1≤i≤n ∑ 1≤i≤nai ai≠∑ i=1n1

(in fact, ∑ ⁡ 1≤i≤n ∑ ⁡ 1≤i≤nai ai = ∑ ⁡ i=1nn = n2).

26. [25] Show that ∏ ⁡ i=0n ∏ ⁡ j=0iaiaj may be expressed in terms of ∏ ⁡ i=0nai by manipulating the ∏ ⁡ -notation as stated in exercise 22.

Proposition. ∏ ⁡ 0≤i≤n ∏ ⁡ 0≤j≤iaiaj = ∏ ⁡ 0≤i≤nai n+2.

Proof. We must show that ∏ ⁡ 0≤i≤n ∏ ⁡ 0≤j≤iaiaj = ∏ ⁡ 0≤i≤nai n+2.

First note that

∏ 0≤i≤n ∏ 0≤j≤iaiaj = ∏ 0≤i≤n ∏ i≤j≤naiaj;

and then that

∏ 0≤i≤n ∏ 0≤j≤iaiaj 2 = ∏ 0≤i≤n ∏ 0≤j≤iaiaj ∏ 0≤i≤n ∏ 0≤j≤iaiaj = ∏ 0≤i≤n ∏ 0≤j≤iaiaj ∏ 0≤j≤iaiaj = ∏ 0≤i≤n ∏ 0≤j≤iaiaj ∏ i≤j≤naiaj = ∏ 0≤i≤n ∏ 0≤j≤naiaj aiai = ∏ 0≤i≤n ∏ 0≤j≤naiaj ∏ 0≤i≤nai2 = ∏ 0≤i≤n ain+1 ∏ 0≤j≤naj ∏ 0≤i≤nai 2 = ∏ 0≤i≤nai n+1 ∏ 0≤i≤n ∏ 0≤j≤naj ∏ 0≤i≤nai 2 = ∏ 0≤i≤nai n+1 ∏ 0≤j≤naj n+1 ∏ 0≤i≤nai 2 = ∏ 0≤i≤nai 2n+4.

Therefore,

∏ 0≤i≤n ∏ 0≤j≤iaiaj = ∏ 0≤i≤nai n+2

as we needed to show. □

27. [M20] Generalize the result of exercise 1.2.1-9 by proving that

∏ j=1n(1 − a j) ≥ 1 −∑ j=1na j,

assuming that 0 < aj < 1.

Proposition. ∏ ⁡ 1≤j≤n(1 − aj) ≥ 1 −∑ ⁡ 1≤j≤naj if 0 < aj < 1.

Proof. Let 0 < aj < 1 for 1 ≤ j ≤ n, n ≥ 0. We must show that

∏ 1≤j≤n(1 − aj) ≥ 1 −∑ 1≤j≤naj.

If n = 0, then clearly ∏ ⁡ 1≤j≤n(1 − aj) = 1 ≥ 1 = 1 −∑ ⁡ 1≤j≤naj. Then, assuming

∏ 1≤j≤k(1 − aj) ≥ 1 −∑ 1≤j≤kaj,

we must show that

∏ 1≤j≤k+1(1 − aj) ≥ 1 −∑ 1≤j≤k+1aj.

But since ak+1 ∑ ⁡ 0≤j≤kak ≥ 1,

∏ 1≤j≤k+1(1 − aj) = ∏ 1≤j≤k(1 − aj)(1 − ak+1) ≥ 1 −∑ 1≤j≤kaj (1 − ak+1) = 1 −∑ 1≤j≤kaj −ak+1 − ak+1 ∑ 0≤j≤kaj = 1 −∑ 0≤j≤kaj − ak+1 + ak+1 ∑ 0≤j≤kak ≥ 1 −∑ 0≤j≤kaj + ak+1 = 1 −∑ 1≤j≤k+1aj

as we needed to show. □

28. [M22] Find a simple formula for ∏ ⁡ j=2n(1 − 1∕j2).

We find that

∏ 2≤j≤n(1 − 1 j2) = ∏ 2≤j≤nj2 − 1 j2 = ∏ 2≤j≤n 1 j2 ∏ 2≤j≤n(j − 1)(j + 1) = ∏ 2≤j≤n1 j2 ∏ 1≤j≤n−1j ∏ 3≤j≤n+1j = 1 n!2(n − 1)!(n + 1)! 2 = n + 1 2n .

▸ 29. [M30] (a) Express ∑ ⁡ i=0n ∑ ⁡ j=0i ∑ ⁡ k=0jaiajak in terms of the multiple-sum notation explained at the end of the section. (b) Express the same sum in terms of ∑ ⁡ i=0nai, ∑ ⁡ i=0nai2, and ∑ ⁡ i=0nai3 [see Eq. (13)].

a.
We can express the sum in terms of the multiple-sum notation as
∑ 0≤i≤n ∑ 0≤j≤i ∑ 0≤k≤jaiajak = ∑ 0≤k≤j≤i≤naiajak.
b.
From Eq. (13) for arbitrary i > 0 we have
∑ 0≤j≤i−1 ∑ 0≤k≤jajak = 1 2 ∑ 0≤j≤i−1aj 2 + ∑ 0≤j≤i−1aj2

or equivalently

∑ 0≤j≤i−1aj 2 = 2∑ 0≤j≤i−1 ∑ 0≤k≤jajak −∑ 0≤j≤i−1aj2.

Also for arbitrary i > 0 we have that

∑ 0≤j≤iaj 3 = ∑ 0≤j≤i−1aj + ai 3 = ∑ 0≤j≤i−1aj 3 + 3a i ∑ 0≤j≤i−1aj 2 + 3a i2 ∑ 0≤j≤i−1aj + ai3 = ∑ 0≤j≤i−1aj 3 + 3a i 2∑ 0≤j≤i−1 ∑ 0≤k≤jajak −∑ 0≤j≤i−1aj2 + 3a i2 ∑ 0≤j≤i−1aj + ai3 = ∑ 0≤j≤i−1aj 3 + 6∑ 0≤j≤i−1 ∑ 0≤k≤jaiajak − 3∑ 0≤j≤i−1aiaj2 + 3∑ 0≤j≤i−1ai2a j + ai3 = ∑ 0≤j≤i−1aj 3 + 6∑ 0≤j≤i−1 ∑ 0≤k≤jaiajak − 3∑ 0≤j≤iaiaj2 + 3∑ 0≤j≤iai2a j + ai3 = ∑ 0≤j≤i−1aj 3 + 6∑ 0≤j≤i ∑ 0≤k≤jaiajak − 6∑ 0≤j≤iai2a j − 3∑ 0≤j≤iaiaj2 + 3∑ 0≤j≤iai2a j + ai3 = ∑ 0≤j≤i−1aj 3 + 6∑ 0≤j≤i ∑ 0≤k≤jaiajak − 3∑ 0≤j≤iaiaj2 − 3∑ 0≤j≤iai2a j + ai3 = ∑ 0≤j≤i−1aj 3 + 6∑ 0≤j≤i ∑ 0≤k≤jaiajak − 3∑ 0≤j≤iaiaj(ai + aj) + ai3

and in the trivial case for i = 0 that

∑ 0≤j≤iaj 3 = a 03 = (0)3 + 6a 03 − 3(2)a 03 + a 03

so that

6∑ 0≤j≤i ∑ 0≤k≤jaiajak = ∑ 0≤j≤iaj 3 −∑ 0≤j≤i−1aj 3 + 3∑ 0≤j≤iaiaj(ai + aj) − ai3

and so that for 0 ≤ i ≤ n we have 6∑ ⁡ 0≤i≤n ∑ ⁡ 0≤j≤i ∑ ⁡ 0≤k≤jaiajak equivalent to

∑ 0≤i≤n ∑ 0≤j≤iaj 3 −∑ 0≤j≤i−1aj 3 + 3∑ 0≤i≤n ∑ 0≤j≤iaiaj(ai + aj) −∑ 0≤i≤nai3.

We may also prove by induction that

∑ 0≤i≤n ∑ 0≤j≤iaj 3 −∑ 0≤j≤i−1aj 3 = ∑ 0≤i≤nai 3.

If n = 0, clearly a03 − (0)3 = a03. Then, assuming

∑ 0≤i≤k ∑ 0≤j≤iaj 3 −∑ 0≤j≤i−1aj 3 = ∑ 0≤i≤kai 3

we must show that

∑ 0≤i≤k+1 ∑ 0≤j≤iaj 3 −∑ 0≤j≤i−1aj 3 = ∑ 0≤i≤k+1ai 3.

But

∑ 0≤i≤k+1 ∑ 0≤j≤iaj 3 −∑ 0≤j≤i−1aj 3 = ∑ 0≤i≤k ∑ 0≤j≤iaj 3 −∑ 0≤j≤i−1aj 3 + ∑ 0≤j≤k+1aj 3 −∑ 0≤j≤kaj 3 = ∑ 0≤i≤kai 3 + ∑ 0≤j≤k+1aj 3 −∑ 0≤j≤kaj 3 = ∑ 0≤j≤k+1aj 3

And so finally, we have that

6∑ 0≤i≤n ∑ 0≤j≤i ∑ 0≤k≤jaiajak = ∑ 0≤i≤nai 3 + 3∑ 0≤i≤n ∑ 0≤j≤iaiaj(ai + aj) −∑ 0≤i≤nai3

and

6∑ 0≤i≤n ∑ 0≤j≤i ∑ 0≤k≤jaiajak = ∑ 0≤i≤nai 3 + 3∑ 0≤i≤n ∑ 0≤j≤iaiaj(ai + aj) −∑ 0≤i≤nai3 = ∑ 0≤i≤nai 3 + ∑ 0≤i≤n ∑ 0≤j≤iaiaj(ai + aj) −∑ 0≤i≤nai3 + 2∑ 0≤i≤n ∑ 0≤j≤iaiaj(ai + aj) = ∑ 0≤i≤nai 3 + ∑ 0≤i≤n 1 2∑ 0≤j≤iaiaj(ai + aj) + 1 2∑ i≤j≤naiaj(ai + aj) −∑ 0≤i≤nai3 + 2∑ 0≤i≤n ∑ 0≤j≤iaiaj(ai + aj) = ∑ 0≤i≤nai 3 + ∑ 0≤i≤n 1 2∑ 0≤j≤naiaj2 + 1 2∑ 0≤j≤nai2a j + ai3 −∑ 0≤i≤nai3 + 2∑ 0≤i≤n ∑ 0≤j≤iaiaj(ai + aj) = ∑ 0≤i≤nai 3 + ∑ 0≤i≤nai ∑ 0≤i≤nai2 + ∑ 0≤i≤nai3 −∑ 0≤i≤nai3 + 2∑ 0≤i≤n ∑ 0≤j≤iaiaj(ai + aj) = ∑ 0≤i≤nai 3 + ∑ 0≤i≤nai ∑ 0≤i≤nai2 + 2∑ 0≤i≤n ∑ 0≤j≤iaiaj(ai + aj) = ∑ 0≤i≤nai 3 + ∑ 0≤i≤nai ∑ 0≤i≤nai2 + ∑ 0≤i≤n ∑ 0≤j≤iaiaj(ai + aj) + ∑ i≤j≤naiaj(ai + aj) = ∑ 0≤i≤nai 3 + ∑ 0≤i≤nai ∑ 0≤i≤nai2 + ∑ 0≤i≤n ∑ 0≤j≤naiaj(ai + aj) + 2ai3 = ∑ 0≤i≤nai 3 + ∑ 0≤i≤nai ∑ 0≤i≤nai2 + ∑ 0≤i≤n ∑ 0≤j≤naiaj2 + ∑ 0≤i≤n ∑ 0≤j≤nai2a j + 2∑ 0≤i≤nai3 = ∑ 0≤i≤nai 3 + ∑ 0≤i≤nai ∑ 0≤i≤nai2 + 2 ∑ 0≤i≤nai ∑ 0≤i≤nai2 + 2∑ 0≤i≤nai3 = ∑ 0≤i≤nai 3 + 3 ∑ 0≤i≤nai ∑ 0≤i≤nai2 + 2∑ 0≤i≤nai3

Therefore,

∑ 0≤i≤n ∑ 0≤j≤i ∑ 0≤k≤jaiajak = 1 6 ∑ 0≤i≤nai 3 + 1 2 ∑ 0≤i≤nai ∑ 0≤i≤nai2 + 1 3∑ 0≤i≤nai3.

▸ 30. [M23] (J. Binet, 1812.) Without using induction, prove the identity

(∑ j=1na jxj)(∑ j=1nb jyj) =(∑ j=1na jyj)(∑ j=1nb jxj)+∑ 1≤j≤k≤n(ajbk−akbj)(xjyk−xkyj).

[An important special case arises when w1,… ⁡,wn,z1,… ⁡,zn are aribtrary complex numbers and we set aj = wj, bj = z̄j, xj = w̄j, yj = zj:

(∑ j=1n|w j|2)(∑ j=1n|z j|2) =|∑ j=1nw jzj|2 + ∑ 1≤j≤k≤n|wjz̄k − wkz̄j|2.

The terms |wjz̄j|2 are nonnegative, so the famous Cauchy-Schwarz inequality

(∑ j=1n|w j|2)(∑ j=1n|z j|2) ≥|∑ j=1nw jzj|2

is a consequence of Binet’s formula.]

Proposition. ∑ ⁡ 1≤j≤najxj ∑ ⁡ 1≤j≤nbjyj = ∑ ⁡ 1≤j≤najyj ∑ ⁡ 1≤j≤nbjxj + ∑ ⁡ 1≤j≤n ∑ ⁡ j<k≤n(ajbk − akbj)(xjyk − xkyj).

Proof. We need to show that

∑ 1≤j≤najxj ∑ 1≤j≤nbjyj = ∑ 1≤j≤najyj ∑ 1≤j≤nbjxj+∑ 1≤j≤n ∑ j<k≤n(ajbk−akbj)(xjyk−xkyj).

But

∑ 1≤j≤n ∑ j<k≤n(ajbk − akbj)(xjyk − xkyj) = ∑ 1≤j≤n ∑ j<k≤najbk(xjyk − xkyj) −∑ 1≤j≤n ∑ j<k≤nakbj(xjyk − xkyj) = ∑ 1≤j≤n ∑ j<k≤najbk(xjyk − xkyj) + ∑ 1≤j≤n ∑ j<k≤nakbj(xkyj − xjyk) = ∑ 1≤j≤n ∑ j<k≤najbk(xjyk − xkyj) + ∑ 1≤k≤n ∑ k<j≤najbk(xjyk − xkyj) = ∑ 1≤j≤n ∑ j<k≤najbk(xjyk − xkyj) + ∑ 1≤j≤n ∑ 1≤k<jajbk(xjyk − xkyj) = ∑ 1≤j≤n ∑ 1≤k<jajbk(xjyk − xkyj) + 0 + ∑ 1≤j≤n ∑ j<k≤najbk(xjyk − xkyj) = ∑ 1≤j≤n ∑ 1≤k<jajbk(xjyk − yjxk) + ∑ 1≤j≤najbj(xjyj − yjxj) + ∑ 1≤j≤n ∑ j<k≤najbk(xjyk − yjxk) = ∑ 1≤j≤n ∑ 1≤k≤najbk(xjyk − yjxk) = ∑ 1≤j≤n ∑ 1≤k≤najxjbkyk −∑ 1≤j≤n ∑ 1≤k≤najyjbkxk = ∑ 1≤j≤najxj ∑ 1≤j≤nbjyj −∑ 1≤j≤najyj ∑ 1≤j≤nbjxj

as we needed to show. □

31. [M20] Use Binet’s formula to express the sum ∑ ⁡ 1≤j≤k≤n(uj − uk)(vj − vk) in terms of ∑ ⁡ j=1nujvj, ∑ ⁡ j=1nuj, and ∑ ⁡ j=1nvj.

We want to find an expression for

∑ 1≤jn ∑ j<k≤n(uj − uk)(vj − vk)

in terms of ∑ ⁡ 1≤j≤nujvj, ∑ ⁡ 1≤j≤nuj, and ∑ ⁡ 1≤j≤nvj.

From Binet’s formula we have that

∑ 1≤j≤n ∑ j<k≤n(ajbk−akbj)(xjyk−xkyj) = ∑ 1≤j≤najxj ∑ 1≤j≤nbjyj−∑ 1≤j≤najyj ∑ 1≤j≤nbjxj .

If we let aj = uj, xj = vj, and bj = yj = 1, we find

∑ 1≤j≤n ∑ j<k≤n(uj − uk)(vj − vk) = ∑ 1≤j≤nujvj ∑ 1≤j≤n1 −∑ 1≤j≤nuj ∑ 1≤j≤nvj = n∑ 1≤j≤nujvj −∑ 1≤j≤nuj ∑ 1≤j≤nvj .

______________________________________________________________________________________________________________________________

[See Soobschch. Mat. Obschch. Khar’kovskom Univ. 4, 2 (1882), 93–98.]

32. [M20] Prove that

∏ j=1n ∑ i=1ma ij = ∑ 1≤i1,…,in≤mai11…ainn.

Proposition. ∏ ⁡ 1≤j≤n ∑ ⁡ 1≤ij≤maijj = ∑ ⁡ 1≤i1,… ⁡,in≤mai11… ⁡ainn.

Proof. We need to show that

∏ 1≤j≤n ∑ 1≤ij≤maijj = ∑ 1≤i1,…,in≤mai11…ainn.

If n = 1, clearly ∑ ⁡ 1≤i1≤mai11 = ∑ ⁡ 1≤i1≤mai11. Then, assuming that

∏ 1≤j≤k ∑ 1≤ij≤maijj = ∑ 1≤i1,…,ik≤mai11…aikk

we must show that

∏ 1≤j≤k+1 ∑ 1≤ij≤maijj = ∑ 1≤i1,…,ik+1≤mai11…aik+1(k+1).

But

∏ 1≤j≤k+1 ∑ 1≤ij≤maijj = ∏ 1≤j≤k ∑ 1≤ij≤maijj ∑ 1≤ik+1≤maik+1(k+1) = ∑ 1≤i1,…,ik≤mai11…aikk ∑ 1≤ik+1≤maik+1(k+1) = ∑ 1≤i1,…,ik+1≤mai11…aik+1(k+1)

as we needed to show. □

▸ 33. [M30] One evening Dr. Matrix discovered some formulas that might even be classed as more remarkable than those of exercise 20:

1 (a − b)(a − c) + 1 (b − a)(b − c) + 1 (c − a)(c − b) = 0,
a (a − b)(a − c) + b (b − a)(b − c) + c (c − a)(c − b) = 0,
a2 (a − b)(a − c) + b2 (b − a)(b − c) + c2 (c − a)(c − b) = 1,
a3 (a − b)(a − c) + b3 (b − a)(b − c) + c3 (c − a)(c − b) = a + b + c.

Prove that these formulas are a special case of a general law; let x1,x2,… ⁡,xn be distinct numbers, and show that

∑ j=1n(x jr ∏ 1≤k≤n k≠j (xj−xk)) = 0, if 0 ≤ r < n − 1, 1, if r = n − 1; ∑ j=1nxj,if r = n.

Proposition. ∑ ⁡ 1≤i≤n xir ∏ ⁡ 1≤j≤n j≠i (xi−xj) = 0 if 0 ≤ r < n − 1 1 if r = n − 1 ∑ ⁡ 1≤i≤nxiif r = n if x1,x2,… ⁡,xn distinct.

Proof. For an arbitrary series of distinct numbers xi, 1 ≤ i ≤ n, and for an arbitrary ι, 1 ≤ ι ≤ n, let P(xι) = xιr for 0 ≤ r ≤ n, and Q(xι) = ∏ ⁡ 1≤i≤n i≠ι (xι−xi). By the fundamental theorem of algebra and the method of partial fractions, since r = deg ⁡ P ≤ deg ⁡ Q + 1 = n, we have that

P(xι) Q(xι) = xιr ∏ 1≤i≤n i≠ι (xι − xi) = D(xι)+∑ 1≤i≤n i≠ι ci xι − xi

for constants ci, where D(xι) is the polynomial divisor with remainder R(xι) such that

P(xι) = D(xι)Q(xι) + R(xι)deg ⁡ R < deg ⁡ Q

By polynomial division, we have

D(xι) = 0 if 0 ≤ r = deg ⁡ P < deg ⁡ Q = n − 1 1 if r = deg ⁡ P = deg ⁡ Q = n − 1 ∑ 1≤i≤nxiif r = deg ⁡ P = deg ⁡ Q + 1 = n

Also, for an arbitrary κ, 1 ≤ κ ≤ n, we have

xιr ∏ 1≤i≤n i≠ι (xι − xi) = D(xι) + ∑ 1≤i≤n i≠ι ci xι − xi ⇔ xιr (xι − xκ)∏ 1≤i≤n i≠ι i≠κ (xι − xi) = D(xι) + ∑ 1≤i≤n i≠ι i≠κ ci xι − xi + cκ xι − xκ ⇔ xιr(xι − xκ) (xι − xκ)∏ 1≤i≤n i≠ι i≠κ (xι − xi) = (xι − xκ)D(xι) + (xι − xκ)∑ 1≤i≤n i≠ι i≠κ ci xι − xi + cκ(xι − xκ) xι − xκ ⇔ xιr ∏ 1≤i≤n i≠ι i≠κ (xι − xi) = (xι − xκ)D(xι) + (xι − xκ)∑ 1≤i≤n i≠ι i≠κ ci xι − xi + cκ ⇔cκ = xιr ∏ 1≤i≤n i≠ι i≠κ (xι − xi) + (xκ − xι)D(xι) + (xκ − xι)∑ 1≤i≤n i≠ι i≠κ ci xι − xi

Letting ι = κ, we find

cκ = xκr ∏ 1≤i≤n i≠κ (xκ − xi) + (xκ − xκ)D(xκ) + (xκ − xκ)∑ 1≤i≤n i≠κ ci xκ − xi = xκr ∏ 1≤i≤n i≠κ (xκ − xi)

And so

xιr ∏ 1≤i≤n i≠ι (xι − xi) = D(xι) + ∑ 1≤i≤n i≠ι ci xι − xi = D(xι) + ∑ 1≤i≤n i≠ι xir (xι − xi)∏ 1≤j≤n j≠i (xi − xj) = D(xι) −∑ 1≤i≤n i≠ι xir (xi − xι)∏ 1≤j≤n j≠i (xi − xj)

or equivalently

D(xι) = xιr ∏ 1≤i≤n i≠ι (xι − xi)+∑ 1≤i≤n i≠ι xir (xi − xι)∏ 1≤j≤n j≠i (xi − xj)

letting ι = n yields

D(xn) = xnr ∏ 1≤i≤n i≠n (xn − xi) + ∑ 1≤i≤n i≠n xir (xi − xn)∏ 1≤j≤n j≠i (xi − xj) = ∑ 1≤i≤n−1 xir (xi − xn)∏ 1≤j≤n−1 j≠i (xi − xj) + xnr ∏ 1≤i≤n−1 (xn − xi) = ∑ 1≤i≤n−1 ∏ 1≤j≤n−1 j≠i ∏ 1≤k≤n−1 k≠j (xj − xk) xir xi−xn ∏ 1≤i≤n−1 ∏ 1≤j≤n−1 j≠i (xi − xj) + xnr ∏ 1≤i≤n−1 (xn − xi) = ∑ 1≤i≤n−1 1 xi−xn ∏ 1≤j≤n−1 j≠i ∏ 1≤k≤n−1 k≠j (xj − xk)xir ∏ 1≤i≤n−1 ∏ 1≤j≤n−1 j≠i (xi − xj) + xnr ∏ 1≤i≤n−1 (xn − xi) = ∑ 1≤i≤n−1 ∏ 1≤j≤n−1 j≠i (xj − xn) ∏ 1≤j≤n−1 j≠i ∏ 1≤k≤n−1 k≠j (xj − xk)xir ∏ 1≤i≤n−1(xi − xn) ∏ 1≤i≤n−1 ∏ 1≤j≤n−1 j≠i (xi − xj) + xnr ∏ 1≤i≤n−1 (xn − xi) = ∑ 1≤i≤n−1 ∏ 1≤j≤n−1 (xn − xj) ∏ 1≤j≤n−1 j≠i (xj − xn) ∏ 1≤j≤n−1 j≠i ∏ 1≤k≤n−1 k≠j (xj − xk)xir ∏ 1≤j≤n−1 (xn − xj) ∏ 1≤i≤n−1(xi − xn) ∏ 1≤i≤n−1 ∏ 1≤j≤n−1 j≠i (xi − xj) + ∏ 1≤i≤n−1 (xi − xn) ∏ 1≤i≤n−1 ∏ 1≤j≤n−1 j≠i (xi − xj)xnr ∏ 1≤j≤n−1 (xn − xj) ∏ 1≤i≤n−1(xi − xn) ∏ 1≤i≤n−1 ∏ 1≤j≤n−1 j≠i (xi − xj) = U V

where

U = ∑ 1≤i≤n−1 ∏ 1≤j≤n−1 (xn − xj) ∏ 1≤j≤n−1 j≠i (xj − xn) ∏ 1≤j≤n−1 j≠i ∏ 1≤k≤n−1 k≠j (xj − xk)xir + ∏ 1≤i≤n−1 (xi − xn) ∏ 1≤i≤n−1 ∏ 1≤j≤n−1 j≠i (xi − xj)xnr = ∑ 1≤i≤n−1 ∏ 1≤j≤n−1 (xn − xj) ∏ 1≤j≤n−1 j≠i (xj − xn) ∏ 1≤j≤n−1 j≠i ∏ 1≤k≤n−1 k≠j (xj − xk)xir + ∏ 1≤i≤n−1 ∏ 1≤j≤n−1 j≠i (xi − xj)(xi − xn)xnr = ∑ 1≤i≤n−1 ∏ 1≤j≤n−1 (xn − xj) ∏ 1≤j≤n−1 j≠i (xj − xn) ∏ 1≤j≤n−1 j≠i ∏ 1≤k≤n−1 k≠j (xj − xk)xir + ∏ 1≤i≤n−1 ∏ 1≤j≤n j≠i (xi − xj)xnr = ∑ 1≤i≤n−1 ∏ 1≤j≤n−1 j≠i ∏ 1≤k≤n−1 k≠j (xj − xk)(xj − xn) ∏ 1≤k≤n k≠n (xn − xk)xir + ∏ 1≤j≤n j≠n ∏ 1≤k≤n k≠j (xj − xk)xnr = ∑ 1≤i≤n−1 ∏ 1≤j≤n−1 j≠i ∏ 1≤k≤n k≠j (xj − xk) ∏ 1≤k≤n k≠n (xn − xk)xir + ∏ 1≤j≤n j≠n ∏ 1≤k≤n k≠j (xj − xk)xnr = ∑ 1≤i≤n−1 ∏ 1≤j≤n j≠i ∏ 1≤k≤n k≠j (xj − xk)xir + ∏ 1≤j≤n j≠n ∏ 1≤k≤n k≠j (xj − xk)xnr = ∑ 1≤i≤n ∏ 1≤j≤n j≠i ∏ 1≤k≤n k≠j (xj − xk)xir

and where

V = ∏ 1≤j≤n−1 (xn − xj) ∏ 1≤i≤n−1(xi − xn) ∏ 1≤i≤n−1 ∏ 1≤j≤n−1 j≠i (xi − xj) = ∏ 1≤i≤n−1 ∏ 1≤j≤n−1 j≠i (xi − xj)(xi − xn) ∏ 1≤j≤n−1 (xn − xj) = ∏ 1≤i≤n−1 ∏ 1≤j≤n j≠i (xi − xj) ∏ 1≤j≤n−1 (xn − xj) = ∏ 1≤i≤n ∏ 1≤j≤n j≠i (xi − xj)

so that

D(xn) = U V = ∑ 1≤i≤n ∏ 1≤j≤n j≠i ∏ 1≤k≤n k≠j (xj − xk)xir ∏ 1≤i≤n ∏ 1≤j≤n j≠i (xi − xj) = ∑ 1≤i≤n xir ∏ 1≤j≤n j≠i (xi − xj).

That is,

∑ 1≤i≤n xir ∏ 1≤j≤n j≠i (xi − xj) = 0 if 0 ≤ r < n − 1 1 if r = n − 1 ∑ 1≤i≤nxiif r = n

as we needed to show. □

______________________________________________________________________________________________________________________________

Notes: Dr. Matrix was anticipated in this discovery by L. Euler, who wrote to Christian Golbach about it on 9 November 1762. See Euler’s Institutionum Calculi Integralis 2 (1769), §1169; and E. Waring, Phil. Trans. 69 (1779), 64–67… ⁡ [J. J. Sylvester, Quart. J. Math. 1 (1857), 141–152.]

34. [M25] Prove that

∑ k=1n ∏ 1≤r≤n, r≠m(x + k − r) ∏ 1≤r≤n, r≠k(k − r) = 1,

provided that 1 ≤ m ≤ n and x is arbitrary. For example, if n = 4 and m = 2, then

x(x − 2)(x − 3) (−1)(−2)(−3) + (x + 1)(x − 1)(x − 2) (1)(−1)(−2) + (x + 2)x(x − 1) (2)(1)(−1) + (x + 3)(x + 1)x (3)(2)(1) = 1.

Proposition. ∑ ⁡ 1≤i≤n∏ ⁡ 1≤j≤n, j≠k(x+i−j) ∏ ⁡ 1≤j≤n, j≠i(i−j) = 1 provided that 1 ≤ k ≤ n and x is arbitrary.

Proof. We may prove

∑ 1≤i≤n ∏ 1≤j≤n, j≠k(x + i − j) ∏ 1≤j≤n, j≠i(i − j) = 1

provided that 1 ≤ k ≤ n and x is arbitrary, but first we shall first prove the more general result

∑ 1≤i≤n ∏ 1≤j≤n,j≠k(yi − zj) ∏ 1≤j≤n,j≠i(yi − yj) = 1.

Let P(y) be the polynomial representation of ∏ ⁡ 1≤j≤n,j≠k(y − zj) where

P(y) = ∏ 1≤j≤n j≠k (y−zj) = ∑ 0≤j≤n−1cjyj

for arbitrary coefficients c0,c1,… ⁡,cn−1. Note that since P(y) = yn−1 + … ⁡, cn−1 = 1. Then, from exercise 33, we find that

∑ 1≤i≤n ∏ 1≤j≤n,j≠k(yi − zj) ∏ 1≤j≤n,j≠i(yi − yj) = ∑ 1≤i≤n P(yi) ∏ 1≤j≤n,j≠i(yi − yj) = ∑ 1≤i≤n ∑ 0≤j≤n−1cjyij ∏ 1≤j≤n,j≠i(yi − yj) = ∑ 1≤j≤n ∑ 0≤i≤n−1ciyji ∏ 1≤k≤n,k≠j(yj − yk) = ∑ 0≤i≤n−1ci ∑ 1≤j≤n yji ∏ 1≤k≤n,k≠j(yj − yk) = ∑ 0≤i<n−1ci ∑ 1≤j≤n yji ∏ 1≤k≤n,k≠j(yj − yk) + cn−1 ∑ 1≤j≤n yjn−1 ∏ 1≤k≤n,k≠j(yj − yk) = ∑ 0≤i<n−1ci 0 + cn−1 1 = cn−1 = 1.

Letting yi = i and zi = i + x for 1 ≤ i ≤ n, 1 ≤ k ≤ n, and x arbitrary, we have that

∑ 1≤i≤n ∏ 1≤j≤n, j≠k(x + i − j) ∏ 1≤j≤n, j≠i(i − j) = 1

as we needed to show. □

35. [HM20] The notation sup ⁡ R(j)aj is used to denote the least upper bound of the elements aj, in a manner analogous to the ∑ ⁡ - and ∏ ⁡ -notations. (When R(j) is satisfied for only finitely many j, the notation max ⁡ R(j)aj is often used to denote the same quantity.) Show how rules (a), (b), (c), and (d) can be adapted for manipulation of this notation. In particular discuss the following analog of rule (a):

(sup ⁡ R(j)ai) + (sup ⁡ S(j)bj) = sup ⁡ R(i)(sup ⁡ S(j)(ai + bj)),

and give a suitable definition for the notation when R(j) is satisfied for no j.

We have the following analogs for the least upper bound: an additive law:

(sup ⁡ R(i)ai) + (sup ⁡ S(j)bj) = sup ⁡ R(i)(sup ⁡ S(j)(ai + bj))

as well as a multiplicative law:

(sup ⁡ R(i)ai)(sup ⁡ S(j)bj) = sup ⁡ R(i)(sup ⁡ S(j)(aibj))

provided ai and aj are nonnegative; change of variable:

sup ⁡ R(i)ai = sup ⁡ R(j) = sup ⁡ R(p(j))ap(j);

interchanging order of bound:

sup ⁡ R(i) sup ⁡ S(j)aij = sup ⁡ S(j) sup ⁡ R(i)aij;

and manipulating the domain:

sup ⁡ (sup ⁡ R(j)aj,sup ⁡ S(j)aj) = sup ⁡ R(j)∨S(j)aj.

A suitable definition for the notation when R(j) is satisfied for no j would be

sup ⁡ j∈∅ = −∞

since −∞ acts as the identity element for the least upper bound, in that sup ⁡ (sup ⁡ j∈∅aj,sup ⁡ R(j)aj) = sup ⁡ R(j)aj.

36. [M23] Show that the determinant of the combinatorial matrix is xn−1(x + ny).

Proposition. det ⁡ [y + δijx]n = xn−1(ny + x).

Proof. For

An = [aij]n = [y+δijx]n = y + x y ⋯ y y y + x⋯ y ⋮ ⁡ ⋮ ⁡ ⋱ ⋮ ⁡ y y ⋯y + x n,

we must show that

det ⁡ [aij]n = xn−1(ny + x).

Let

aij′ = ai1 if j = 1 aij − ai1if 2 ≤ j ≤ n

and

aij″ ⁡ = a1j′ + ∑ 2≤k≤nakj′if i = 1 aij′ if 2 ≤ i ≤ n

so that det ⁡ [aij]n = det ⁡ [aij″ ⁡ ]n where

aij″ ⁡ = ny + xif i = 1, j = 1 0 if i = 1, 2 ≤ j ≤ n y if 2 ≤ i ≤ n, j = 1 δijx if 2 ≤ i ≤ n, 2 ≤ j ≤ n

since:

a11′ = a 11 = y + x i = 1,j = 1 a1j′ = −x i = 1,2 ≤ j ≤ n ai1′ = a i1 = y 2 ≤ i ≤ n,j = 1 aij′ = δ ijx 2 ≤ i ≤ n,2 ≤ j ≤ n a11″ ⁡ = a 11′ + ∑ 2≤k≤nak1′ = y + x + ∑ 2≤k≤ny = ny + x i = 1,j = 1 a1j″ ⁡ = a 1j′ + ∑ 2≤k≤nakj′ = −x + ∑ 2≤k≤nδkjx = 0 i = 1,2 ≤ j ≤ n ai1″ ⁡ = a i1′ = y 2 ≤ i ≤ n,j = 1 aij″ ⁡ = a ij′ = δ ijx 2 ≤ i ≤ n,2 ≤ j ≤ n

Then, since [aij″ ⁡]n is a triangular matrix (aij″ ⁡ = 0 whenever i < j), we have that

det ⁡ [aij]n = det ⁡ [aij″ ⁡ ] n = ∏ 1≤k≤nakk″ = a11″ ⁡ ∏ 2≤k≤nakk″ = (ny + x)∏ 2≤k≤nδkkx = (ny + x)∏ 2≤k≤nx = (ny + x)xn−1 = xn−1(ny + x)

as we needed to show. □

▸ 37. [M24] Show that the determinant of Vandermonde’s matrix is

∏ 1≤j≤nxj ∏ 1≤i<j≤n(xj − xi).

Proposition. det ⁡ [xji]n = ∏ ⁡ 1≤i≤nxi ∏ ⁡ 1≤j<i≤n(xi − xj).

Proof. For

An = [aij]n = [xji] n = x1 x2 ⋯ xn x12x22⋯xn2 ⋮ ⁡ ⋮ ⁡ ⋱ ⋮ ⁡ x 1nx 2n⋯x nn n,

we must show that

det ⁡ [aij]n = ∏ 1≤i≤nxi ∏ 1≤j<i≤n(xi − xj).

Let

aij′ = ai1 if j = 1 aij − ai1if 2 ≤ j ≤ n

and

aij″ ⁡ = a11′ if i = 1,j = 1 a1j′ if i = 1,2 ≤ j ≤ n ai1′− x1a(i−1)1′if 2 ≤ i ≤ n,j = 1, aij′− x1a(i−1)j′if 2 ≤ i ≤ n,2 ≤ j ≤ n,

so that det ⁡ [aij] = det ⁡ [aij″ ⁡ ] where

aij″ ⁡ = x1 if i = 1, j = 1 xj − x1 if i = 1, 2 ≤ j ≤ n 0 if 2 ≤ i ≤ n, j = 1 xji−1(xj − x1)if 2 ≤ i ≤ n, 2 ≤ j ≤ n

since:

a11′ = a 11 = x1 i = 1,j = 1 a1j′ = a 1j − a11 = xj − x1 i = 1,2 ≤ j ≤ n ai1′ = a i1 = x1i 2 ≤ i ≤ n,j = 1 aij′ = a ij − ai1 = xji − x 1i 2 ≤ i ≤ n,2 ≤ j ≤ n a11″ ⁡ = a 11′ = a 11 = x1 i = 1,j = 1 a1j″ ⁡ = a 1j′ = a 1j − a11 = xj − x1 i = 1,2 ≤ j ≤ n ai1″ ⁡ = a i1′− x 1a(i−1)1′ = a i1 − x1a(i−1)1) = 0 2 ≤ i ≤ n,j = 1 aij″ ⁡ = a ij′− x 1a(i−1)j′ = a ij − ai1 − x1(a(i−1)j − a(i−1)1) = xji−1(x j − x1)2 ≤ i ≤ n,2 ≤ j ≤ n

Let

qij = xj+1 − x1if i = j,1 ≤ i,j ≤ n − 1 0 otherwise

so that minor ⁡ ([a″ ⁡ ]n,1,1) = ([qij]n−1 minor ⁡ ([a]n,1,1)T)T and det ⁡ [qij]n−1 = ∏ ⁡ 1≤k≤n−1(xk+1 − x1). Then

det ⁡ [aij]n = det ⁡ [aij″ ⁡ ] n = ∑ 1≤i≤nai1″ cofactor ⁡ (a i1″ ⁡ ) = a11″ ⁡ cofactor ⁡ (a 11″ ⁡ ) + ∑ 2≤i≤nai1″ cofactor ⁡ (a i1″ ⁡ ) = x1(−1)1+1 det ⁡ minor ⁡ ([a″ ⁡ ] n,1,1) + 0 = x1 det ⁡ minor ⁡ ([a″ ⁡ ]n,1,1) = x1 det ⁡ (([qij]n−1 minor ⁡ ([a]n,1,1)T)T) = x1 det ⁡ ([qij]n−1 minor ⁡ ([a]n,1,1)T) = x1 det ⁡ [qij]n−1 det ⁡ (minor ⁡ ([a]n,1,1)T) = x1 det ⁡ [qij]n−1 det ⁡ minor ⁡ ([a]n,1,1) = x1 ∏ 1≤k≤n−1(xk+1 − x1)det ⁡ minor ⁡ ([a]n,1,1).

We shall finally use this recursive identity to give a proof by mathematical induction on n.

If n = 1, then clearly

det ⁡ [aij]1 = det ⁡ [a11]1 = x1 = ∏ 1≤i≤nxi ∏ 1≤j<i≤1(xi − xj).

Then, assuming that

det ⁡ [aij]k = ∏ 1≤i≤kxi ∏ 1≤j<i≤k(xi − xj)

or equivalently that

det ⁡ minor ⁡ ([a]k+1,1,1) = ∏ 2≤i≤k+1xi ∏ 2≤j<i≤k+1(xi − xj)

we must show that

det ⁡ [aij]k+1 = ∏ 1≤i≤k+1xi ∏ 1≤j<i≤k+1(xi − xj).

But

det ⁡ [aij]k+1 = x1 ∏ 1≤i≤k(xi+1 − x1)det ⁡ [a(i+1)(j+1)]k = x1 ∏ 1≤i≤k(xi+1 − x1) ∏ 2≤i≤k+1xi ∏ 2≤j<i≤k+1(xi − xj) = x1 ∏ 2≤i≤k+1(xi − x1) ∏ 2≤i≤k+1xi ∏ 2≤j<i≤k+1(xi − xj) = ∏ 1≤i≤1xi ∏ 1≤j<i≤k+1(xi − xj) ∏ 2≤i≤k+1xi ∏ 2≤j<i≤k+1(xi − xj) = ∏ 1≤i≤k+1xi ∏ 1≤j<i≤k+1(xi − xj)

as we needed to show. □

▸ 38. [M25] Show that the determinant of Cauchy’s matrix is

∏ 1≤i<j≤n(xj − xi)(yj − yi)∏ 1≤i,j≤n(xi + yi).

Proposition. det ⁡ [1∕(xi + yj)]n = ∏ ⁡ 1≤i<j≤n(xj − xi)(yj − yi)∏ ⁡ 1≤i,j≤n(xi + yj).

Proof. For

An = [aij]n = [1∕(xi+yi)]n = 1∕(x1 + y1)1∕(x1 + y2)⋯1∕(x1 + yn) 1∕(x2 + y1)1∕(x2 + y2)⋯1∕(x2 + yn) ⋮ ⁡ ⋮ ⁡ ⋱ ⋮ ⁡ 1∕(xn + y1)1∕(xn + y2)⋯1∕(xn + yn) n,

we must show that

det ⁡ [aij]n = ∏ 1≤i<j≤n(xj − xi)(yj − yi)∏ 1≤i,j≤n(xi + yj).

Let

aij′ = ai1 if j = 1 aij − ai1if 2 ≤ j ≤ n

so that det ⁡ [aij] = det ⁡ [aij′] where

aij′ = 1∕(x1 + y1) if i = 1, j = 1 ((y1 − yj)∕(x1 + y1))(1∕(x1 + yj))if i = 1, 2 ≤ j ≤ n 1∕(x1 + y1) if 2 ≤ i ≤ n, j = 1 ((y1 − yj)∕(xi + y1))(1∕(xi + yj)) if 2 ≤ i ≤ n, 2 ≤ j ≤ n

since:

a11′ = a 11 = 1∕(x1 + y1) i = 1,j = 1 a1j′ = a 1j − a11 = ((y1 − yj)∕(x1 + y1))(1∕(x1 + yj)) i = 1,2 ≤ j ≤ n ai1′ = a i1 = 1∕(xi + y1) 2 ≤ i ≤ n,j = 1 aij′ = a ij − ai1 = ((y1 − yj)∕(xi + y1))(1∕(xi + yj)) 2 ≤ i ≤ n,2 ≤ j ≤ n

Let

bij = 1 if j = 1 1∕(xi + yj)otherwise,

pij = 1∕(xi + y1)if i = j,1 ≤ i,j ≤ n 0 otherwise, and qij = 1 if i = j = 1 y1 − yjif i = j,2 ≤ i,j ≤ n 0 otherwise

so that

[aij′] n = ([qij]n([pij]n[bij]n)T)T

and:

det ⁡ [pij]n = ∏ 1≤k≤n 1 xk + y1 det ⁡ [qij]n = ∏ 2≤k≤ny1 − yk

Also let

bij′ = b1j if i = 1 bij − b1jif 2 ≤ i ≤ n

so that det ⁡ [bij] = det ⁡ [bij′] where

bij′ = 1 if i = 1, j = 1 1∕(xi + yj) if i = 1, 2 ≤ j ≤ n 0 if 2 ≤ i ≤ n, j = 1 ((x1 − xi)∕(x1 + yj))(1∕(xi + yj))if 2 ≤ i ≤ n, 2 ≤ j ≤ n

since:

b11′ = b 11 = 1 i = 1,j = 1 b1j′ = b 1j = 1∕(xi + yj) i = 1,2 ≤ j ≤ n bi1′ = b i1 − b11 = 0 2 ≤ i ≤ n,j = 1 bij′ = b ij − b1j = ((x1 − xi)∕(x1 + yj))(1∕(xi + yj)) 2 ≤ i ≤ n,2 ≤ j ≤ n

Also let

rij = 1 if i = j = 1 x1 − xiif i = j,2 ≤ i,j ≤ n 0 otherwise

and

sij = 1 if i = j = 1 1∕(x1 + yj)if i = j,2 ≤ i,j ≤ n 0 otherwise

so that

minor ⁡ ([b′] n,1,1) = ([sij]n−1([rij]n−1 minor ⁡ ([a]n,1,1))T)T

and:

det ⁡ [rij]n−1 = ∏ 2≤k≤nx1 − xk det ⁡ [sij]n−1 = ∏ 2≤k≤n 1 x1 + yj

Then

det ⁡ [aij] = det ⁡ [aij″ ⁡ ] n = det ⁡ (([qij]n([pij]n[bij]n)T)T) = det ⁡ ([qij]n)det ⁡ ([pij]n[bij]n) = det ⁡ ([qij]n)det ⁡ ([pij]n)det ⁡ ([bij]n) = ∏ 2≤k≤ny1 − yk det ⁡ ([pij]n)det ⁡ ([bij]n) = ∏ 2≤k≤ny1 − yk ∏ 1≤k≤n 1 xk + y1 det ⁡ ([bij]n) = ∏ 2≤k≤ny1 − yk ∏ 1≤k≤n 1 xk + y1 det ⁡ ([bij′] n) = ∏ 2≤k≤ny1 − yk ∏ 1≤k≤n 1 xk + y1 b11′cofactor ⁡ (b 11′) + ∑ 2≤i≤nbi1′cofactor ⁡ (b i1′) = ∏ 2≤k≤ny1 − yk ∏ 1≤k≤n 1 xk + y1 b11′cofactor ⁡ (b 11′) + 0 = ∏ 2≤k≤ny1 − yk ∏ 1≤k≤n 1 xk + y1 det ⁡ (minor ⁡ ([b′] n,1,1)) + 0 = ∏ 2≤k≤ny1 − yk ∏ 1≤k≤n 1 xk + y1 det ⁡ (([sij]n−1([rij]n−1 minor ⁡ ([a]n,1,1))T)T) = ∏ 2≤k≤ny1 − yk ∏ 1≤k≤n 1 xk + y1 det ⁡ ([sij]n−1)det ⁡ ([rij]n−1 minor ⁡ ([a]n,1,1)) = ∏ 2≤k≤ny1 − yk ∏ 1≤k≤n 1 xk + y1 det ⁡ ([sij]n−1)det ⁡ ([rij]n−1)det ⁡ (minor ⁡ ([a]n,1,1)) = ∏ 2≤k≤ny1 − yk ∏ 1≤k≤n 1 xk + y1 ∏ 2≤k≤n 1 x1 + yj det ⁡ ([rij]n−1)det ⁡ (minor ⁡ ([a]n,1,1)) = ∏ 2≤k≤ny1 − yk ∏ 1≤k≤n 1 xk + y1 ∏ 2≤k≤n 1 x1 + yj ∏ 2≤k≤nx1 − xk det ⁡ (minor ⁡ ([a]n,1,1)) = ∏ 2≤i≤n(xi − x1)(yi − y1) ∏ 1≤i,j≤n(xi + y1)(x1 + yj)det ⁡ (minor ⁡ ([a]n,1,1)).

We shall finally use this recursive identity to give a proof by mathematical induction on n.

If n = 1, then clearly

det ⁡ [aij]1 = det ⁡ [a11]1 = 1∕(x1 + y1) = ∏ 1≤i<j≤1(xj − xi)(yj − yi)∏ 1≤i,j≤1(xi + yj).

Then, assuming that

det ⁡ [aij]k = ∏ 1≤i<j≤k(xj − xi)(yj − yi)∏ 1≤i,j≤k(xi + yj)

or equivalently that

det ⁡ minor ⁡ ([a]k+1,1,1) = ∏ 2≤i<j≤k+1(xj − xi)(yj − yi)∏ 2≤i,j≤k+1(xi + yj)

we must show that

det ⁡ [aij]k+1 = ∏ 1≤i<j≤k+1(xj − xi)(yj − yi)∏ 1≤i,j≤k+1(xi + yj).

But

det ⁡ [aij]k+1 = ∏ 2≤i≤k+1(xi − x1)(yi − y1) ∏ 1≤i,j≤k+1(xi + y1)(x1 + yj)det ⁡ minor ⁡ ([a]k+1,1,1) = ∏ 2≤i≤k+1(xi − x1)(yi − y1) ∏ 1≤i,j≤k+1(xi + y1)(x1 + yj) ∏ 2≤i<j≤k+1(xj − xi)(yj − yi) ∏ 2≤i,j≤k+1(xi + yj) = ∏ 2≤j≤k+1(xj − x1)(yj − y1) ∏ 2≤i<j≤k+1(xj − xi)(yj − yi) ∏ 1≤i,j≤k+1(xi + y1)(x1 + yj) ∏ 2≤i,j≤k+1(xi + yj) = ∏ 1≤i<j≤k+1(xj − xi)(yj − yi) ∏ 2≤i<j≤k+1(xj − xi)(yj − yi) ∏ 1≤i,j≤k+1(xi + y1)(x1 + yj) ∏ 2≤i,j≤k+1(xi + yj) = ∏ 1≤i<j≤k+1(xj − xi)(yj − yi)∏ 1≤i,j≤k+1(xi + yj)

as we needed to show. □

39. [M23] Show that the inverse of a combinatorial matrix is a combinatorial matrix with the entries bij = (−y + δij(x + ny))∕x(x + ny).

Proposition. [y + δijx]n−1 = −y+δij(x+ny) x(x+ny) n.

Proof. For

An = [aij]n = [y+δijx]n = y + x y ⋯ y y y + x⋯ y ⋮ ⁡ ⋮ ⁡ ⋱ ⋮ ⁡ y y ⋯y + x n,

we must show that

[aij]n−1 = − y + δij(x + ny) x(x + ny) n.

Let In = [δij]n and Jn = [1]n so that An = yJn + xIn, and note that Jn2 = nJn. Then

An(−yJn + xIn) = x2I n − ny2J n ⇔ An(−yJn + xIn) + ny2J n = x2I n ⇔ An(−yJn + xIn) + nyInyJn = x2I n ⇔ An(−yJn + xIn) + nyInyJn + xnyIn = x2I n + xnyIn ⇔ An(−yJn + xIn) + nyIn(yJn + xIn) = x2I n + xnyIn ⇔ An(−yJn + xIn) + nyInAn = x(x + ny)In ⇔ An(−yJn + xIn + nyIn) = x(x + ny)In ⇔ An((x + ny)In − yJn) = x(x + ny)In ⇔ An (x + ny)In − yJn x(x + ny) = In.

That is, for

Bn = (x + ny)In − yJn x(x + ny) = (x + ny)[δij]n − y[1]n x(x + ny) = (x + ny)δij − y x(x + ny) n = − y + δij(x + ny) x(x + ny) n,

AnBn = In, or equivalently, that

An−1 = − y + δij(x + ny) x(x + ny) n

as we needed to show. □

40. [M24] Show that the inverse of Vandermonde’s matrix is given by

bij =(∑ 1≤k1<⋯<kn−j≤n k1,…,kn−j≠i (−1)j−1x k1…xkn−j)xi ∏ 1≤k≤n k≠i (xk−xi).

Don’t be dismayed by the complicated sum in the numerator—it is just the coefficient of xj−1 in the polynomial (x1 − x)… ⁡(xn − x)∕(xi − x).

Proposition. [xji]n−1 = ∑ ⁡ 1≤k1<⋯<kn−j≤n k1,… ⁡,kn−j≠i (−1)j−1xk1… ⁡xkn−j xi ∏ ⁡ 1≤k≤n k≠i (xk − xi).

Proof. Let

An = [aij]n = [xji] n = x1 x2 ⋯ xn x12x22⋯xn2 ⋮ ⁡ ⋮ ⁡ ⋱ ⋮ ⁡ x 1nx 2n⋯x nn n.

We must show that

[aij]n−1 = ∑ 1≤k1<…<kn−j≤n k1,…,kn−j≠i (−1)j−1x k1⋯xkn−jxi ∏ 1≤k≤n k≠i (xk − xi)n.

Let [bij]n = [aij]n−1, so that by the definition of inverse and matrix multiplication, we have

[bij]n[aij]n = ∑ 1≤k≤nbikakj n = ∑ 1≤k≤nbikxjk n = [δij]n.

We require

∑ 1≤k≤nbikxjk = δ ij,

or equivalently, given a polynomial Pi(x) = ∑ ⁡ 1≤k≤nbikxk for arbitrary x, we require

Pi(x) = δij

with given data points (δi1,x1),(δi2,x2),… ⁡,(δin,xn). By polynomial interpolation, expanding the polynomial to an initial but trivial (n + 1)th variable x0 = 0, with bi0 = b0j = 0, in order to obtain a complete set of differences, we have

Pi(x) = ∑ 0≤k≤n δik ∏ 0≤m≤n m≠k x − xm xk − xm = ∑ 0≤k≤n k=i δik ∏ 0≤m≤n m≠k x − xm xk − xm + ∑ 0≤k≤n k≠i δik ∏ 0≤m≤n m≠k x − xm xk − xm = δii ∏ 0≤m≤n m≠i x − xm xi − xm + 0 = ∏ 0≤k≤n k≠i x − xk xi − xk = x − x0 xi − x0 ∏ 1≤k≤n k≠i x − xk xi − xk = x xi ∏ 1≤k≤n k≠i xk − x xk − xi.

For x = xj, we have

∑ 1≤k≤nbikxjk = P i(xj) = xj xi ∏ 1≤k≤n k≠i xk − xj xk − xi = δij.

What is left is to find the coefficients bij. By de Moivre’s work1 , we may do so.

The real and different roots of Pi(x) are exactly those xr where r≠i, since in such a case, we have

Pi(xr) = xr xi ∏ 1≤k≤n k≠i xk − xr xk − xi = (xr−xr)xr xi ∏ 1≤k≤n k≠i k≠r xk − xr xk − xi = 0.

Let these roots be denoted by xr1,xr2,… ⁡,xrn.

Since matrix multiplication with inverse is commutative, We also have

∑ 1≤k≤nbikxjk = ∑ 1≤k≤nxkib kj = δij.

By de Moivre’s identities, with our trivially expanded polynomial,

δij = ∑ 0≤k≤nxkib kj ⇔ bkj = ∑ 1≤m≤n(−1)m ∑ 1≤r1<⋯<rm≤n r1,…,rm≠i xr1⋯xrmδ(n−m)jxk ∏ 1≤m≤n m≠k (xk − xm)

we then have, since δ(n−k)j requires j = n − k and since n − j and j − 1 have opposite parity, that

bij = ∑ 1≤k≤n(−1)k ∑ 1≤r1<⋯<rk≤n r1,…,rk≠i xr1⋯xrkδ(n−k)jxi ∏ 1≤k≤n k≠i (xi − xk) = (−1)n−j ∑ 1≤r1<⋯<rn−j≤n r1,…,rn−j≠i xr1⋯xrn−jxi ∏ 1≤k≤n k≠i (xi − xk) = (−1)(−1)j−1 ∑ 1≤r1<⋯<rn−j≤n r1,…,rn−j≠i xr1…xrn−j(−1)xi ∏ 1≤k≤n k≠i (xk − xi) = ∑ 1≤r1<⋯<rn−j≤n r1,…,rn−j≠i (−1)j−1x r1⋯xrn−jxi ∏ 1≤k≤n k≠i (xk − xi).

Therefore

[aij]n−1 = ∑ 1≤k1<…<kn−j≤n k1,…,kn−j≠i (−1)j−1x k1⋯xkn−jxi ∏ 1≤k≤n k≠i (xk − xi)n

as we needed to show. □

______________________________________________________________________________________________________________________________

[A. de Moivre, The Doctrine of Chances, 2nd edition (London: 1738), 197–199.]

41. [M26] Show that the inverse of Cauchy’s matrix is given by

bij =(∏ 1≤k≤n(xj+yk)(xk+yi))(xj+yi)(∏ 1≤k≤n k≠j (xj−xk))(∏ 1≤k≤n k≠i (yi−yk)).

Let

An = [aij]n = [1∕(xi+yi)]n = 1∕(x1 + y1)1∕(x1 + y2)⋯1∕(x1 + yn) 1∕(x2 + y1)1∕(x2 + y2)⋯1∕(x2 + yn) ⋮ ⁡ ⋮ ⁡ ⋱ ⋮ ⁡ 1∕(xn + y1)1∕(xn + y2)⋯1∕(xn + yn) n,

We must show that

[aij]n−1 = ∏ 1≤k≤n(xj + yk)(xk + yi)(xj + yi) ∏ 1≤k≤n k≠j (xj − xk) ∏ 1≤k≤n k≠i (yi − yk)n.

But, by the definition of inverse,

[aij]n−1 = [cofactor ⁡ (aij)]nT det ⁡ [aij]n = [cofactor ⁡ (aji)]n ∏ 1≤u<v≤n(xv − xu)(yv − yu)∏ 1≤u,v≤n(xu + yv) = (−1)j+i det ⁡ minor ⁡ ([a]n,j,i) ∏ 1≤u<v≤n(xv − xu)(yv − yu)∏ 1≤u,v≤n(xu + yv) = (−1)j+i ∏ 1≤u<v≤n u≠j,u≠i v≠j,v≠i (xv − xu)(yv − yu)∏ 1≤u,v≤n u≠j v≠i (xu + yv) ∏ 1≤u<v≤n(xv − xu)(yv − yu)∏ 1≤u,v≤n(xu + yv) = (−1)j+i ∏ 1≤u<v≤n u≠j v≠j (xu − xv) ∏ 1≤u<v≤n u≠i v≠i (yu − yv)∏ 1≤u,v≤n u≠j v≠i (xu + yv) ∏ 1≤u<v≤n(xu − xv) ∏ 1≤u<v≤n(yu − yv)∏ 1≤u,v≤n(xu + yv) = (−1)j+i ∏ 1≤u,v≤n(xu + yv) ∏ 1≤u,v≤n u≠j v≠i (xu + yv) ∏ 1≤u<v≤n u≠j v≠j (xu − xv) ∏ 1≤u<v≤n(xu − xv) ∏ 1≤u<v≤n u≠i v≠i (yu − yv) ∏ 1≤u<v≤n(yu − yv) = (−1)j+i ∏ 1≤u,v≤n u=j∨v=i (xu + yv) 1 ∏ 1≤u<v≤n u=j∨v=j (xu − xv) 1 ∏ 1≤u<v≤n u=i∨v=i (yu − yv) = (−1)j+i(x j + yi)∏ 1≤u≤n u≠j (xu + yi)∏ 1≤v≤n v≠i (xj + yv) 1 ∏ 1≤u≤j−1(xu − xj) ∏ j+1≤v≤n(xj − xv) 1 ∏ 1≤u≤i−1(yu − yi) ∏ i+1≤v≤n(yi − yv) = (−1)j+i(x j + yi)∏ 1≤u≤n u≠j (xu + yi)∏ 1≤v≤n v≠i (xj + yv) (−1)j−1 ∏ 1≤u≤n u≠j (xj − xu) (−1)i−1 ∏ 1≤v≤n v≠i (yi − yv) = (−1)2(j+i−1)(x j + yi) ∏ 1≤u≤n u≠j (xu + yi) ∏ 1≤u≤n u≠j (xj − xu) ∏ 1≤v≤n v≠i (xj + yv) ∏ 1≤v≤n v≠i (yi − yv) = (xj + yi)2 ∏ 1≤u≤n u≠j (xu + yi) ∏ 1≤v≤n v≠i (xj + yv) (xj + yi) ∏ 1≤u≤n u≠j (xj − xu) ∏ 1≤v≤n v≠i (yi − yv) = ∏ 1≤u≤n(xu + yi) ∏ 1≤v≤n(xj + yv) (xj + yi) ∏ 1≤u≤n u≠j (xj − xu) ∏ 1≤v≤n v≠i (yi − yv) = ∏ 1≤k≤n(xj + yk)(xk + yi) (xj + yi) ∏ 1≤k≤n k≠j (xj − xk) ∏ 1≤k≤n k≠i (yi − yk)

as we needed to show.

42. [M18] What is the sum of all n2 elements in the inverse of the combinatorial matrix?

For the inverse of the combinatorial matrix

1 y + δijxn−1 = − y + δij(x + ny) x(x + ny) n

the sum of all n2 elements is simply

∑ 1≤i≤n ∑ 1≤j≤n − y + δij(x + ny) x(x + ny) = ∑ 1≤i≤n ∑ 1≤j≤n j≠i − y − y + x + ny x(x + ny) = ∑ 1≤i≤n − y(n − 1) − y + x + ny x(x + ny) = ∑ 1≤i≤n − ny + y − y + x + ny x(x + ny) = ∑ 1≤i≤n x x(x + ny) = ∑ 1≤i≤n 1 x + ny = n x + ny.

43. [M24] What is the sum of all n2 elements in the inverse of Vandermonde’s matrix? [Hint: Use exercise 33.]

For the inverse of the Vandermonde matrix

xji n−1 = ∑ 1≤k1<…<kn−j≤n k1,…,kn−j≠i xk1…xkn−j xi ∏ 1≤k≤n k≠i (xk − xi)n = bij n,

we want to find the sum of all n2 elements ∑ ⁡ 1≤i,j≤nbij.

In the case that any xκ = 1, 1 ≤ κ ≤ n, from exercise 40, we have that

[bij]n[xji] n = ∑ 1≤k≤nbikxjk n = [δij]n

or equivalently that

∑ 1≤j≤nbij = δiκ ⇒∑ 1≤i≤n ∑ 1≤j≤nbij = 1 ⇒∑ 1≤i,j≤nbij = 1 − (1 − 1∕xκ) ⇒∑ 1≤i,j≤nbij = 1 −∏ 1≤k≤n(1 − 1∕xk).

Otherwise, in the case that xκ≠1, 1 ≤ κ ≤ n, again from exercise 40, given Pi(x) = ∑ ⁡ 1≤k≤nbikxk for x = 1 we have that

∑ 1≤i,j≤nbij = ∑ 1≤i≤nPi(x) = ∑ 1≤i≤n x xi ∏ 1≤k≤n k≠i xk − x xk − xi

or equivalently, for x = 1, that

∑ 1≤i≤nPi(1) = ∑ 1≤i≤n 1 xi ∏ 1≤k≤n k≠i xk − 1 xk − xi = ∑ 1≤i≤n xi − 1 xi − 1 1 xi ∏ 1≤k≤n k≠i xk − 1 xk − xi = ∑ 1≤i≤n (xi − 1)∏ 1≤k≤n k≠i (xk − 1) xi(xi − 1)∏ 1≤k≤n k≠i (xk − xi) = ∑ 1≤i≤n ∏ 1≤k≤n(xk − 1) xi(xi − 1)∏ 1≤k≤n k≠i (xk − xi) = ∏ 1≤k≤n(xk − 1)∑ 1≤i≤n 1 xi(xi − 1)∏ 1≤k≤n k≠i (xk − xi) = ∏ 1≤k≤n(1 − xk)∑ 1≤i≤n 1 xi(xi − 1)∏ 1≤k≤n k≠i (xi − xk) = ∏ 1≤k≤n(xk − 1) ∑ 1≤i≤n 1 xi ∏ 1≤k≤n k≠i (xi − xk) −∑ 1≤i≤n 1 (xi − 1)∏ 1≤k≤n k≠i (xi − xk).

Also, from exercise 33, we have for an extrapolated and distinct x0 = 1 and for r = 0 that

∑ 0≤i≤n 1 ∏ 0≤j≤n j≠i (xi − xj) = 1 ∏ 1≤j≤n(1 − xj)+∑ 1≤i≤n 1 (xi − 1)∏ 1≤j≤n j≠i (xi − xj) = 1

or equivalently that

∑ 1≤i≤n 1 (xi − 1)∏ 1≤j≤n j≠i (xi − xj) = 1− 1 ∏ 1≤j≤n(1 − xj).

Similarly, if x0 = 0

∑ 1≤i≤n 1 xi ∏ 1≤j≤n j≠i (xi − xj) = 1− 1 ∏ 1≤j≤n(−xj).

This yields

∑ 1≤i,j≤nbij = ∑ 1≤i≤nPi(1) = ∏ 1≤k≤n(xk − 1) ∑ 1≤i≤n 1 xi ∏ 1≤k≤n k≠i (xi − xk) −∑ 1≤i≤n 1 (xi − 1)∏ 1≤k≤n k≠i (xi − xk) = ∏ 1≤k≤n(xk − 1) 1 − 1 ∏ 1≤j≤n(−xj) −1 − 1 ∏ 1≤j≤n(1 − xj) = ∏ 1≤k≤n(xk − 1) 1 ∏ 1≤j≤n(−xj) − 1 ∏ 1≤j≤n(1 − xj) = ∏ 1≤k≤n(xk − 1) ∏ 1≤j≤n(−xj) −∏ 1≤k≤n(xk − 1) ∏ 1≤j≤n(1 − xj) = ∏ 1≤k≤n(1 − xk) ∏ 1≤j≤n(1 − xj) −∏ 1≤k≤n(xk − 1) ∏ 1≤j≤n(xj) = 1 −∏ 1≤k≤n(xk − 1) ∏ 1≤j≤n(xj) = 1 −∏ 1≤k≤nxk − 1 xk = 1 −∏ 1≤k≤n(1 − 1∕xk).

Therefore, in all cases,

∑ 1≤i,j≤nbij = 1 −∏ 1≤k≤n(1 − 1∕xk).

▸ 44. [M26] What is the sum of all n2 elements in the inverse of Cauchy’s matrix?

For the inverse of the Cauchy matrix

[aij]n−1 = [b ij]n = ∏ 1≤k≤n(xj + yk)(xk + yi) (xj + yi) ∏ 1≤k≤n k≠j (xj − xk) ∏ 1≤k≤n k≠i (yi − yk)n

we want to find the sum of all n2 elements ∑ ⁡ 1≤i,j≤nbij.

But

∑ 1≤i≤nbij = ∑ 1≤i≤n ∏ 1≤k≤n(xj + yk)(xk + yi) (xj + yi) ∏ 1≤k≤n k≠j (xj − xk) ∏ 1≤k≤n k≠i (yi − yk) = ∑ 1≤i≤n ∏ 1≤k≤n(xj + yk) ∏ 1≤k≤n(xk + yi) (xj + yi) ∏ 1≤k≤n k≠j (xj − xk) ∏ 1≤k≤n k≠i (yi − yk) = ∏ 1≤k≤n(xj + yk) ∏ 1≤k≤n k≠j (xj − xk)∑ 1≤i≤n ∏ 1≤k≤n(xk + yi) (xj + yi)∏ 1≤k≤n k≠i (yi − yk) = ∏ 1≤k≤n(xj + yk) ∏ 1≤k≤n k≠j (xj − xk)∑ 1≤i≤n ∏ 1≤k≤n k≠i (xk + yi) ∏ 1≤k≤n k≠i (yi − yk) = ∏ 1≤k≤n(xj + yk) ∏ 1≤k≤n k≠j (xj − xk)

since in general

∑ 1≤i≤n ∏ 1≤k≤n−1(yi − zk) ∏ 1≤k≤n k≠i (yi − yk) = 1

from exercise 34.

Then, for some polynomial P(x) of order n − 2 and by exercise 33

∑ 1≤i,j≤nbij = ∑ 1≤i≤n ∏ 1≤k≤n(xj + yk) ∏ 1≤k≤n k≠j (xj − xk) = ∑ 1≤i≤n xjn + ∑ 1≤k≤nyk xjn−1 + P(xj) ∏ 1≤k≤n k≠j (xj − xk) = ∑ 1≤i≤n xjn ∏ 1≤k≤n k≠j (xj − xk) + ∑ 1≤k≤nyk ∑ 1≤i≤n xjn−1 ∏ 1≤k≤n k≠j (xj − xk) + ∑ 1≤i≤n P(xj) ∏ 1≤k≤n k≠j (xj − xk) = ∑ 1≤k≤nxk + ∑ 1≤k≤nyk + 0 = ∑ 1≤k≤n(xk + yk).

▸ 45. [M25] A Hilbert matrix, sometimes called an n × n segment of the (infinite) Hilbert matrix, is a matrix for which aij = 1∕(i + j − 1). Show that this is a special case of Cauchy’s matrix, find its inverse, show that each element of the inverse is an integer, and show that the sum of all elements of the inverse is n2. [Note: Hilbert matrices have often been used to test various matrix manipulation algorithms, because they are numerically unstable, and they have known inverses. However, it is a mistake to compare the known inverse, given in this exercise, to the computed inverse of a Hilbert matrix, since the matrix to be inverted must be expressed in rounded numbers beforehand; the inverse of an approximate Hilbert matrix will be somewhat different from the inverse of an exact one, due to the instability present. Since the elements of the inverse are integers, and since the inverse matrix is just as unstable as the original, the inverse can be specified exactly, and one should try to invert the inverse. The integers that appear in the inverse are, however, quite large.] The solution to this problem requires an elementary knowledge of factorials and binomial coefficients, which are discussed in Sections 1.2.5 and 1.2.6.

The Hilbert matrix, defined as

An = [aij]n = 1 i + j − 1n = 1 1 2 … ⁡ 1 n 1 2 1 3 … ⁡ 1 n+1 ⋮ ⁡ ⋮ ⁡ ⋱ ⋮ ⁡ 1 n 1 n+1… ⁡ 1 2n−1 n

is a special case of a Cauchy matrix Cn = [cij]n = 1 xi+yj n with xi = i and yj = j − 1.

As such, its inverse [bij]n = 1 i+j−1 n−1 has elements given by

bij = ∏ 1≤k≤n(j + k − 1)(k + i − 1) (j + i − 1) ∏ 1≤k≤n k≠j (j − k) ∏ 1≤k≤n k≠i (i − k) = ∏ 1≤k≤n(k + i − 1)∏ 1≤k≤n(k + i − 1) (j + i − 1) ∏ 1≤k≤j k≠j (j − k) ∏ j≤k≤n k≠j (−1)(k − j) ∏ 1≤k≤i k≠i (i − k) ∏ i≤k≤n k≠i (−1)(k − i) = (j + n − 1)! (j − 1)! (i + n − 1)! (i − 1)! (j + i − 1)(j − 1)! 1! (−1)n−j+1(n − j)! 1! (i − 1)! 1! (−1)n−i+1(n − i)! 1! = (j + n − 1)!(i + n − 1)! (−1)2n−i−j+2(j + i − 1)(j − 1)!(i − 1)!(j − 1)!(n − j)!(i − 1)!(n − i)! = (j + n − 1)!(i + n − 1)! (−1)−(i+j)(j + i − 1)((j − 1)!(i − 1)!)2(n − j)!(n − i)! = (−1)i+j(i + n − 1)!(j + n − 1)! (i + j − 1)(i − 1)!2(j − 1)!2(n − i)!(n − j)!.

It is clear that bij is an integer, as it may be cast into binomial coefficients as

bij = (−1)i+j(i + n − 1)!(j + n − 1)! (i + j − 1)(i − 1)!2(j − 1)!2(n − i)!(n − j)! = (−1)i+j 1 (j − 1)!(i − 1)! (i + n − 1)! (i − 1)! (j + n − 1)! (j + i − 1)(n − i)! 1 (n − j)!(j − 1)! = (−1)i+jj (i + j − 2)! (j − 1)!(i − 1)! (i + n − 1)! n!(i − 1)! (j + n − 1)! (j + i − 1)!(n − i)! n! (n − j)!j! = (−1)i+jj (i + j − 2)! ((i + j − 2) − (i − 1))!(i − 1)! (i + n − 1)! ((i + n − 1) − (i − 1))!(i − 1)! (j + n − 1)! ((j + n − 1) − (n − i))!(n − i)! n! (n − j)!j! = (−1)i+jji + j − 2 i − 1 i + n − 1 i − 1 j + n − 1 n − i n j .

From exercise 44, the sum of the elements of the inverse is simply

∑ 1≤i,j≤nbij = ∑ 1≤k≤n(k + k − 1) = 2 ∑ 1≤k≤nk −∑ 1≤k≤n1 = 2n(n + 1) 2 − n = n2 + n − n = n2.

______________________________________________________________________________________________________________________________

[For further information, see J. Todd, J. Research Nat. Bur. Stand. 65 (1961), 19–22; A. Cauchy, Exercices d’analyse et de physique mathématique 2 (1841), 151–159.]

▸ 46. [M30] Let A be an m × n matrix, and let B be an n × m matrix. Given that 1 ≤ j1,j2,… ⁡,jm ≤ n, let Aj1j2… ⁡jm denote the m × m matrix consisting of columns j1,… ⁡,jm of A, and let Bj1j2… ⁡jm denote the m × m matrix consisting of rows j1,… ⁡,jm of B. Prove the Binet-Cauchy identity

det ⁡ (AB) = ∑ 1≤j1<j2<⋯<jm≤n det ⁡ (Aj1j2… ⁡ jm)det ⁡ (Bj1j2… ⁡ jm).

(Note the special cases: (i) m = n, (ii) m = 1, (iii) B = AT, (iv) m > n, (v) m = 2.)

Proposition. det ⁡ (AB) = ∑ ⁡ 1≤j1<j2<⋯<jm≤n det ⁡ (Aj1j2… ⁡ jm)det ⁡ (Bj1j2… ⁡ jm).

Proof. Let A be an m × n matrix, and let B be an n × m matrix. Given that 1 ≤ j1,j2,… ⁡,jm ≤ n, let Aj1j2… ⁡jm denote the m × m matrix consisting of columns j1,… ⁡,jm of A, and let Bj1j2… ⁡jm denote the m × m matrix consisting of rows j1,… ⁡,jm of B. We must show that

det ⁡ (AB) = ∑ 1≤j1<j2<⋯<jm≤n det ⁡ (Aj1j2… ⁡ jm)det ⁡ (Bj1j2… ⁡ jm).

Let

𝜖(k1,… ⁡,km) = sign ⁡ (∏ 1≤i<j≤m(kj − ki))

for

sign ⁡ (x) = [x > 0] − [x < 0]

be the Levi-Civita function so that in general

det ⁡ ([cij]n) = ∑ 1≤i1,…,in≤n𝜖(i1,…,in)∏ 1≤i≤nai,ii;

and so that if (k1,… ⁡,km) and (l1,… ⁡,lm) are identical except that ki = lj and kj = li, so that 𝜖(k1,… ⁡,km) and − 𝜖(l1,… ⁡,lm), then in general

det ⁡ (Bk1… ⁡ km) = 𝜖(k1,… ⁡ ,km)det ⁡ (Bj1… ⁡ jm)

if j1 ≤⋯ ≤ jm are the numbers k1,… ⁡,km rearranged into nondecreasing order.

Then

det ⁡ (AB) = ∑ 1≤l1,…,lm≤m𝜖(l1,…,lm) ∑ 1≤k≤na1kbkl1 ⋯ ∑ 1≤k≤namkbklm = ∑ 1≤k1,…,km≤na1k1…amkm ∑ 1≤l1,…,lm≤m𝜖(l1,…,lm)bk1l1…bkmlm = ∑ 1≤k1,…,km≤na1k1…amkm det ⁡ (Bk1… ⁡ km) = ∑ 1≤k1,…,km≤n𝜖(k1,…,km)a1k1…amkm det ⁡ (Bj1… ⁡ jm) = ∑ 1≤j1≤⋯≤jm≤n det ⁡ (Aj1… ⁡ jm)det ⁡ (Bj1… ⁡ jm)

as we needed to show.

Note that if m = n, we have the usual identity for square matrices

det ⁡ (AB) = det ⁡ (A)det ⁡ (B);

if m = 1, we have the dot product

det ⁡ (AB) = ∑ 1≤k≤na1kbk1 = A ⋅ B;

if B = AT, we have the square

det ⁡ (AB) = det ⁡ (AAT) = ∑ 1≤j1≤⋯≤jm≤n det ⁡ (Aj1… ⁡ jm)2;

and if m = 2, we have the first nontrivial case of the identity

det ⁡ (AB) = ∑ 1≤j1≤j2≤n det ⁡ (Aj1j2)det ⁡ (Bj1j2).
□

______________________________________________________________________________________________________________________________

[J. de l’École Polytechnique 9 (1813), 280–354; 10 (1815), 29–112. Binet and Cauchy presented their papers on the same day in 1812.]

47. [M27] (C. Krattenthaler.) Prove that

det ⁡ (x + q2)(x + q3)(x + p1)(x + q3)(x + p1)(x + p2) (y + q2)(y + q3)(y + p1)(y + q3)(y + p1)(y + p2) (z + q2)(z + q3)(z + p1)(z + q3)(z + p1)(z + p2) = (x − y)(x − z)(y − z)(p1 − q2)(p1 − q3)(p2 − q3).

and generalize this equation to an identity for an n × n determinant in 3n − 2 variables x1,… ⁡,xn,p1,… ⁡,pn−1,q2,… ⁡,qn. Compare your formula to the result of exercise 38.

We may prove the more general equation.

Proposition. det ⁡ ∏ ⁡ 1≤k≤j−1(xi + pk) ∏ ⁡ j+1≤k≤n(xi + qk)n = ∏ ⁡ 1≤i<j≤n(xi − xj)(pi − qj).

Proof. Let

An = [aij]n = ∏ 1≤k≤j−1(xi + pk) ∏ j+1≤k≤n(xi + qk)n = ∏ 2≤k≤n(x1 + qk)(x1 + p1)∏ 3≤k≤n(x1 + qk)… ⁡∏ 1≤k≤n−1(x1 + pk) ∏ 2≤k≤n(x2 + qk)(x2 + p1)∏ 3≤k≤n(x2 + qk)… ⁡∏ 1≤k≤n−1(x2 + pk) ⋮ ⁡ ⋮ ⁡ ⋱ ⋮ ⁡ ∏ 2≤k≤n(xn + qk)(xn + p1)∏ 3≤k≤n(xn + qk)… ⁡∏ 1≤k≤n−1(xn + pk) n.

We must show that

det ⁡ [aij]n = ∏ 1≤i<j≤n(xi − xj)(pi − qj).

Let

aij′ = ai1 if j = 1 aij − ai,j−1if 2 ≤ j ≤ n

so that det ⁡ [aij]n = det ⁡ [aij′]n where

aij′ = ∏ 1≤k≤j−1(xi + pk) ∏ j+1≤k≤n(xi + qk) if j = 1 (pj−1 − qj) ∏ 1≤k≤j−2(xi + pk) ∏ j+1≤k≤n(xi + qk)if 2 ≤ j ≤ n

since for 2 ≤ j ≤ n

aij′ = ∏ 1≤k≤j−1(xi + pk) ∏ j+1≤k≤n(xi + qk) −∏ 1≤k≤j−2(xi + pk) ∏ j≤k≤n(xi + qk) = (xi + pj−1)∏ 1≤k≤j−2(xi + pk) ∏ j+1≤k≤n(xi + qk) −∏ 1≤k≤j−2(xi + pk) (xi + qj)∏ j+1≤k≤n(xi + qk) = ((xi + pj−1) − (xi + qj)) ∏ 1≤k≤j−2(xi + pk) ∏ j+1≤k≤n(xi + qk) = (pj−1 − qj) ∏ 1≤k≤j−2(xi + pk) ∏ j+1≤k≤n(xi + qk).

Let

qij = 1 if j = 1 pj−1 − qjif 2 ≤ j ≤ n

so that

[aij′] n = ([qij]n[bij]nT)T

and that

det ⁡ [aij′] n = det ⁡ (([qij]n[bij]nT)T) = det ⁡ ([qij]n[bij]nT) = det ⁡ [qij]n det ⁡ ([bij]nT) = ∏ 2≤k≤n(pk−1 − qk)det ⁡ ([bij]nT) = ∏ 2≤k≤n(pk−1 − qk)det ⁡ [bij]n

where

[bij] = ∏ 1≤k≤j−2(xi + pk) ∏ j+1≤k≤n(xi + qk)n.

Repeating this transformation, each time over fewer columns to factor out (pj−2 − qj−1), (pj−3 − qj−2), etc., we eventually have

det ⁡ [aij]n = ∏ 1≤i<j≤n(pi − qj)det ⁡ [bij]n

for

[bij]n = ∏ j+1≤k≤n(xi + qk) n.

Then let

bij′ = bij − (qj+1)bi,j+1if 1 ≤ j ≤ n − 1 b in if j = n

so that det ⁡ [bij]n = det ⁡ [bij′]n where

bij′ = xi ∏ j+2≤k≤n(xi + qk)if 1 ≤ j ≤ n − 1 xi if j = n

since for 1 ≤ j ≤ n − 1

bij′ = ∏ j+1≤k≤n(xi + qk) − (qj+1)∏ j+2≤k≤n(xi + qk) = (xi + qj+1)∏ j+2≤k≤n(xi + qk) − (qj+1)∏ j+2≤k≤n(xi + qk) = ((xi + qj+1) − (qj+1))∏ j+2≤k≤n(xi + qk) = xi ∏ j+2≤k≤n(xi + qk).

Repeating this transformation also, each time over fewer columns to factor out qj+2, qj+3, etc., we eventually have

det ⁡ [aij]n = ∏ 1≤i<j≤n(pi − qj)det ⁡ [bij]n = ∏ 1≤i<j≤n(pi − qj)det ⁡ [bij′] n

for

[bij′] n = xin−j n.

Let

cij′ = b1j′ if i = 1 bij′− b1j′if 2 ≤ i ≤ n

and

cij″ ⁡ = c1j′− x1c1(j+1)′if i = 1,1 ≤ j ≤ n − 1 c1n′ if i = 1,j = n cij′− x1ci(j+1)′if 2 ≤ i ≤ n,1 ≤ j ≤ n − 1 cin′ if 2 ≤ i ≤ n,j = n ≤ n,

so that det ⁡ [bij′] = det ⁡ [cij″ ⁡ ] where

cij″ ⁡ = 0 if i = 1,1 ≤ j ≤ n − 1 1 if i = 1,j = n xin−j−1(xi − x1)if 2 ≤ i ≤ n,1 ≤ j ≤ n − 1 0 if 2 ≤ i ≤ n,j = n ≤ n,

since:

c1j′ = b 1j′ = x 1n−j i = 1,1 ≤ j ≤ n − 1 c1n′ = b 1n′ = x 10 = 1 i = 1,j = n cij′ = b ij′− b 1j′ = x in−j − x 1n−j 2 ≤ i ≤ n,1 ≤ j ≤ n − 1 cin′ = b ij′− b 1j′ = x i0 − x 10 = 0 2 ≤ i ≤ n,j = n ≤ n c1j″ ⁡ = x 1n−j − x 1x1n−j−1 = 0 i = 1,1 ≤ j ≤ n − 1 c1n″ ⁡ = x 10 = 1 i = 1,j = n cij″ ⁡ = x in−j − x 1n−j − x 1(xin−j−1 − x 1n−j−1) = x in−j−1(x i − x1)2 ≤ i ≤ n,1 ≤ j ≤ n − 1 cin″ ⁡ = x i0 − x 10 = 0 2 ≤ i ≤ n,j = n ≤ n

Let

rij = xi+1 − x1if i = j,1 ≤ i,j ≤ n − 1 0 otherwise

so that minor ⁡ ([b′]n,1,n) = [rij]n−1 minor ⁡ ([c″ ⁡ ]n,1,n) and det ⁡ [rij]n−1 = ∏ ⁡ 1≤k≤n−1(xk+1 − x1). Then

det ⁡ [cij″ ⁡ ] n = ∑ 1≤i≤ncin″ cofactor ⁡ (c in″ ⁡ ) = c1n″ ⁡ cofactor ⁡ (c 1n″ ⁡ ) + ∑ 2≤i≤ncin″ cofactor ⁡ (c in″ ⁡ ) = (1)(−1)1+n det ⁡ minor ⁡ ([c″ ⁡ ] n,1,n) + 0 = (−1)n+1 det ⁡ minor ⁡ ([c″ ⁡ ] n,1,n) = (−1)n+1 det ⁡ [r ij]n−1 det ⁡ (minor ⁡ ([b′] n,1,n)) = (−1)n+1 det ⁡ [r ij]n−1 det ⁡ minor ⁡ ([b′] n,1,n) = (−1)n+1 ∏ 1≤k≤n−1(xk+1 − x1)det ⁡ minor ⁡ ([b′] n,1,n) = ∏ 1≤k≤n−1(xk+1 − x1)det ⁡ minor ⁡ ([b′] n,1,n)

since a product of n + 1 signs will cancel with n − 1 (i.e., (−1)n+1(−1)n−1 = (−1)2n = 1). This recursive identity allows us to prove by mathematical induction on n that

det ⁡ [bij′] k = ∏ 1≤i<j≤k(xi − xj).

Therefore

det ⁡ An = det ⁡ [aij]n = det ⁡ ∏ 1≤k≤j−1(xi + pk) ∏ j+1≤k≤n(xi + qk)n = ∏ 1≤i<j≤n(pi − qj)det ⁡ [bij]n = ∏ 1≤i<j≤n(pi − qj)det ⁡ [bij′] n = ∏ 1≤i<j≤n(pi − qj) xin−j n = ∏ 1≤i<j≤n(pi − qj)∏ 1≤i<j≤k(xi − xj) = ∏ 1≤i<j≤n(xi − xj)(pi − qj)

as we needed to show. □

______________________________________________________________________________________________________________________________

[Manuscripta Math. 69 (1990), 177–178.]