Tri de Shell

Tri de Shell barres de couleur de l'algorithme

Le tri de Shell ou Modèle:Langue en anglais est un algorithme de tri. C'est une amélioration notable du tri par insertion au niveau de la vitesse d'exécution, mais ce tri n'est pas stable. Le principe de l'algorithme est simple mais l'étude du temps d'exécution est très complexe, et plusieurs problèmes sont toujours ouverts à ce sujet.

Le nom vient de son inventeur Modèle:Lien (1924-2015) qui publia l'algorithme dans le numéro de Modèle:Date de Communications of the ACM.

Principe

Amélioration du tri par insertion

Le tri de Shell est une amélioration du tri par insertion en observant deux choses :

• Le tri par insertion est efficace si la liste est à peu près triée (1) ;
• Le tri par insertion est inefficace en moyenne car il ne déplace les valeurs que d'une position par instruction (2).

Le tri de Shell trie chaque liste d'éléments séparés de n positions chacun avec le tri par insertion. L'algorithme effectue plusieurs fois cette opération en diminuant n jusqu'à n=1 ce qui équivaut à trier tous les éléments ensemble.

Le fait de commencer avec des éléments espacés permet de pallier l'inconvénient (2), tandis que lorsque l'on fait à la fin avec un espacement de 1, ce qui est en fait un tri par insertion ordinaire, on tire parti de l'avantage (1).

Code Python

gaps = [701, 301, 132, 57, 23, 10, 4, 1]

def shell_sort(tab):
  n = len(tab)
  for m in gaps:
    for r in range(m):
      # tri par insertion des positions de la forme k * m + r
      for i in range (r + m, n, m):
        j = i
        x = tab[i]
        while j > r and tab[j-m] > x:
          tab[j] = tab[j-m]
          j = j - m
        tab[j] = x

Gap ou espacement entre les éléments

Les premiers espacements optimaux (empiriquement trouvés) sont les suivants : 1, 4, 10, 23, 57, 132, 301, 701.

On remarque que le facteur entre ces valeurs, exception faite des deux premières, est d'environ 2,3. On peut effectivement prolonger cette liste avec ce facteur si les dimensions du tableau dépassent environ 1600 éléments. Par exemple pour étendre la liste gaps jusqu'à la taille nécessaire:

gap = gaps[0]
while gap<length(liste):
  gap = round(gap*2.3);
  gaps = [gap] + gaps

Des espacements de la forme $2^p{3}^q$, dans l'ordre croissant, garantissent quant à eux la meilleure complexité théorique prouvée aujourd'hui, $O(n{\log }^{2}n)$.

Complexité

Sur des tableaux de moins d'une centaine d'éléments, ce tri est aussi rapide qu'un tri rapide simple. Mais plutôt que d'être en compétition avec l'algorithme quicksort, il peut être utilisé pour son optimisation quand les sous-listes à traiter deviennent petites.

Le choix de la suite des espacements entre les éléments qu'on trie à chaque étape (gap) est très important. Il peut faire varier la complexité dans le pire cas de $O\left(n^2\right)$ à $O\left(n{\log }^{2}n\right)$. Il est également possible que la complexité en moyenne puisse être $O(n\log n)$ avec un bon choix d'espacements (problème ouvert).

Des bornes inférieures ont aussi été publiées, on peut citer une borne de $\mathrm{\Omega }\left(n{\left(\frac{\log n}{\log \log n}\right)}^2\right)$ sur la complexité dans le pire cas quels que soient les espacements.

Le tableau suivant compare les gaps publiés jusqu'à aujourd'hui :

OEIS	Terme général (k ≥ 1)	Gaps ou espacements concrets	Complexité dans le pire des cas	Auteur et année de publication
	$\lfloor \frac{N}{2^k}\rfloor$	$\lfloor \frac{N}{2}\rfloor ,\lfloor \frac{N}{4}\rfloor ,\dots ,1$	$\mathrm{\Theta }\left(N^2\right)$ [i.e quand N = 2p]	Modèle:Lien, 1959
	$2\lfloor \frac{N}{2^{k+1}}\rfloor +1$	$2\lfloor \frac{N}{4}\rfloor +1,\dots ,3,1$	$\mathrm{\Theta }\left(N^{\frac{3}{2}}\right)$	Frank & Lazarus, 1960
A168604	$2^k-1$	$1,3,7,15,31,63,\dots$	$\mathrm{\Theta }\left(N^{\frac{3}{2}}\right)$	Modèle:Lien, 1963
A083318	$2^k+1$, préfixé avec 1	$1,3,5,9,17,33,65,\dots$	$\mathrm{\Theta }\left(N^{\frac{3}{2}}\right)$	Papernov & Stasevich, 1965
A003586	Nombres successifs de la forme $2^p{3}^q$ (entier friable 3-lisse)	$1,2,3,4,6,8,9,12,\dots$	$\mathrm{\Theta }\left(N{\log }^{2}N\right)$	Modèle:Lien, 1971
A003462	$\frac{3^k-1}{2}$, plus petit que $\lceil \frac{N}{3}\rceil$	$1,4,13,40,121,\dots$	$\mathrm{\Theta }\left(N^{\frac{3}{2}}\right)$	Modèle:Lien, 1971
A036569	$\begin{array}{cc}& \prod _I{a}_q,\text{où}\\ a_q=& \min \{n\in \mathbb{N}:n\ge {\left(\frac{5}{2}\right)}^{q+1},\mathrm{\forall }p:0\le p<q\Rightarrow pgcd(a_p,n)=1\}\\ I=& \{0\le q<r\mid q\ne \frac{1}{2}(r^2+r)-k\}\\ r=& \lfloor \sqrt{2k+\sqrt{2k}}\rfloor \end{array}$	$1,3,7,21,48,112,\dots$	$O\left(N^{1+\sqrt{\frac{8\ln (5/2)}{\ln (N)}}}\right)$	Incerpi & Robert Sedgewick, 1985, Knuth
A036562	$4^k+3\cdot 2^{k-1}+1$, préfixé par 1	$1,8,23,77,281,\dots$	$O\left(N^{\frac{4}{3}}\right)$	Sedgewick, 1986
A033622	$\{\begin{array}{cc}9(2^k-2^{\frac{k}{2}})+1& k\text{ pair},\\ 8\cdot 2^k-6\cdot 2^{(k+1)/2}+1& k\text{ impair}\end{array}$	$1,5,19,41,109,\dots$	$O\left(N^{\frac{4}{3}}\right)$	Sedgewick, 1986
	$h_k=\max \{\lfloor \frac{5h_{k-1}}{11}\rfloor ,1\},h_0=N$	$\lfloor \frac{5N}{11}\rfloor ,\lfloor \frac{5}{11}\lfloor \frac{5N}{11}\rfloor \rfloor ,\dots ,1$	Inconnue	Modèle:Lien & Modèle:Lien, 1991
A108870	$\lceil \frac{1}{5}(9\cdot {\left(\frac{9}{4}\right)}^{k-1}-4)\rceil$	$1,4,9,20,46,103,\dots$	Inconnue	Tokuda, 1992
A102549	Inconnue (trouvé expérimentalement)	$1,4,10,23,57,132,301,701$	Inconnue	Ciura, 2001

Liens externes

• Shell Sort en C, Java, Pascal, C#
• Modèle:En Illustration dynamique du tri shell