最长公共子序列(poj1458)

本题来自poj1458,题目大意:给定两个字符串X,Y,求X,Y的最长公共子序列,不要求子序列连续,但子序列在X,Y中的位置是严格单调递增的。

例如 abcfbc与abfcab的最长公共子序列就是abcd,其长度为4;programming和contest的最长公共子序列就是on,其长度为2;abcd和mnp的最长公共子序列为空

可以采用动态规划的方法来解该题,计算字符串A、B的最长公共子序列

  令L(i,j)表示A中前i个字符和B中前j个字符的最长公共子序列
  则L(i+1,j+1)等于多少?
  1.当A[i+1]=B[j+1]时,L(i+1,j+1) = L(i,j)+1
  2.如果A[i+1] != B[j+1],则L(i+1,j+1) = max(L(i+1,j), L(i,j+1))

在上边界和右边界上的值为零

#include <stdio.h>
#include <string.h>

int L[1000][1000];
/*计算字符串A、B的最长公共子序列,采用动态规划的解法
  L(i,j)表示A中前i个字符和B中前j个字符的最长公共子序列
  则L(i+1,j+1)等于多少?
  1.当A[i+1]=B[j+1]时,L(i+1,j+1) = L(i,j)+1
  2.如果A[i+1] != B[j+1],则L(i+1,j+1) = max(L(i+1,j), L(i,j+1))
*/
int TheCommonSubsequence(char *strA, int lengthA, char *strB, int lengthB)
{
	int i, j;

	for (i = 0; i < lengthA; i++)
	{
		L[0][i] = 0;
		L[i][0] = 0;
	}

	for (i = 1; i < lengthA; i++)
	{
		for (j = 1; j < lengthB; j++)
		{
			if (strA[i] == strB[j])
				L[i][j] = L[i - 1][j - 1] + 1;
			else
			{
				L[i][j] = (L[i - 1][j] > L[i][j - 1]) ? L[i - 1][j] : L[i][j - 1];
			}
		}
	}
	return L[lengthA - 1][lengthB - 1];
}

int main()
{
	char strA[1000], strB[1000];
	
	while (EOF != scanf("%s %s", strA + 1, strB + 1))
	{
		printf("%d\n", TheCommonSubsequence(strA, strlen(strA), strB, strlen(strB)));
	}
	return 0;
}


  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
以下是Java解决POJ3233—矩阵幂序列问题的代码和解释: ```java import java.util.Scanner; public class Main { static int n, k, m; static int[][] A, E; public static void main(String[] args) { Scanner sc = new Scanner(System.in); n = sc.nextInt(); k = sc.nextInt(); m = sc.nextInt(); A = new int[n][n]; E = new int[n][n]; for (int i = 0; i < n; i++) { for (int j = 0; j < n; j++) { A[i][j] = sc.nextInt() % m; E[i][j] = (i == j) ? 1 : 0; } } int[][] res = matrixPow(A, k); int[][] ans = matrixAdd(res, E); printMatrix(ans); } // 矩阵乘法 public static int[][] matrixMul(int[][] a, int[][] b) { int[][] c = new int[n][n]; for (int i = 0; i < n; i++) { for (int j = 0; j < n; j++) { for (int k = 0; k < n; k++) { c[i][j] = (c[i][j] + a[i][k] * b[k][j]) % m; } } } return c; } // 矩阵快速幂 public static int[][] matrixPow(int[][] a, int b) { int[][] res = E; while (b > 0) { if ((b & 1) == 1) { res = matrixMul(res, a); } a = matrixMul(a, a); b >>= 1; } return res; } // 矩阵加法 public static int[][] matrixAdd(int[][] a, int[][] b) { int[][] c = new int[n][n]; for (int i = 0; i < n; i++) { for (int j = 0; j < n; j++) { c[i][j] = (a[i][j] + b[i][j]) % m; } } return c; } // 输出矩阵 public static void printMatrix(int[][] a) { for (int i = 0; i < n; i++) { for (int j = 0; j < n; j++) { System.out.print(a[i][j] + " "); } System.out.println(); } } } ``` 解释: 1. 首先读入输入的n、k、m和矩阵A,同时初始化单位矩阵E。 2. 然后调用matrixPow函数求出A的k次幂矩阵res。 3. 最后将res和E相加得到结果ans,并输出。 4. matrixMul函数实现矩阵乘法,matrixPow函数实现矩阵快速幂,matrixAdd函数实现矩阵加法,printMatrix函数实现输出矩阵。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值