Plus longue sous-séquence commune

Il est possible de ramener le problème de recherche de plus longue sous séquence commune (PLSC) entre deux chaînes données à une recherche entre deux chaînes de taille inférieure grâce au théorème suivant (où ${\displaystyle X_{l}}$ désigne les ${\displaystyle l}$ premiers caractères de la séquence ${\displaystyle X}$)[2]:

Les trois cas ${\displaystyle \left\langle x_{m}=y_{n}\right\rangle }$, ${\displaystyle \left\langle x_{m}\neq y_{n}~{\mathsf {ET}}~z_{k}\neq x_{m}\right\rangle }$ et ${\displaystyle \left\langle x_{m}\neq y_{n}~{\mathsf {ET}}~z_{k}\neq y_{n}\right\rangle }$ sont exhaustifs, ce qui permet bien de se ramener à un problème de taille inférieure.

On crée un tableau à deux dimensions ${\displaystyle c[1\cdot \cdot m][1\cdot \cdot n]}$ dans lequel chaque case ${\displaystyle c[i][j]}$ est destiné à contenir la longueur des PLSCs entre ${\displaystyle X_{i}}$ et ${\displaystyle Y_{j}}$. On peut alors calculer de proche en proche ${\displaystyle c[i][j]}$ pour chaque couple d'indice ${\displaystyle i}$ et ${\displaystyle j}$. Du théorème précédent découle en effet la formule[2]:

${\displaystyle c[i][j]={\begin{cases}0&{\text{ si }}i=0{\text{ ou }}j=0,\\c[i-1][j-1]+1&{\text{ si }}i,j>0{\text{ et }}x_{i}=y_{j},\\{\mathsf {max}}(c[i][j-1],c[i-1][j])&{\text{ si }}i,j>0{\text{ et }}x_{i}\neq y_{j}.\end{cases}}}$

Le calcul du contenu des cases de ${\displaystyle c}$ peut être effectué avec une complexité ${\displaystyle {\mathcal {O}}(mn)}$, car le contenu de chaque case est calculable à partir des cases précédente en ${\displaystyle {\mathcal {O}}(1)}$[2].

La formule précédente permet de calculer de proche en proche les cases de ${\displaystyle c}$. On peut reconstituer une plus longue sous-séquence commune grâce à lui.

Pour cela on effectue un parcours depuis ${\displaystyle c[m][n]}$ suivant la règle suivante

Depuis une case ${\displaystyle c[i][j]}$ de valeur ${\displaystyle \alpha }$:

• Si ${\displaystyle x_{i}=y_{j}}$, on passe à la case ${\displaystyle c[i-1][j-1]}$ de valeur ${\displaystyle \alpha -1}$ et on ajoute ce caractère (${\displaystyle x_{i}=y_{j}}$) au début de la PLSC en construction.

• Si ${\displaystyle x_{i}\neq y_{j}}$,
  • Si ${\displaystyle c[i][j-1]=c[i-1][j]=\alpha }$, on passe indifféremment à la case ${\displaystyle c[i-1][j]}$ ou ${\displaystyle c[i][j-1]}$.
  • Si ${\displaystyle c[i][j-1]=\alpha >c[i-1][j]}$, on passe à la case ${\displaystyle c[i][j-1]}$
  • Si ${\displaystyle c[i-1][j]=\alpha >c[i][j-1]}$, on passe à la case ${\displaystyle c[i-1][j]}$

Un exemple de parcours est donné par le tableau suivant, grâce auquel on déduit que MJAU est une plus longue sous-séquence commune à MZJAWXU et XMJYAUZ :

Le calcul du contenu des cases de ${\displaystyle c}$ peut être effectué avec une complexité ${\displaystyle {\mathcal {O}}(mn)}$, car le contenu de chaque case est calculable à partir des cases précédente en ${\displaystyle {\mathcal {O}}(1)}$[2].

Une fois ${\displaystyle c}$ connu, l'obtention d'une PLSC a une complexité ${\displaystyle {\mathcal {O}}(m+n)}$[2].