Pell number

From Wikipedia, the free encyclopedia

In mathematics, the Pell numbers and companion Pell numbers (Pell-Lucas numbers) are both sequences of integers that have been known since ancient times. They are defined by a recurrence relation similar to that for the Fibonacci numbers, and grow exponentially, proportionally to powers of the silver ratio. Pell numbers arise in the approximation of the square root of 2, in the definition of square triangular numbers, in the construction of nearly-isosceles integer right triangles, and in certain combinatorial enumeration problems.[1]

As with Pell's equation, the name of the Pell numbers stems from Leonhard Euler's mistaken attribution of the equation and the numbers derived from it to John Pell. The Pell-Lucas numbers are also named after Edouard Lucas, who studied sequences defined by recurrences of this type; the Pell and companion Pell numbers are Lucas sequences.

Contents

[edit] Pell numbers

The Pell numbers are defined by the recurrence relation

P_n=\begin{cases}0&\mbox{if }n=0;\\1&\mbox{if }n=1;\\2P_{n-1}+P_{n-2}&\mbox{otherwise.}\end{cases}

In words, the sequence of Pell numbers starts with 0 and 1, and then each Pell number is the sum of twice the previous Pell number and the Pell number before that. The first few terms of the sequence are

0, 1, 2, 5, 12, 29, 70, 169, 408, 985, 2378... (sequence A000129 in OEIS).

The Pell numbers can also be expressed by the closed form formula

P_n=\frac{(1+\sqrt2)^n-(1-\sqrt2)^n}{2\sqrt2}.

For large values of n, the \scriptstyle (1+\sqrt 2)^n term dominates this expression, so the Pell numbers are approximately proportional to powers of the silver ratio \scriptstyle (1+\sqrt 2), analogous to the growth rate of Fibonacci numbers as powers of the golden ratio.

A third definition is possible, from the matrix formula

\begin{pmatrix} P_{n+1} & P_n \\ P_n & P_{n-1} \end{pmatrix} = \begin{pmatrix} 2 & 1 \\ 1 & 0 \end{pmatrix}^n.

Many identities can be derived or proven from these definitions; for instance an identity analogous to Cassini's identity for Fibonacci numbers,

P_{n+1}P_{n-1}-P_n^2 = (-1)^n,

is an immediate consequence of the matrix formula (found by considering determinants).[2]

[edit] Approximation to the square root of two

Rational approximations to regular octagons, with coordinates derived from the Pell numbers.
Rational approximations to regular octagons, with coordinates derived from the Pell numbers.

Pell numbers arise historically and most notably in the rational approximation to the square root of 2. If two large integers x and y form a solution to the Pell equation

\displaystyle x^2-2y^2=\pm 1,

then their ratio \tfrac{x}{y} provides a close approximation to \scriptstyle\sqrt 2. The sequence of approximations of this form is

1, \frac32, \frac75, \frac{17}{12}, \frac{41}{29}, \frac{99}{70}, \dots

where the denominator of each fraction is a Pell number and the numerator is the sum of a Pell number and its predecessor in the sequence. That is, the solutions have the form \tfrac{P_{n-1}+P_n}{P_n}. The approximation

\sqrt 2\approx\frac{577}{408}

of this type was known to Indian mathematicians in the third or fourth century B.C.[3] The Greek mathematicians of the fifth century B.C. also knew of this sequence of approximations[4]; they called the denominators and numerators of this sequence side and diameter numbers and the numerators were also known as rational diagonals or rational diameters.[5]

These approximations can be derived from the continued fraction expansion of \scriptstyle\sqrt 2:

\sqrt 2 = 1 + \cfrac{1}{2 + \cfrac{1}{2 + \cfrac{1}{2 + \cfrac{1}{2 + \cfrac{1}{\ddots\,}}}}}.

Truncating this expansion to any number of terms produces one of the Pell-number-based approximations in this sequence; for instance,

\frac{577}{408} = 1 + \cfrac{1}{2 + \cfrac{1}{2 + \cfrac{1}{2 + \cfrac{1}{2 + \cfrac{1}{2 + \cfrac{1}{2 + \cfrac{1}{2}}}}}}}.

As Knuth (1994) describes, the fact that Pell numbers approximate \scriptstyle\sqrt 2 allows them to be used for accurate rational approximations to a regular octagon with vertex coordinates (\pm P_i,\pm P_{i+1}) and (\pm P_{i+1},\pm P_i). All vertices are equally distant from the origin, and form nearly uniform angles around the origin. Alternatively, the points (\pm(P_i+P_{i-1}),0), (0,\pm(P_i+P_{i-1})), and (\pm P_i,\pm P_i) form approximate octagons in which the vertices are nearly equally distant from the origin and form uniform angles.

[edit] Primes and squares

A Pell prime is a Pell number that is prime. The first few Pell primes are

2, 5, 29, 5741, ... (sequence A086383 in OEIS).

As with the Fibonacci numbers, a Pell number Pn can only be prime if n itself is prime.

The only Pell numbers that are squares, cubes, or any higher power of another integer are 0, 1, and 169 = 132.[6]

However, despite having so few squares or other powers, Pell numbers have a close connection to square triangular numbers.[7] Specifically, these numbers arise from the following identity of Pell numbers:

\bigl((P_{k-1}+P_k)\cdot P_k\bigr)^2 = \frac{(P_{k-1}+P_k)^2\cdot\left((P_{k-1}+P_k)^2-(-1)^k\right)}{2}.

The left side of this identity describes a square number, while the right side describes a triangular number, so the result is a square triangular number.

Santana and Diaz-Barrero (2006) prove another identity relating Pell numbers to squares and showing that the sum of the Pell numbers up to P4n + 1 is always a square:

\sum_{i=0}^{4n+1} P_i = \left(\sum_{r=0}^n 2^r{2n+1\choose 2r}\right)^2 = (P_{2n}+P_{2n+1})^2.

For instance, the sum of the Pell numbers up to P5, 0 + 1 + 2 + 5 + 12 + 29 = 49, is the square of P2 + P3 = 2 + 5 = 7. The numbers P2n + P2n + 1 forming the square roots of these sums,

1, 7, 41, 239, 1393, 8119, 47321, ... (sequence A002315 in OEIS),

are known as the NSW numbers.

[edit] Pythagorean triples

Integer right triangles with nearly equal legs, derived from the Pell numbers.
Integer right triangles with nearly equal legs, derived from the Pell numbers.

If a right triangle has integer side lengths a, b, c (necessarily satisfying the Pythagorean theorem a2+b2=c2), then (a,b,c) is known as a Pythagorean triple. As Martin (1875) describes, the Pell numbers can be used to form Pythagorean triples in which a and b are one unit apart, corresponding to right triangles that are nearly isosceles. Each such triple has the form

(2P_{n}P_{n+1}, P_{n+1}^2 - P_{n}^2, P_{n+1}^2 + P_{n}^2=P_{2n+1}).

The sequence of Pythagorean triples formed in this way is

(4,3,5), (20,21,29), (120,119,169), (696,697,985), ...

[edit] Companion Pell numbers (Pell-Lucas numbers)

The companion Pell numbers or Pell-Lucas numbers are defined by the recurrence relation

Q_n=\begin{cases}2&\mbox{if }n=0;\\2&\mbox{if }n=1;\\2Q_{n-1}+Q_{n-2}&\mbox{otherwise.}\end{cases}

In words: the first two numbers in the sequence are both 2, and each successive number is formed by adding twice the previous Pell-Lucas number to the Pell-Lucas number before that. The first few terms of the sequence are (sequence A002203 in OEIS): 2, 2, 6, 14, 34, 82, 198, 478...

The companion Pell numbers can be expressed by the closed form formula

Q_n=(1+\sqrt 2)^n+(1-\sqrt 2)^n.

These numbers are all even; each such number is twice the numerator in one of the rational approximations to \scriptstyle\sqrt 2 discussed above.

[edit] Notes

  1. ^ For instance, Sellers (2002) proves that the number of perfect matchings in the Cartesian product of a path graph and the graph K4-e can be calculated as the product of a Pell number with the corresponding Fibonacci number.
  2. ^ For the matrix formula and its consequences see Ercolano (1979) and Kilic and Tasci (2005). Additional identities for the Pell numbers are listed by Horadam (1971) and Bicknell (1975).
  3. ^ As recorded in the Shulba Sutras; see e.g. Dutka (1986), who cites Thibaut (1875) for this information.
  4. ^ See Knorr (1976) for the fifth century date, which matches Proclus' claim that the side and diameter numbers were discovered by the Pythagoreans. For more detailed exploration of later Greek knowledge of these numbers see Thompson (1929), Vedova (1951), Ridenhour (1986), Knorr (1998), and Filep (1999).
  5. ^ For instance, as several of the references from the previous note observe, in Plato's Republic there is a reference to the "rational diameter of 5", by which Plato means 7, the numerator of the approximation 7/5 of which 5 is the denominator.
  6. ^ Pethő (1992); Cohn (1996). Although the Fibonacci numbers are defined by a very similar recurrence to the Pell numbers, Cohn writes that an analogous result for the Fibonacci numbers seems much more difficult to prove.
  7. ^ Sesskin (1962). See the square triangular number article for a more detailed derivation.

[edit] References

  • Cohn, J. H. E. (1996). "Perfect Pell powers". Glasgow Mathematical Journal 38 (1): 19–20. MR1373953. 
  • Dutka, Jacques (1986). "On square roots and their representations". Archive for History of Exact Sciences 36 (1): 21–39. doi:10.1007/BF00357439. MR0863340. 
  • Kilic, Emrah; Tasci, Dursun (2005). "The linear algebra of the Pell matrix". Boletín de la Sociedad Matemática Mexicana, Tercera Serie 11 (2): 163–174. MR2207722. 
  • Knorr, Wilbur (1976). "Archimedes and the measurement of the circle: A new interpretation". Archive for History of Exact Sciences 15 (2): 115–140. doi:10.1007/BF00348496. MR0497462. 
  • Pethő, A. (1992). "The Pell sequence contains only trivial perfect powers". Sets, graphs, and numbers (Budapest, 1991): 561–568, Colloq. Math. Soc. János Bolyai, 60, North-Holland. MR1218218. 
  • Thibaut, George (1875). "On the Súlvasútras". Journal of the Royal Asiatic Society of Bengal 44: 227–275. 

[edit] External links