Suite de Sylvester

Représentation graphique de la convergence vers 1 de la somme 1/2 + 1/3 + 1/7 + 1/43 +... Chaque rang dencarrés de côté 1/na une aire totale de 1/n;l'ensemble des rangs recouvre exactement un carré plus grand, d'aire 1. (Les carrés de côté 1/1807 ou moins sont trop petits pour être dessinés sur le graphique.)

Enthéorie des nombres,lasuite de Sylvesterest unesuited'entierstelle que chaque terme est le produit de tous les termes précédents augmenté de 1, en partant d'un terme initial égal à 2. Les premiers termes de la suite sont:

2;3;7;43;1 807; 3 263 443; 10 650 056 950 807; 113 423 713 055 421 844 361 000 443 (Voir la suiteA000058de l'OEIS).

En hommage à la démonstration par Euclide de l'infinitude desnombres premiers,les termes de cette suite sont aussi parfois appelés "nombres d'Euclide".

La suite de Sylvester doit son nom àJames Joseph Sylvesterqui, le premier, étudia ses propriétés dans les années 1880. Ses termes présentent unecroissance exponentielle double.La série formée de la somme des inverses de cette suite converge vers 1, plus vite que toute autre série somme infinie d'inverses d'entiers convergeant vers 1.

La relation de récurrence qui définit les termes de la suite permet defactoriserceux-ci plus facilement que toute autre série de croissance comparable, mais, du fait de la croissance rapide de la série, la décomposition en nombres premiers n'est connue que pour quelques termes. Des valeurs extraites de cette suite ont été utilisées pour construire des représentations de 1 sous forme dedéveloppement en fractions égyptiennes,et intervient dans l'étude desvariétés d'Einsteinsasakiennes(en).

Définitions

La suite de Sylvester est définie parrécurrenceforte:

$s_{0}=2,\ \forall n\geqslant 0,\,s_{n+1}=1+\prod _{k=0}^{n}s_{k}$

On vérifie alors que$s_{n+1}-1=\prod _{k=0}^{n}s_{k}=(s_{n}-1)s_{n}$,si bien qu'on peut également la définir par récurrence simple:

$s_{0}=2,\ \forall n\geqslant 0,\,s_{n+1}=s_{n}(s_{n}-1)+1\quad [1]$

On en déduit aussi unedéfinition par récurrencedouble:

$s_{0}=2,s_{1}=3,\,\forall n\geqslant 0\,\,s_{n+2}=s_{n+1}^{2}-s_{n}(s_{n}-1)$

On obtient une suite strictement croissante d'entiers, donc de limite infinie.

Série des inverses et lien avec les fractions égyptiennes

Comme${\frac {s_{n+1}}{s_{n}}}=s_{n}-1+{1 \over {s_{n}}}\sim s_{n}\longrightarrow +\infty$,la série de terme général$1/s_{n}$est convergente d'après lecritère de D'Alembert.

Un résultat remarquable est que la somme de cette série est égale à 1.

En effet, il résulte de la relation de récurrence [1] que:

${\frac {1}{s_{n}-1}}-{\frac {1}{s_{n+1}-1}}={\frac {1}{s_{n}}}$;

La somme des termes de$1/s_{0}$à$1/s_{N}$se réduit alors, partéléscopage,à:

$\sum _{n=0}^{N}{\frac {1}{s_{n}}}=\sum _{n=0}^{N}\left({\frac {1}{s_{n}-1}}-{\frac {1}{s_{n+1}-1}}\right)={\frac {1}{s_{0}-1}}-{\frac {1}{s_{N+1}-1}}=1-{\frac {1}{s_{N+1}-1}}\quad [2].$

Puisque la suite$(s_{n})$tend vers l'infini,$\sum _{n=0}^{N}{\frac {1}{s_{n}}}$tend bien vers 1.

On peut donc écrire:

$1=\sum _{n=0}^{\infty }{\frac {1}{s_{n}}}={\frac {1}{2}}+{\frac {1}{3}}+{\frac {1}{7}}+{\frac {1}{43}}+{\frac {1}{1807}}+\cdots$

ce qui donne une représentation de 1 sous forme de somme infinie defractions égyptiennes.

La même formule écrite sous la forme${\frac {1}{2}}={\frac {1}{3}}+{\frac {1}{7}}+{\frac {1}{43}}+{\frac {1}{1807}}+\cdots$fournit un développement de1/2en somme infinie de fractions égyptiennes de dénominateurs impairs.

De plus, la formule [2] écrite sous la forme:$1=\sum _{n=0}^{N}{\frac {1}{s_{n}}}+{\frac {1}{s_{N+1}-1}}$donne une représentation de 1 sous forme de somme de fractions égyptiennes de longueur quelconque (tronquer la somme infinie après un nombre arbitraire de termes et soustraire 1 du dénominateur de la dernière fraction); par exemple:

$1={\tfrac {1}{2}}+{\tfrac {1}{3}}+{\tfrac {1}{6}},\quad 1={\tfrac {1}{2}}+{\tfrac {1}{3}}+{\tfrac {1}{7}}+{\tfrac {1}{42}},\quad 1={\tfrac {1}{2}}+{\tfrac {1}{3}}+{\tfrac {1}{7}}+{\tfrac {1}{43}}+{\tfrac {1}{1806}},\quad \dots$

La somme des$k$premiers termes de la suite$\left({1 \over s_{n}}\right)$constitue la meilleure approximation possible par défaut de 1 à l'aide de$k$fractions égyptiennes.Par exemple, la somme des quatre premiers termes de la série vaut 1805/1806 et par conséquent, pour représenter tout nombre dans l'intervalle ouvert ]1805/1806, 1[, une somme de fractions égyptiennes doit avoir au moins cinq termes. Pour cette raison, toute série d'inverses d'entiers convergeant vers 1 converge moins vite que la série des inverses de la suite de Sylvester.

On peut interpréter la suite de Sylvester comme le résultat de l'algorithme gloutonpour la décomposition du nombre 1 en somme de fractions égyptiennes, algorithme qui, à chaque étape, choisit la plus grande fraction égyptienne dont la somme avec les précédentes est strictement inférieure à 1.

Plus précisément si on pose$s_{0}=2,s_{n}=\min {\left\{m\in \mathbb {N} /\sum _{k=0}^{n-1}{1 \over s_{k}}+{1 \over m}<1\right\}}$,on obtient de nouveau la suite de Sylvester.

Étude asymptotique

Comme$s_{n+1}\sim s_{n}^{2}$,la suite de Sylvester présente une croissancedoublement exponentielle.

La suite$\left(s_{n}^{1/2^{n+1}}\right)$tend en décroissant vers un nombreE= 1,264084735305...appelé constante de Vardi (suiteA076393de l'OEIS)et on montre que non seulement$s_{n}\sim E^{2^{n+1}}$mais qu'on peut calculer exactement les termes de la suite de Sylvester par la formule:

$s_{n}=\left\lfloor E^{2^{n+1}}\right\rceil$

où$\left\lfloor x\right\rceil$désigne l'entier le plus prochede$x$.

La croissance doublement exponentielle de la suite de Sylvester est comparable à celle de la suite desnombres de Fermat$F_{n}$.Ceux-ci sont habituellement définis par l'expression en double exponentielle:$2^{2^{n}}+1$,mais ils sont aussi définis par une récurrence très voisine de celle définissant la suite de Sylvester:

$F_{n+1}=2+\prod _{k=0}^{n}{F_{k}}$,alors que$s_{n+1}=1+\prod _{k=0}^{n}s_{k}$.

Divisibilité et factorisations

Si$i<j$,il résulte de la définition que$s_{j}\equiv 1\mod s_{i}$.Par conséquent, deux termes quelconques de la suite de Sylvester sontpremiers entre eux.La suite peut donc servir de preuve à l'assertion "il existe une infinité de nombres premiers", puisqu'un nombre premier donné ne peut diviser qu'un terme de la suite au plus.

Plusieurs travaux ont été consacrés à la factorisation des termes de la suite de Sylvester en nombres premiers, mais il demeure beaucoup d'incertitudes à ce sujet. Par exemple, on ne sait pas si tous les termes de la suite sontsans facteurs carrés,bien que tous les termes connus le soient.

CommeVardi (1991)l'indique, il est facile de déterminer de quel terme de la suite de Sylvester (s'il existe) un nombre premier$p$donné est diviseur: il suffit de calculer les termes de la suite modulo$p$à l'aide de la définition par récurrence (en répétant "remplacer$x$par$x^{2}-x+1\mod p$") jusqu'à trouver un terme nul. Et si l'on obtient un terme non nul qui avait déjà été trouvé, on est assuré que$p$ne divise aucun terme de la suite. À l'aide de cette technique, il a obtenu qu'exactement 1166 des trois premiers millions de nombres premiers sont des diviseurs des nombres de Sylvesteret qu'aucun d'entre eux n'était élevé au carré.

Un résultat deJones (2006)conclut que la densité dans l'ensemble des nombres premiers des diviseurs premiers des termes de la suite de Sylvester est nulle.

D'autre part Odonia montré que tous les diviseurs premiers des termes de la suite de Sylvester (excepté 2 et 3) sont de la forme 6k+1.

La table ci-après présente la factorisation des termes de la suite de Sylvester (à l'exception des quatre premiers termes qui sont tous des nombres premiers):

(la notation Pn représente un nombre de$n$chiffres dont on sait qu'il est premier, et Cn un nombre de$n$chiffres dont on sait qu'il est composé, mais dont la décomposition est inconnue)

n	Facteurs desn
4	13 × 139
5	3263443, premier
6	547 × 607 × 1033 × 31051
7	29881 × 67003 × 9 119 521 × 6 212 157 481
8	5 295 435 634 831 × 31 401 519 357 481 261 × 77 366 930 214 021 991 992 277
9	181 × 1987 × 112 374 829 138 729 × 114 152 531 605 972 711 × P68
10	2287 × 2 271 427 × 21 430 986 826 194 127 130 578 627 950 810 640 891 005 487 × P156
11	73 × C416
12	2 589 377 038 614 498 251 653 × 2 872 413 602 289 671 035 947 763 837 × C785
13	52387 × 5 020 387 × 5 783 021 473 × 401 472 621 488 821 859 737 × 287 001 545 675 964 617 409 598 279 × C1600
14	13999 × 74203 × 9 638 659 × 57 218 683 × 10 861 631 274 478 494 529 × C3293
15	17881 × 97 822 786 011 310 111 × 54 062 008 753 544 850 522 999 875 710 411 × C6618
16	128551 × C13335
17	635263 × 1 286 773 × 21 269 959 × C26661
18	50 201 023 123 × 139 263 586 549 × 60 466 397 701 555 612 333 765 567 × C53313
19	C106721
20	352867 × 6 210 298 470 888 313 × C213419
21	387 347 773 × 1 620 516 511 × C426863
22	91 798 039 513 × C853750

La liste croissante des nombres premiers qui divisent un terme de la suite de Sylvester est répertoriéeA007996dans l'OEIS.

Généralisation: développement en série de Sylvester d'un réel quelconque

On peut généraliser l'algorithme gloutonde décomposition en somme infinie defractions égyptiennesà un réel strictement positif$x$autre que 1: cet algorithme choisit, à chaque étape, la plus grande fraction égyptienne dont la somme avec les précédentes est strictement inférieure à$x$.

On pose donc$s_{n}=\min {\left\{m\in \mathbb {N} /\sum _{k=0}^{n-1}{1 \over s_{k}}+{1 \over m}<x\right\}}$,ou, de façon plus pratique:

$x_{-1}=x$puis pour tout entier naturel$n$:$\left\{{\begin{matrix}{s_{n}=\lfloor {\frac {1}{x_{n-1}}}\rfloor +1}\\{x_{n}=x_{n-1}-{1 \over s_{n}}}\end{matrix}}\right.$.

Ledéveloppement en série deSylvesterde$x$s'écrit alors:

$x={\frac {1}{s_{0}}}+{\frac {1}{s_{1}}}+{\frac {1}{s_{2}}}+\dotsb$que l'on notera$]s_{0},s_{1},...[$.

La suite d'entiers naturels non nuls$(s_{n})_{n\geqslant 0}$est alors définie de façon unique par le fait que$s_{n+1}\geqslant s_{n}(s_{n}-1)+1$pour tout$n$et le réel$x$est rationnel si et seulement si$s_{n+1}=s_{n}(s_{n}-1)+1$à partir d'un certain rang.

Le développement du nombre 1 correspond bien sûr à la suite de Sylvester vue ci-dessus:$1=]2,3,7,43,1807,...[$.

Exemples irrationnels:

• ${\sqrt {2}}=1+]3,13,253,218201,...[$,voir la suiteA006487de l'OEIS;
• $e=2+]2,5,55,9999,3620211523,...[$,voir la suiteA006525de l'OEIS;
• $\pi =3+]8,61,5020,128541455,...[$,voir la suiteA001466de l'OEIS.

Si au lieu de prendre la plus grande fraction égyptienne dont la somme avec les précédentes eststrictementinférieure à$x$,on prend celle qui est inférieureou égaleà$x$,l'algorithme est identique si$x$est irrationnel, mais il s'arrête si$x$est rationnel et on obtient ledéveloppement égyptien glouton finide$x$.

Par exemple,${3 \over 16}={1 \over 6}+{1 \over 48}=]6,49,2353=49\times 48+1,5534257=2353\times 2352+1,...[$.

Le développement en série de Sylvester a des liens aveccelui de Engel.

Unicité des séries à croissance rapide ayant une limite rationnelle

Sylvester observa lui-même que la suite qui porte maintenant son nom semblait posséder la propriété unique d'avoir une croissance extrêmement rapide, tandis que la série somme de ses inverses convergeait vers un rationnel.

Plus précisément, il résulte des travaux deBadea (1993)que, si une suite d'entiers$(a_{n})$croît suffisamment pour qu'à partir d'un certain rang$n$:

$a_{n+1}\geqslant a_{n}(a_{n}-1)+1$

et si la série:$\sum {\frac {1}{a_{n}}}$converge vers unrationnel$A$,alors, pour tout$n$au-delà d'une certaine valeur, la suite peut être définie par la même relation de récurrence que la suite de Sylvester:

$a_{n+1}=a_{n}(a_{n}-1)+1$

En 1980,Erdősconjecturaque pour les résultats de ce type, l'inégalité conditionnant la croissance de la suite pouvait être remplacée par la condition plus faible:

$a_{n+1}\sim a_{n}^{2}$.

Applications

Boyer et al. (2005)utilisent les propriétés de la suite de Sylvester pour spécifier le grand nombre devariétés d'Einsteinsasakiennes(en)possédant latopologie différentielledesphèresde dimension impaire ou d'autressphères exotiques.Ils démontrent que le nombre demétriques riemanniennesdes variétés d'Einstein sasakiennes sur unesphère topologiquede dimension$2n-1$est au moins proportionnelle à$s_{n}$et croît donc selon une double exponentielle de$n$.

SelonGalambos et Woeginger (1995),Brown (1979)etLiang (1980)ont utilisé la suite de Sylvester pour construire un algorithme séquentiel minimisant leproblème de bin packing.Seiden et Woeginger (2005)ont aussi utilisé cette suite pour élaborer un algorithme minimisant leproblème de bin packingà deux dimensions.

Leproblème de Známconsiste à trouver un ensemble de nombres tel que chacun divise le produit de tous les autres plus 1, sans être égal à cette valeur. Si ce n'était cette dernière condition, les ensembles formés des premiers termes de la suite de Sylvester seraient des solutions à ce problème. Avec cette condition, les solutions constituent une suite dont la définition est similaire à celle de la suite de Sylvester. Les solutions du problème de Znám ont des applications dans la classification des singularités des surfacesBrenton et Hill (1988)et dans la théorie desautomates finis non déterministes(Domaratzkiet al.2005).

Curtiss (1922)présente une approximation de 1 par la somme de$k$fractions unitairescomme approximation inférieure du nombre de diviseurs de toutnombre parfaitetMiller (1919)utilise la même propriété pour minimiser la taille de certains groupes.

Notes et références

(en)Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé«Sylvester's sequence»(voir la liste des auteurs).

Voir aussi

Bibliographie

• (en)CatalinBadea,«A theorem on irrationality of infinite series and applications»,Acta Arithmetica,vol.63,‎1993,p.313-323
• (en)Catalin Badea, «On some criteria for irrationality for series of positive rationals: a survey»,1995
• (en)Charles P.Boyer,KrzysztofGalickietJános Kollár,«Einstein metrics on spheres»,Ann. of Math.,vol.162,n^o1,‎2005,p.557-580(DOI10.4007/annals.2005.162.557,lire en ligne)
• (en)LawrenceBrentonet RichardHill,«On the Diophantine equation 1=Σ1/n_i+ 1/Πn_iand a class of homologically trivial complex surface singularities»,Pacific J. Math.,vol.133,n^o1,‎1988,p.41-67(lire en ligne)
• (en)D. J.Brown,A lower bound for on-line one-dimensional bin packing algorithms; Version = Tech. Rep. R-864,Coordinated Science Lab.,Univ. Illinois, Urbana-Champaign,1979
• (en)D. R.Curtiss,«On Kellogg's diophantine problem»,Amer. Math. Monthly,vol.29,‎1922,p.380-387(DOI10.2307/2299023,JSTOR2299023)
• (en)MichaelDomaratzki,KeithEllul,Jeffrey Shallitet Ming-WeiWang,«Non-uniqueness and radius of cyclic unary NFAs»,International Journal of Foundations of Computer Science,vol.16,n^o5,‎2005,p.883-896(DOI10.1142/S0129054105003352,lire en ligne)
• (en)PaulErdősetRonald Graham,Old and New Problems and Results in Combinatorial Number Theory,Univ. de Genève,coll.« Monographies deL'Enseignement mathématique» (n^o28),1980
• (en)GáborGalambosetGerhard J.Woeginger,«On-line bin packing — A restricted survey»,Mathematical Methods of Operations Research,vol.42,n^o1,‎1995,p.25(DOI10.1007/BF01415672)
• (en)Solomon Golomb,«On certain nonlinear recurring sequences»,Amer. Math. Monthly,vol.70,n^o4,‎1963,p.403-405(DOI10.2307/2311857,JSTOR2311857)
• (en)Ronald Graham,Donald KnuthetOren Patashnik,Concrete Mathematics,Addison-Wesley,1989(ISBN0-201-55802-5),Exercise4.37
• (en)Rafe Jones, «The density of prime divisors in the arithmetic dynamics of quadratic polynomials»,2006(arXivmath/0612415v1)
• (en)Frank M. Liang, «A lower bound for on-line bin packing»,Information Processing Letters,vol.10,n^o2,‎1980,p.76-79(DOI10.1016/S0020-0190(80)90077-0)
• (en)G. A. Miller(en),«Groups possessing a small number of sets of conjugate operators»,Trans. Amer. Math. Soc.,vol.20,n^o3,‎1919,p.260-270(DOI10.2307/1988867,JSTOR1988867)
• (en)Martin Rosenman, «Problem 3536»,Amer. Math. Monthly,vol.40,n^o3,‎1933,p.180(JSTOR2301036)
• (en)H. E. Salzer,«The approximation of numbers as sums of reciprocals»,Amer. Math. Monthly,vol.54,n^o3,‎1947,p.135-142(DOI10.2307/2305906,JSTOR2305906)
• (en)Steven S.SeidenetGerhard J.Woeginger,«The two-dimensional cutting stock problem revisited»,Math. Prog.,vol.102,n^o3,‎2005,p.519-530(DOI10.1007/s10107-004-0548-1)
• (en)K. Soundararajan,«Approximating 1 from below usingnEgyptian fractions»,2005(arXivmath/0502247v1)
• (en)J. J. Sylvester,«On a point in the theory of vulgar fractions»,Amer. J. Math.,vol.3,n^o4,‎1880,p.332-335(DOI10.2307/2369261,JSTOR2369261)
• (en)IlanVardi,Computational Recreations in Mathematica,Redwood City (Calif.),Addison-Wesley,1991,286p.(ISBN0-201-52989-0),p.82-89

Articles connexes

• Constante de Cahen
• Développement en série de Engel
• Nombre pseudoparfait primaire

Liens externes

• (en)Irrationality of Quadratic Sums,MathPagesde K. S. Brown.
• (en)Eric W. Weisstein,«Sylvester's Sequence», surMathWorld