Fonction de compte des nombres premiers - Définition

Source: Wikipédia sous licence CC-BY-SA 3.0.
La liste des auteurs de cet article est disponible ici.

Introduction

En mathématiques, la fonction de compte des nombres premiers est la fonction comptant le nombre de nombres premiers inférieurs ou égaux à un nombre réel x. Elle est notée \pi(x)\, (à ne pas confondre avec la constante π).

Les 60 premières valeurs de π(n)

Historique

Depuis Euclide, il est connu qu'il existe des nombres premiers en quantité infinie. Pour affiner la connaissance de ces nombres, la théorie des nombres s'est attelée à en déterminer le taux de croissance. À la fin du XVIIIe siècle, Gauss et Legendre ont conjecturé que cette quantité était proche de  x/\operatorname{ln}(x) .

De manière équivalente, on peut l'écrire

\lim_{x\rightarrow +\infty}\frac{\pi(x)}{x/\operatorname{ln}(x)}=1

Cette affirmation constitue le théorème des nombres premiers, prouvé indépendamment par Hadamard et de La Vallée Poussin, en 1896, grâce à la fonction zêta de Riemann. Une assertion équivalente est:

\lim_{x\rightarrow +\infty}\pi(x) / \operatorname{Li}(x)=1 ,

\operatorname{Li}\, est la fonction Logarithme intégral.

Des estimateurs plus précis de \pi(x)\, sont aujourd'hui connus. Par exemple:

\pi(x) = \operatorname{li}(x) + O \left( x \exp \left( -\frac{\sqrt{\ln(x)}}{15} \right) \right) ,

Des preuves du théorème des nombres premiers n'utilisant pas l'analyse complexe furent proposée en 1948 par Atle Selberg et Paul Erdős.

Autres fonctions de compte des nombres premiers

D'autres fonctions de compte des nombres premiers sont aussi utilisées car elles sont plus pratiques pour travailler. Une d'elles est la fonction de compte des nombres premiers de Riemann, notée \Pi(x)\, ou J(x)\, . Celle-ci possède des sauts de 1/n pour les puissances de nombres premiers pn, qui prennent une valeur à mi-chemin entre les deux côtés des discontinuités. Elle peut être aussi définie comme une transformation de Mellin inverse. Formellement, nous pouvons définir J par

J(x) = \frac12 \bigg(\sum_{p^n < x} \frac1n\ + \sum_{p^n \le x} \frac1n\bigg)

p est un nombre premier.

Nous pouvons aussi écrire

J(x) = \sum_{n=1}^\infty \frac1n \pi(x^{1/n}) et
 \ln( \zeta (s))=s\int_{0}^{\infty}dxJ(x) x^{-s+1}

et en connaissant la relation entre la fonction log de la fonction Riemann et la fonction de von Mangoldt  \Lambda\, , et la formule de Perron, nous avons :

 J(x)= \sum_{2}^{x} \frac{ \Lambda (n)}{ln(n)}

excepté où se trouvent les discontinuités aux puissances de nombres premiers, et ainsi, \pi(x)\, peut être retrouvé à partir de J par l'inversion de Möbius.

La fonction de Tchebychev pèse les nombres premiers ou les puissances de nombres premiers pn par ln p:

\theta(x)=\sum_{p\le x}\ln p,
\psi(x) = \frac12 \bigg(\sum_{p^n < x} \ln p\ + \sum_{p^n \le x} \ln p\bigg).

Excepté aux discontinuités des puissances de nombres premiers, nous avons

\psi(x)=\sum_{n=1}^\infty\theta(x^{1/n})=\sum_{n\le x}\Lambda(n),

\Lambda(n)\, est la fonction de von Mangoldt.

Algorithmes d'évaluation de π(x)

Une façon simple de calculer \pi(x)\, , si x\, est un nombre petit, est d'utiliser le crible d'Ératosthène de manière à trouver tous les premiers inférieurs à x\, et ensuite de les compter.

Une manière plus élaborée pour trouver \pi(x)\, a été inventée par Legendre: étant donné x\, , si p_1\, p_2\, , …,  p_k\, sont des nombres premiers distincts, alors le nombre d'entiers inférieurs ou égaux à x qui ne sont divisibles par aucun p_i\, est

\lfloor x\rfloor - \sum_{i}\left\lfloor\frac{x}{p_i}\right\rfloor + \sum_{i<j}\left\lfloor\frac{x}{p_ip_j}\right\rfloor - \sum_{i<j<k}\left\lfloor\frac{x}{p_ip_jp_k}\right\rfloor + \cdots , (où \lfloor\cdot\rfloor désigne la fonction Partie entière).

Ce nombre est donc égal à : \pi(x)-\pi\left(\sqrt{x}\right)+1\, où les nombres p_1\, p_2\, , …,  p_k\, sont les nombres premiers inférieurs ou égaux à la racine carrée de x\, .

Dans une série d'articles publiés entre 1870 et 1885, Ernst Meissel décrivit (et utilisa) une manière combinatoire pratique pour évaluer \pi(x)\, . Soit p_1\, p_2\, , …,  p_n\, les premiers n\, nombres premiers et notons par \Phi(m,n)\, le nombre de nombres naturels inférieurs à m\, qui ne sont divisibles par aucun p_i\, . Alors

\Phi(m,n)=\Phi(m,n-1)-\Phi\left(\left[\frac{m}{p_n}\right],n-1\right),\,

Soit un nombre naturel donné m\, , si n=\pi\left(\sqrt[3]{m}\right)\, et si \mu=\pi\left(\sqrt{m}\right)-n\, , alors

\pi(m)=\Phi(m,n)+n(\mu+1)+\frac{\mu^2-\mu}{2}-1-\sum_{k=1}^\mu\pi\left(\frac{m}{p_{n+k}}\right).\,

En utilisant cette approche, Meissel calcula \pi(x)\, , pour x\, égal à 5x105, 106, 107, et 108.

En 1959, Derrick Lehmer étendit et simplifia la méthode de Meissel. Définissons, pour un réel m et pour des nombres naturels n, et k, Pk(m,n) comme le nombre de nombres inférieurs à m avec exactement k facteurs premiers, tous plus grands que pn. De plus, fixons P_0(m,n)=1\, . Alors

\Phi(m,n)=\sum_{k=0}^{+\infty}P_k(m,n),\,

où la somme actuelle possède seulement de manière finie plusieurs termes différents de zéro. Soit y désignant un entier tel que \sqrt[3]{m}\le y\le\sqrt{m}\, , et fixons n=\pi(y)\, . Alors P_1(m,n)=\pi(m)-n\, et P_k(m,n)=0\, lorsque k ≥ 3. Par conséquent

\pi(m)=\Phi(m,n)+n-1-P_2(m,n)\, .

Le calcul de P_2(m,n)\, peut être obtenu de cette manière :

P_2(m,n)=\sum_{y<p\le\sqrt{m}}\left(\pi\left(\frac mp\right)-\pi(p)+1\right).\,

D'un autre côté, le calcul de Φ(m,n) peut être fait en utilisant les règles suivantes :

  1. \Phi(m,0)=\lfloor m\rfloor;\,
  2. \Phi(m,b)=\Phi(m,b-1)-\Phi\left(\frac m{p_b},b-1\right).\,

En utilisant cette méthode et un IBM 701, Lehmer a été capable de calculer \pi\left(10^{10}\right) .

Le mathématicien chinois Hwang Cheng, dans une conférence sur les fonctions de nombres premiers à l'université de Bordeaux utilisa les identités suivantes :

e^{(a-1)\Theta}f(x)=f(ax)\,
J(x)=\sum_{n=1}^{\infty}\frac{\pi(x^{1/n})}{n}\,

et en notant x=e^t\, , en appliquant la transformée de Laplace aux deux côtés et en appliquant une somme géométrique sur e^{n\Theta}\, , on obtient l'expression :

 \frac{1}{2{\pi}i}\int_{c-i\infty}^{c+i\infty}g(s)t^{s}=\pi(t)
 (1-e^{\Theta})\frac{ln\zeta(s)}{s}=g(s)
 \Theta(s)=s\frac{d}{ds}.
Page générée en 0.231 seconde(s) - site hébergé chez Contabo
Ce site fait l'objet d'une déclaration à la CNIL sous le numéro de dossier 1037632
A propos - Informations légales | Partenaire: HD-Numérique
Version anglaise | Version allemande | Version espagnole | Version portugaise