最长公共子序列及其输出

给出两个字符串,找到最长公共子序列(LCS),返回LCS的长度。
说明
最长公共子序列的定义:
最长公共子序列问题是在一组序列(通常2个)中找到最长公共子序列(注意:不同于子串,LCS不需要是连续的子串)。该问题是典型的计算机科学问题,是文件差异比较程序的基础,在生物信息学中也有所应用。
https://en.wikipedia.org/wiki/Longest_common_subsequence_problem
样例
给出"ABCD" 和 "EDCA",这个LCS是 "A" (或 D或C),返回1


给出 "ABCD" 和 "EACB",这个LCS是"AC"返回 2

package com.lianxi.dp;

import java.util.Scanner;

//最长公共子序列
public class LCS {

	
	 public static void backtrack(int i,int j,int [][]b, String s)
	 {
		 if(i==0||j==0)
		  return;
		 else if(b[i][j]==0)
		 {
			 backtrack(i-1,j-1,b,s);
			 System.out.print(s.charAt(i-1));
		 }
		 else if(b[i][j]==-1)
		 {
			 backtrack(i-1,j,b,s);
		 }
		 else
		 {
			 backtrack(i,j-1,b,s);
		 }
		 
		 
	 }

	public static void main(String[] args) {
	  Scanner sc=new Scanner(System.in);
	  String s1=sc.nextLine();
	  String s2=sc.nextLine();
	  int d1=s1.length();
	  int d2=s2.length();
	  int [][]c=new int[s1.length()+1][s2.length()+1];
	  int [][]b=new int[s1.length()+1][s2.length()+1];//记录获得c[i][j]的路径分别为0,-1,1,便于回溯公共子序列
	  for(int s=0;s<s1.length()+1;s++)
		  c[s][0]=0;
	  for(int t=0;t<s2.length()+1;t++)
		  c[0][t]=0;
	  for(int i=1;i<s1.length()+1;i++)
		  for(int j=1;j<s2.length()+1;j++)
		  {
			  if(s1.charAt(i-1)==s2.charAt(j-1))
			  {
				  c[i][j]=c[i-1][j-1]+1;
				  b[i][j]=0;
			  }
			  else if(c[i-1][j]>c[i][j-1])
			  {    c[i][j]=c[i-1][j];
			       b[i][j]=-1;
			  }
			  else
			  {    c[i][j]=c[i][j-1];
			       b[i][j]=1;
			  }
			  
		  }
	     backtrack(d1,d2,b,s1);
	   System.out.println();
	  System.out.println(c[s1.length()][s2.length()]);

	}

}



以下是C语言代码实现: ``` #include <stdio.h> #include <string.h> #define MAX_LEN 100 int max(int a, int b) { return (a > b) ? a : b; } int lcs(char *str1, char *str2, int m, int n, char *result) { // 创建动态规划表 int lcs_table[m+1][n+1]; int i, j; for (i = 0; i <= m; i++) { for (j = 0; j <= n; j++) { if (i == 0 || j == 0) lcs_table[i][j] = 0; else if (str1[i-1] == str2[j-1]) lcs_table[i][j] = lcs_table[i-1][j-1] + 1; else lcs_table[i][j] = max(lcs_table[i-1][j], lcs_table[i][j-1]); } } // 反向追踪lcs i = m, j = n; int index = lcs_table[m][n]; result[index] = '\0'; while (i > 0 && j > 0) { if (str1[i-1] == str2[j-1]) { result[index-1] = str1[i-1]; i--; j--; index--; } else if (lcs_table[i-1][j] > lcs_table[i][j-1]) i--; else j--; } return lcs_table[m][n]; } int main() { char str1[MAX_LEN], str2[MAX_LEN], result[MAX_LEN]; printf("Enter string1: "); scanf("%s", str1); printf("Enter string2: "); scanf("%s", str2); int m = strlen(str1); int n = strlen(str2); int lcs_len = lcs(str1, str2, m, n, result); printf("Length of LCS: %d\n", lcs_len); printf("LCS: %s\n", result); return 0; } ``` 其中,`lcs`函数用于计算最长公共子序列,`max`函数用于返回两个整数中的较大值。首先创建一个大小为`(m+1) x (n+1)`的动态规划表`lcs_table`,其中`lcs_table[i][j]`表示`str1[0:i-1]`和`str2[0:j-1]`的最长公共子序列的长度。然后,通过填充表中的值,可以找到最长公共子序列的长度。最后,通过反向跟踪表,可以找到最长公共子序列本身。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值