Numero di Perrin

In matematica, i numeri di Perrin sono definiti dalla relazione di ricorrenza

P ( 0 ) = 3 , P ( 1 ) = 0 , P ( 2 ) = 2 {\displaystyle P(0)=3,\,P(1)=0,\,P(2)=2} ,

e

P ( n ) = P ( n 2 ) + P ( n 3 ) {\displaystyle P(n)=P(n-2)+P(n-3)} per n > 2 {\displaystyle n>2} .

La sequenza dei numeri di Perrin inizia con

3, 0, 2, 3, 2, 5, 5, 7, 10, 12, 17, 22, 29, 39 ...[1]

Il numero dei diversi insiemi indipendenti massimali in un grafo ciclo con n {\displaystyle n} vertici è conteggiato dal numero Perrin n {\displaystyle n} -esimo per n > 1 {\displaystyle n>1} .[2]

Storia

Questa sequenza era già stata menzionata implicitamente da Édouard Lucas (1876). Nel 1899, la stessa sequenza fu menzionata esplicitamente da François Olivier Raoul Perrin.[3] La trattazione più estesa di questa sequenza fu fatta da Adams e Shanks (1982).

Proprietà

Funzione generatrice

La funzione generatrice della sequenza di Perrin è

G ( P ( n ) ; x ) = 3 x 2 1 x 2 x 3 . {\displaystyle G(P(n);x)={\frac {3-x^{2}}{1-x^{2}-x^{3}}}.}

Formula della matrice

( 0 1 0 0 0 1 1 1 0 ) n ( 3 0 2 ) = ( P ( n ) P ( n + 1 ) P ( n + 2 ) ) {\displaystyle {\begin{pmatrix}0&1&0\\0&0&1\\1&1&0\end{pmatrix}}^{n}{\begin{pmatrix}3\\0\\2\end{pmatrix}}={\begin{pmatrix}P\left(n\right)\\P\left(n+1\right)\\P\left(n+2\right)\end{pmatrix}}}

Formula simile di Binet

La sequenza dei numeri di Perrin può essere scritta in termini di potenze delle radici dell'equazione

x 3 x 1 = 0. {\displaystyle x^{3}-x-1=0.}

Questa equazione ha 3 radici; una radice reale p {\displaystyle p} (nota come il numero plastico) e due radici complesse coniugate q {\displaystyle q} ed r {\displaystyle r} . Date queste tre radici, l'analogo per la sequenza di Perrin della formula di Binet per la sequenza di Lucas è

P ( n ) = p n + q n + r n . {\displaystyle P\left(n\right)={p^{n}}+{q^{n}}+{r^{n}}.}

Poiché le grandezze delle radici complesse q {\displaystyle q} ed r {\displaystyle r} sono entrambe minori di 1 {\displaystyle 1} , le potenze di queste radici si avvicinano a 0 {\displaystyle 0} per n {\displaystyle n} grande. Per n {\displaystyle n} grande la formula si riduce a:

P ( n ) p n {\displaystyle P\left(n\right)\approx {p^{n}}}

Questa formula può essere usata per calcolare rapidamente i valori della sequenza di Perrin per n {\displaystyle n} grande. Il rapporto di termini successivi nella sequenza di Perrin si avvicina a p {\displaystyle p} , ossia al numero plastico, che ha un valore approssimativamente di 1,324718. Questa costante ha con la sequenza di Perrin la stessa relazione che la sezione aurea ha con la sequenza di Lucas. Collegamenti simili esistono anche tra p {\displaystyle p} e la sequenza di Padovan, tra la sezione aurea e la successione di Fibonacci, e tra la sezione argentea e i numeri di Pell.

Formula della moltiplicazione

Dalla formula di Binet, possiamo ottenere una formula per G ( k n ) {\displaystyle G(kn)} in termini di G ( n 1 ) , G ( n ) {\displaystyle G(n-1),\,G(n)} e G ( n + 1 ) {\displaystyle G(n+1)} ; sappiamo che:

G ( n 1 ) = p 1 p n + q 1 q n + r 1 r n G ( n ) = p n + q n + r n G ( n + 1 ) = p p n + q q n + r r n {\displaystyle {\begin{matrix}G(n-1)&=&p^{-1}p^{n}+&q^{-1}q^{n}+&r^{-1}r^{n}\\G(n)&=&p^{n}+&q^{n}+&r^{n}\\G(n+1)&=&pp^{n}+&qq^{n}+&rr^{n}\end{matrix}}}

che ci dà tre equazioni lineari con i coefficienti sul campo di spezzamento di x 3 x 1 {\displaystyle x^{3}-x-1} ; invertendo una matrice possiamo risolvere per p n , q n , r n {\displaystyle p^{n},q^{n},r^{n}} e poi possiamo elevarle alla k {\displaystyle k} -esima e computare la somma.

Esempio in codice magma:

P<x> := PolynomialRing(Rationals());
S<t> := SplittingField(x^3-x-1); 
P2<y> := PolynomialRing(S);
p,q,r := Explode([r[1] : r in Roots(y^3-y-1)]);
Mi:=Matrix([[1/p,1/q,1/r],[1,1,1],[p,q,r]])^(-1);
T<u,v,w> := PolynomialRing(S,3);
v1 := ChangeRing(Mi,T) *Matrix([[u],[v],[w]]);
[p^i*v1[1,1]^3 + q^i*v1[2,1]^3 + r^i*v1[3,1]^3 : i in [-1..1]];

con il risultato che, se abbiamo u = G ( n 1 ) , v = G ( n ) , w = G ( n + 1 ) {\displaystyle u=G(n-1),v=G(n),w=G(n+1)} , allora

23 G ( 2 n 1 ) = 4 u 2 + 3 v 2 + 9 w 2 + 18 u v 12 u w 4 v w 23 G ( 2 n ) = 6 u 2 + 7 v 2 2 w 2 4 u v + 18 u w + 6 v w 23 G ( 2 n + 1 ) = 9 u 2 + v 2 + 3 w 2 + 6 u v 4 u w + 14 v w 23 G ( 3 n 1 ) = ( 4 u 3 + 2 v 3 w 3 + 9 ( u v 2 + v w 2 + w u 2 ) + 3 v 2 w + 6 u v w ) 23 G ( 3 n ) = ( 3 u 3 + 2 v 3 + 3 w 3 3 ( u v 2 + u w 2 + v w 2 + v u 2 ) + 6 v 2 w + 18 u v w ) 23 G ( 3 n + 1 ) = ( v 3 w 3 + 6 u v 2 + 9 u w 2 + 6 v w 2 + 9 v u 2 3 w u 2 + 6 w v 2 6 u v w ) {\displaystyle {\begin{matrix}23G(2n-1)&=&4u^{2}+3v^{2}+9w^{2}+18uv-12uw-4vw\\23G(2n)&=&-6u^{2}+7v^{2}-2w^{2}-4uv+18uw+6vw\\23G(2n+1)&=&9u^{2}+v^{2}+3w^{2}+6uv-4uw+14vw\\23G(3n-1)&=&\left(-4u^{3}+2v^{3}-w^{3}+9(uv^{2}+vw^{2}+wu^{2})+3v^{2}w+6uvw\right)\\23G(3n)&=&\left(3u^{3}+2v^{3}+3w^{3}-3(uv^{2}+uw^{2}+vw^{2}+vu^{2})+6v^{2}w+18uvw\right)\\23G(3n+1)&=&\left(v^{3}-w^{3}+6uv^{2}+9uw^{2}+6vw^{2}+9vu^{2}-3wu^{2}+6wv^{2}-6uvw\right)\end{matrix}}}

Il numero 23 qui sorge dalla discriminante del polinomio definitore della sequenza.

Questo permette di computare il numero di Perrin n {\displaystyle n} -esimo usando l'aritmetica degli interi in O ( log n ) {\displaystyle O(\log n)} multipli.

Primi e divisibilità

Pseudoprimi di Perrin

È stato dimostrato che per tutti i numeri primi p {\displaystyle p} , p {\displaystyle p} divide P ( p ) {\displaystyle P(p)} . Tuttavia, l'inverso non è vero: per alcuni numeri composti n {\displaystyle n} , n {\displaystyle n} può dividere P ( n ) {\displaystyle P(n)} . Se n {\displaystyle n} ha questa proprietà, è chiamato pseudoprimo di Perrin.

La questione dell'esistenza degli pseudoprimi di Perrin fu considerata dallo stesso Perrin, ma non si sapeva se esistessero finché Adams e Shanks (1982) non scoprirono il più piccolo, 271441 = 5212; il successivo più piccolo è 904631 = 7 x 13 x 9941. Ce ne sono diciassette minori di un miliardo;[4] Jon Grantham ha dimostrato[5] che ci sono infiniti pseudoprimi di Perrin.

Primi di Perrin

Un primo di Perrin è un numero di Perrin che è un numero primo. I primi iniziali di Perrin sono:

2, 3, 5, 7, 17, 29, 277, 367, 853, 14197, 43721, 1442968193, 792606555396977, 187278659180417234321, 66241160488780141071579864797 [6]

Eric W. Weisstein trovò un probabile primo di Perrin a 32.147 cifre, P(263226), nel maggio 2006.

Note

  1. ^ (EN) Sequenza A001608, su On-Line Encyclopedia of Integer Sequences, The OEIS Foundation.
  2. ^ Füredi (1987)
  3. ^ Knuth (2011)
  4. ^ (EN) Sequenza A013998, su On-Line Encyclopedia of Integer Sequences, The OEIS Foundation.
  5. ^ (EN) Jon Grantham, There are infinitely many Perrin pseudoprimes (PDF), in Journal of Number Theory, vol. 130, n. 5, 2010, pp. 1117–1128, DOI:10.1016/j.jnt.2009.11.008.
  6. ^ (EN) Sequenza A074788, su On-Line Encyclopedia of Integer Sequences, The OEIS Foundation.

Bibliografia

  • (EN) Adams, William; Shanks, Daniel, Strong primality tests that are not sufficient, in Mathematics of Computation, vol. 39, n. 159, American Mathematical Society, 1982, pp. 255–300, DOI:10.2307/2007637, JSTOR 2007637, MR 0658231.
  • (EN) Füredi, Zoltán, The number of maximal independent sets in connected graphs, in Journal of Graph Theory, vol. 11, n. 4, 1987, pp. 463–470, DOI:10.1002/jgt.3190110403.
  • (EN) Knuth, Donald E., The Art of Computer Programming, Volume 4A: Combinatorial Algorithms, Part 1, Addison-Wesley, 2011, ISBN 0-201-03804-8.
  • (FR) Lucas, Édouard, Théorie des fonctions numériques simplement périodiques, in American Journal of Mathematics, vol. 1, n. 3, The Johns Hopkins University Press, 1878, pp. 197–240, DOI:10.2307/2369311, JSTOR 2369311.
  • (FR) Perrin, R., Query 1484, in L'Intermédiaire des Mathématiciens, vol. 6, 1899, p. 76.

Collegamenti esterni

  • (EN) Eric W. Weisstein, Numero di Perrin, su MathWorld, Wolfram Research. Modifica su Wikidata
  • (DE) Zentrum für Hirnforschung Institut für Medizinische Kybernetik und Artificial Intelligence, su ai.univie.ac.at. URL consultato il 23 marzo 2014 (archiviato dall'url originale l'8 novembre 2005).
  • (EN) MathPages - Lucas Pseudoprimes, su mathpages.com.
  • (EN) MathPages - Perrin's Sequence, su mathpages.com.
  • (EN) Perrin-like sequence, su oeis.org.
  Portale Matematica: accedi alle voci di Wikipedia che trattano di Matematica