În teoria numerelor , numerele Perrin sunt membri ai unei secvențe liniare recurente :
P (0)=3, P (1)=0, P (2)=2,și
P ( n ) = P ( n − 2) + P ( n − 3) pentru n > 2.Secvența numerelor Perrin începe cu
3 , 0 , 2 , 3, 2, 5 , 5, 7 , 10 , 12 , 17 , 22 , 29 , 39 , ... ( secvența OEIS A001608 )Această secvență a fost menționată de Édouard Lucas în 1876. În 1899, Perrin a folosit în mod explicit aceeași secvență. Cel mai aprofundat studiu al acestei secvențe a fost făcut de Adams și Shanks (1982).
Funcția generatoare a numerelor Perrin este
Secvența numerelor Perrin poate fi scrisă în funcție de gradul rădăcinilor ecuației caracteristice
Această ecuație are trei rădăcini. Una dintre aceste rădăcini p este reală (cunoscută ca număr plastic ). Folosind-o și două rădăcini complexe conjugate q și r , se poate exprima numărul Perrin într-un mod similar cu formula lui Binet pentru numerele Lucas :
Deoarece valorile absolute ale rădăcinilor complexe q și r sunt mai mici decât 1, puterile acestor rădăcini vor tinde spre 0 pe măsură ce n crește . Pentru n mare , formula se simplifică la
Această formulă poate fi utilizată pentru a calcula rapid numerele Perrin pentru n mare . Raportul membrilor succesivi ai secvenței Perrin tinde spre p ≈ 1,324718. Această constantă joacă același rol pentru secvența Perrin ca și raportul de aur pentru numerele Lucas . O relație similară există între p și șirul Padovan , între raportul de aur și numerele Fibonacci și între raportul de argint și numerele Pell .
Din formulele lui Binet, putem obține formule pentru G ( kn ) în termeni de G ( n −1), G ( n ) și G ( n +1). Noi stim aia
Ceea ce ne oferă un sistem de trei ecuații liniare cu coeficienți din câmpul de expansiune al polinomului . Calculând matricea inversă, putem rezolva ecuațiile și obținem . Apoi putem ridica toate cele trei valori obținute la puterea k și calcula suma.
Un exemplu de program în sistemul Magma :
P<x> := PolynomialRing(Rationals()); S<t> := SplittingField(x^3-x-1); P2<y> := PolynomialRing(S); p,q,r := Explode([r[1] : r în Roots(y^3-y-1)]); Mi:=Matrice([[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 în [-1..1]] ;Lasă . Ca urmare a rezolvării sistemului, obținem
Numărul 23 apare în aceste formule ca discriminant al polinomului care definește șirul ( ).
Acest lucru permite calcularea celui de-al n-lea număr Perrin în aritmetica întregului folosind înmulțiri.
S-a dovedit că toate numerele prime p împart P ( p ). Reversul, însă, nu este adevărat - unele numere compuse n pot împărți P ( n ), astfel de numere sunt numite numere Perrin pseudo-prime .
Existența pseudoprimelor Perrin a fost considerată de Perrin însuși, dar nu s-a știut dacă au existat sau nu până când Adams și Shanks (1982) l-au descoperit pe cel mai mic dintre ele, 271441 = 521 2 . Următorul a fost . Sunt cunoscute șaptesprezece numere Perrin pseudo-prime, mai puțin de un miliard; [1] Jon Grantham a demonstrat [2] că există o infinitate de pseudoprime Perrin.
Numerele Perrin, care sunt prime , formează șirul:
2 3 5 7 17 29 277 367 853 14197 43721 1442968193 792606555396977 187278659180417234321 662410681948196555396977Weisstein a găsit un prim probabil Perrin P (263226) cu 32147 de cifre în mai 2006 [3] .