最长公共子序列LCS和最长回文子序列的动态规划算法

最新推荐文章于 2024-02-24 10:30:00 发布

NicoVon

最新推荐文章于 2024-02-24 10:30:00 发布

阅读量308

点赞数

分类专栏：数据结构和算法文章标签： php 算法 LCS

本文链接：https://blog.csdn.net/ray0354315/article/details/53536343

版权

数据结构和算法专栏收录该内容

14 篇文章 0 订阅

订阅专栏

c[i][j] 用来表示 x[i] 和y[j] 的LCS 长度

对c[i][j]，

若i = 0或j = 0，则c[i][j] = 0，两个原序列有一个为空，则没有LCS；

若i,j > 0 且x[i] = y[j]，则c[i][j] = c[i-1][j-1] + 1；

若i,j > 0 且x[i] != y[j]，则c[i][j] = max(c[i,j-1], c[i-1,j])

<?php
$x = array('A','B','C','B','D','A','B');
$y = array('B','D','C','A','B','A');
$c = lcsLength($x, $y);
//var_dump($c);
printLcs($c, $x, $y, count($x), count($y));

function lcsLength($x, $y){
    $m = count($x);
    $n = count($y);
    $c = array();
    
    for ($i = 0; $i <= $m; $i++){
        $c[$i][0] = 0;
    }
    for ($j = 0; $j <= $n; $j++){
        $c[0][$j] = 0;
    }
    
    $i = 1;
    while ($i <= $m) {
    	$j = 1;
        while ($j <= $n) {
            if($x[$i-1] == $y[$j-1]){
                $c[$i][$j] = $c[$i - 1][$j - 1] + 1;
            }else if ($c[$i - 1][$j] >= $c[$i][$j - 1]){
            	$c[$i][$j] = $c[$i - 1][$j];
			}else{
				$c[$i][$j] = $c[$i][$j - 1];
			}
			$j++;
        }
        $i++;
    }
    return $c;
}

function printLcs($c, $x, $y, $i, $j){
	if($i == 0 || $j == 0)
		return;
	if($x[$i - 1] == $y[$j - 1]){
		printLcs($c, $x, $y, $i-1, $j-1);
		echo $x[$i - 1];
	}elseif ($c[$i - 1][$j] >= $c[$i][$j - 1]){
		printLcs($c, $x, $y, $i-1, $j);
	}else{
		printLcs($c, $x, $y, $i, $j-1);
	}
}

拓展：最长回文子序列的动态规划算法

实际上求最长回文子序列即是求字符串与其倒序的LCS，状态转移方程式为：

对给定x[0...n-1] 有dp[i][j]表示其最长回文子序列的值，i指向字符串头，j指向尾

若i=j，则dp[i][j] = 1;

若x[i] = x[j]，dp[i][j] = dp[i+1][j-1] + 2;

若x[i] != x[j]，dp[i][j] = max(dp[i-1][j], dp[i][j-1]);