线性DP-入门篇

目录

数字三角形：

最长上升子序列：

魔族密码： 

编辑距离： 

---

线性动态规划的主要特点是状态转移的推导是按照问题规模 从小到大依次推导，较大规模的问题的解依赖较小规模的问题的解。

数字三角形：

[USACO1.5][IOI1994]数字三角形 Number Triangles - 洛谷https://www.luogu.com.cn/problem/P1216我们来看一道经典的问题数字三角形问题，这个问题应该是每一个学DP的人都会遇到的问题。

闫式DP分析法：

 算法分析：

我们用两维来表示状态，将f[i][j]分为从左上方走下来还是从右上方走下来,然后我们取一个max 最后加上该点的权值即可。

注意：

如果我们遇到f[i][j]的值存在负数的情况的时候，边界上的点（不存在从左上方的点和不存在右上方的点）在计算的时候可能出错，我们需要对f[i][j]进行初始化为负无穷即可。

#include
#include
#include
#include

using namespace std;

const int N = 1010,INF = 1e8;

int n,res;

int a[N][N],f[N][N];

int main()
{
	scanf("%d",&n);
	
	for(int i =1 ;i <= n ;i ++)
		for(int j = 1;j <= i ; j++)
			scanf("%d",&a[i][j]); 
	
	
	for(int i =1 ;i<= n ;i ++)
	{
		for(int j = 1 ;j <= i ; j++)
		{
			f[i][j] = max(f[i-1][j],f[i-1][j-1]) + a[i][j] ;
		}
	}
	
	
	for(int i =1 ;i<= n ;i ++)
	{
		res = max(res,f[n][i]);
	}
	
	printf("%d",res);
	
	return 0;
}

---

最长上升子序列：

最长上升子序列 - 洛谷https://www.luogu.com.cn/problem/B3637最长上升子序列也是一个经典的DP问题了，我们还是从闫式DP分析法开始分析。

闫式DP分析法：

 算法分析：

对于 f[i]我们将其分为0 ~ i-1类，对于a[0] ~a[i-1](0 <= j <= i-1)，当存在a[j] < a[i] 时，说明以j结尾的序列可以接上 i ，进行状态转移。  

#include
#include
#include
#include

using namespace std;

const int N = 1000010,INF = 1e8;

int n,res;
int a[N],f[N];

int main()
{
	scanf("%d",&n);
	for(int i = 1 ;i <= n ;i ++) cin >> a[i];
	
	for(int i = 1 ;i <= n ;i ++)
	{
		f[i] = 1; 
		for(int j = 1 ;j<= i-1 ;j ++)
		{
			if(a[j] < a[i])
			{
				f[i] = max(f[i],f[j] + 1);	
			} 
		}
		
		res = max(res,f[i]);
	}
	
	
	printf("%d",res);
	
	return 0;
} 

魔族密码： 

魔族密码 - 洛谷https://www.luogu.com.cn/problem/P1481

算法分析：

魔族密码是最长上升子序列的一个变种题目，换汤不换药，在最长上升子序列中我们的一个判断条件是是否在0 ~ i - 1中存在a[j] < a[i]，进行状态转移；我们这个题目和最长上升子序列的转移方程一模一样，只不过条件判断变成了0 ~ i - 1中的字符串是否是当前字符串的前缀。

#include
#include
#include
#include

using namespace std;

const int N = 2010,INF = 1e8;

int n,f[N],res;

string str[N];

int main()
{
	cin >> n;
	
	for(int i = 1 ; i <= n ;i ++) cin >> str[i];
	
	for(int i = 1 ;i <= n ;i ++)
	{
		f[i] = 1;
		for(int j = 1;j <= i - 1 ; j ++)
		{
			int len = str[j].length();
			if(str[i].substr(0,len) == str[j])
			{
				f[i] = max(f[i],f[j] + 1);
			} 
		}
		res = max(res,f[i]);
	}
	
	cout << res<

编辑距离： 

编辑距离 - 洛谷https://www.luogu.com.cn/problem/P2758

算法分析：

根据闫式DP分析法，我们考虑 i , j 删除、添加、修改；

当考虑删除第 i 个字符与前 j 个字符匹配时，；

当考虑添加第 i 个字符和前 j - 1 个字符匹配，；添加的第 i 个字符其实就是 b[j];

当考虑修改第 i 个字符和前 j 个字符，有两种情况 a[i] 和 b[j] 已经相等，考虑前i - 1 和前j - 1 个字符，如果不相等还需要一次操作将 a[i] 修改成 b[j] 。；

转移方程：

#include
#include
#include
#include

using namespace std;

const int N = 2010;


char a[N],b[N];

int f[N][N],lena,lenb;

int main()
{
	scanf("%s%s",a + 1 , b + 1);
	
	lena = strlen(a + 1);
	lenb = strlen(b + 1);
	
	for(int i = 0; i <= lena;i ++)f[i][0] = i;
	for(int j = 0; j <= lenb;j ++)f[0][j] = j;
	
	for(int i =1; i<= lena; i++)
	{
		for(int j =1; j<= lenb ;j ++)
		{
			f[i][j] = min(f[i-1][j] + 1, f[i][j-1] + 1);
			f[i][j] = min(f[i][j],f[i-1][j-1] + (a[i]!=b[j]));
		}
	}
	
	printf("%d",f[lena][lenb]);
	
	return 0;
}