AtCoder&Luogu[ABC130E] Common Subsequence 题解

无名程序员就是我

已于 2022-12-31 11:03:49 修改

阅读量112

点赞数

分类专栏： C++ LUOGU Atcoder 文章标签：算法 c++ 动态规划

于 2022-12-31 10:57:52 首次发布

本文链接：https://blog.csdn.net/m0_66786097/article/details/128503566

版权

C++ 同时被 3 个专栏收录

20 篇文章 1 订阅

订阅专栏

LUOGU

8 篇文章 0 订阅

订阅专栏

Atcoder

7 篇文章 0 订阅

订阅专栏

这篇博客介绍了如何使用动态规划解决最长公共子序列（Longest Common Subsequence, LCS）问题。作者详细阐述了状态转移方程，并给出了C++代码实现。通过容斥原理解释了减去dp[i-1][j-1]的原因，确保不重复计数。代码中使用了模运算避免溢出。

摘要由CSDN通过智能技术生成

思路

题目要求我们求最长公共子序列，那么这就是一个 线性dp 了。

我们不妨确定一个设数，如 $dp_{i,j}$ 表示在序列 $a$ 中前 $i$ 个数，与数列 $b$ 中前 $j$ 个数所组成的最长序列的长度。那就有以下 $2$ 种情况。

若 $a_i=b_i$ ，那么： $dp_{i,j}=dp_{i-1,j}+dp_{i,j-1}$
若 $a_i\ne b_i$ ，那么： $dp_{i,j}=dp_{i-1,j}+dp_{i,j-1}-dp_{i-1,j-1}$

你可能要问：为什么要减去 $dp_{i-1,j-1}$ 呢？

回答：根据容斥原理，多的部分是要减掉的

代码如下

#include<iostream>
using namespace std;
int a[10010],b[10010],dp[10010][10010];
const int mod=1000000007;
int main(){
	int n,m;
	cin>>n>>m;
	for(int i=1;i<=n;i++) cin>>a[i];
	for(int i=1;i<=m;i++) cin>>b[i];
	dp[0][0]=1;
	for(int i=1;i<=n;i++)dp[i][0]=1;
	for(int j=1;j<=m;j++)dp[0][j]=1; 
	for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++){
			if(a[i]==b[j]){
				dp[i][j]=(dp[i-1][j]+dp[i][j-1])%mod;
			}
			else{
				dp[i][j]=(dp[i-1][j]+dp[i][j-1]-dp[i-1][j-1])%mod;
			}
		}
	}
	cout<<(dp[n][m]%mod+mod)%mod;
	return 0;
}