线性动态规划问题

一、LIS（最长上升子序列）问题

题目描述

给定一个长度为 N 的数列 (w[N]) ，求数值严格单调递增的子序列的长度最长是多少。

输入格式

第一行包含整数 N 。第二行包含 N 个整数，表示完整序列。

输出格式

输出一个整数，表示最大长度。数据范围 1 ≤ N ≤ 1000 10 9 ≤ w [ i ] ≤ 10 9

Sample Input

3 1 2 1 8 5 6

Sample Output

4

题目分析

这道题是最基础的线性dp问题。

设置一个数组 f[i]表示以第 i 个数结尾的最长上升子序列长度。f[i]的大小等于在它前面的所有上升子序列中（并且 a[i]要大于前面的结尾数）的最大数+1。也就是 if (a[i] > a[j]) f[i] = max(f[i], f[j] + 1)，那么所求的最长上升子序列就是 f 数组中的最大值。求完 f[i]数组后用 mx 记录 f[i]的最大值。最后输出 mx 即可。

状态转移方程为：(a[j] < a[i])f[i] = max(f[i], f[j]+1)

Code:

#include<bits/stdc++.h>
using namespace std;
const int maxn = 1010;
int w[maxn];
int f[maxn];
int main()
{
	int n;
	cin >> n;
	for (int i = 0; i < n; i++) cin >> w[i];
	int mx = 1;// 找出所计算的f[i]之中的最大值，边算边找
	for (int i = 0; i < n; i++) {
		f[i] = 1;  //f[i]默认为1，即自己
		for (int j = 0; j < i; j++) {
			if (w[j] < w[i]) {
				f[i] = max(f[i], f[j] + 1);
			}
		}
		mx = max(mx, f[i]);
	}
	cout << mx << endl;
	return 0;
}

二、 LCS(最长公共子序列)问题

问题描述

字符序列的子序列:是指从给定字符序列中随意地（不一定连续）去掉若干个字符（可能一个也不去掉）后所形成的字符序列。令给定的字符序列X=“x0，x1，…，xm-1”，序列Y=“y0，y1，…，yk-1”是X的子序列，存在X的一个严格递增下标序列<i0，i1，…，ik-1>，使得对所有的j=0，1，…，k-1，有xij=yj。例如，X=“ABCBDAB”，Y=“BCDB”是X的一个子序列。

输入两个字符串，输出它们的最长公共子串长度，若不存在公共子串则输出 No。

其中字符串长度不超过 2e5。

 问题分析详见http://t.csdn.cn/Zq5IG

 Sample Input 1

alsdfkjfjkdsal
fdjskalajfkdsla

Sample Output 1

8

Sample Input 2

ABCPDSFJGODIHJOFDIUSHGD
OSDIHGKODGHBLKSJBHKAGHI

Sample Output 2 

9

Code(Copy) :

z-k

#include <iostream>
#include <string>
#include <stack>
using namespace std;
void LCS(string s1, string s2)
{
    int m = s1.length() + 1;
    int n = s2.length() + 1;
    int** c;
    int** b;
    c = new int* [m];
    b = new int* [m];
    for (int i = 0; i < m; i++)
    {
        c[i] = new int[n];
        b[i] = new int[n];
        for (int j = 0; j < n; j++)
            b[i][j] = 0;
    }
    for (int i = 0; i < m; i++)
        c[i][0] = 0;
    for (int i = 0; i < n; i++)
        c[0][i] = 0;
    for (int i = 0; i < m - 1; i++)
    {
        for (int j = 0; j < n - 1; j++)
        {
            if (s1[i] == s2[j])
            {
                c[i + 1][j + 1] = c[i][j] + 1;
                b[i + 1][j + 1] = 1;          //1表示箭头为  左上
            }
            else if (c[i][j + 1] >= c[i + 1][j])
            {
                c[i + 1][j + 1] = c[i][j + 1];
                b[i + 1][j + 1] = 2;          //2表示箭头向  上
            }
            else
            {
                c[i + 1][j + 1] = c[i + 1][j];
                b[i + 1][j + 1] = 3;          //3表示箭头向  左
            }
        }
    }
    for (int i = 0; i < m; i++)                //输出c数组
    {
        for (int j = 0; j < n; j++)
        {
            cout << c[i][j] << ' ';
        }
        cout << endl;
    }
    stack<char> same;                   //存LCS字符
    stack<int> same1, same2;             //存LCS字符在字符串1和字符串2中对应的下标，方便显示出来
    for (int i = m - 1, j = n - 1; i >= 0 && j >= 0; )
    {
        if (b[i][j] == 1)
        {
            i--;
            j--;
            same.push(s1[i]);
            same1.push(i);
            same2.push(j);
        }
        else if (b[i][j] == 2)
            i--;
        else
            j--;
    }
    cout << s1 << endl;                     //输出字符串1
    for (int i = 0; i < m && !same1.empty(); i++)      //输出字符串1的标记
    {
        if (i == same1.top())
        {
            cout << 1;
            same1.pop();
        }
        else
            cout << ' ';
    }
    cout << endl << s2 << endl;                //输出字符串2
    for (int i = 0; i < n && !same2.empty(); i++)      //输出字符串2的标记
    {
        if (i == same2.top())
        {
            cout << 1;
            same2.pop();
        }
        else
            cout << ' ';
    }
    cout << endl << "最长公共子序列为：";
    while (!same.empty())
    {
        cout << same.top();
        same.pop();
    }
    cout << endl << "长度为：" << c[m - 1][n - 1] << endl;
    for (int i = 0; i < m; i++)
    {
        delete[] c[i];
        delete[] b[i];
    }
    delete[]c;
    delete[]b;
}
int main()
{
    string s1 = "ABCPDSFJGODIHJOFDIUSHGD";
    string s2 = "OSDIHGKODGHBLKSJBHKAGHI";
    LCS(s1, s2);
    return 0;
}