Perrin Sequence

DOWNLOAD Mathematica Notebook Contribute to this entry

The integer sequence defined by the recurrence

 P(n)=P(n-2)+P(n-3)
(1)

with the initial conditions P(0)=3, P(1)=0, P(2)=2. This recurrence relation is the same as that for the Padovan sequence but with different initial conditions. The first few terms for n=0, 1, ..., are 3, 0, 2, 3, 2, 5, 5, 7, 10, 12, 17, ... (OEIS A001608).

FoxTrot by Bill Amend

The above cartoon (Amend 2005) shows an unconventional sports application of the Perrin sequence (right panel). (The left two panels instead apply the Fibonacci numbers).

P(n) is the solution of a third-order linear homogeneous recurrence equation having characteristic equation

 x^3-x-1=0.
(2)

Denoting the roots of this equation by alpha, beta, and gamma, with alpha the unique real root, the solution is then

 P(n)=alpha^n+beta^n+gamma^n.
(3)

Here,

 alpha=(x^3-x-1)_1
(4)

is the plastic constant P, which is also given by the limit

 lim_(n->infty)(P(n))/(P(n-1))=P.
(5)

The asymptotic behavior of P(n) is

 P(n)∼alpha^n.
(6)

The first few primes in this sequence are 2, 3, 2, 5, 5, 7, 17, 29, 277, 367, 853, ... (OEIS A074788), which occur for terms n=2, 3, 4, 5, 6, 7, 10, 12, 20, 21, 24, 34, 38, 75, 122, 166, 236, 355, 356, 930, 1042, 1214, 1461, 1622, 4430, 5802, 9092, 16260, 18926, 23698, 40059, 45003, 73807, 91405, 263226, 316872, 321874, 324098, ... (OEIS A112881), the largest of which are probable primes and a number of which are summarized in the following table.

ndecimal digitsdiscovererdate
9140511163E. W. WeissteinOct. 6, 2005
26322632147E. W. WeissteinMay 4, 2006
31687238698E. W. WeissteinFeb. 4, 2007
32187439309E. W. WeissteinFeb. 19, 2007
32409839580E. W. WeissteinFeb. 25, 2007
58113270970E. W. WeissteinFeb. 15, 2011

Perrin (1899) investigated the sequence and noticed that if n is prime, then n|P(n) (i.e., n divides P(n)). The first statement of this fact is attributed to É. Lucas in 1876 by Stewart (1996). Perrin also searched for but did not find any composite number n in the sequence such that n|P(n). Such numbers are now known as Perrin pseudoprimes. Malo (1900), Escot (1901), and Jarden (1966) subsequently investigated the series and also found no Perrin pseudoprimes. Adams and Shanks (1982) subsequently found that 271441 is such a number.

Wolfram Web Resources

Mathematica »

The #1 tool for creating Demonstrations and anything technical.

Wolfram|Alpha »

Explore anything with the first computational knowledge engine.

Wolfram Demonstrations Project »

Explore thousands of free applications across science, mathematics, engineering, technology, business, art, finance, social sciences, and more.

Computerbasedmath.org »

Join the initiative for modernizing math education.

Online Integral Calculator »

Solve integrals with Wolfram|Alpha.

Step-by-step Solutions »

Walk through homework problems step-by-step from beginning to end. Hints help you try the next step on your own.

Wolfram Problem Generator »

Unlimited random practice problems and answers with built-in Step-by-step solutions. Practice online or make a printable study sheet.

Wolfram Education Portal »

Collection of teaching and learning tools built by Wolfram education experts: dynamic textbook, lesson plans, widgets, interactive Demonstrations, and more.

Wolfram Language »

Knowledge-based programming for everyone.