程序设计第十周 B - LIS & LCS

题目描述

东东有两个序列A和B。

他想要知道序列A的LIS和序列AB的LCS的长度。

注意，LIS为严格递增的，即a1<a2<…<ak(ai<=1,000,000,000)。

Input
第一行两个数n，m（1<=n<=5,000,1<=m<=5,000）
第二行n个数，表示序列A
第三行m个数，表示序列B

Output
输出一行数据ans1和ans2，分别代表序列A的LIS和序列AB的LCS的长度

Simple Input

5 5
1 3 2 5 4
2 4 3 1 5

Simple Output

3 2

想法

求一下a的最小上升子序列的动态规划方法， 加状态f[i] 表示以 Ai 为结尾的最长上升序列的方程。
⚫ 初始化：f[1] = 1
⚫ 转移过程 :f[i]=max{f[j] | j < i &&a[j]<a[i] }+1;
⚫ 输出答案：max{f[i], i=1…n}
求公共子序列时，定义状态g[i],[j] 为 A1, A2, …, Ai 和 B1, B2, …, Bj 的 LCS 长度
⚫ 初始化：初始 g[1][0] = g[0][1] = g[0][0] = 0
⚫ 转移方程：当 Ai == Bj 时，g[i][j] = g[i-1][j-1] + 1
⚫ 否则 g[i][j] = max(g[i-1][j], g[i][j-1])

代码

#include <iostream>
#include<algorithm>
using namespace std;
int n, m, mma = 0;
int a[5005], b[5005],f[5005],g[5005][5005]={0};

void fun(int x)
{
    int t = 0;
    for (int i = 1; i < x;i++)
    {
        if(a[i]<a[x]&&f[i]>t)
            t = f[i];
    }
    f[x] =t + 1;
}
void fun2(int i,int j)
{
    if(a[i]==b[j])
        g[i][j] = g[i - 1][j - 1] + 1;
    else                  
        g[i][j] = max(g[i - 1][j], g[i][j - 1]);
}

int main()
{
        cin >> n >> m;
    for (int i = 1; i <=n;i++)
        cin >> a[i];
    for (int i = 1; i <= m; i++)
        cin >> b[i];
    f[1] = 1;
    for (int i = 2; i <= n;i++)
    {
        fun(i);
    } 
    for (int j = 1; j <= m;j++)
        for (int i = 1; i <= n; i++)
        {
            fun2(i,j);
        }
    sort(f+1, f + n+1);
    cout << f[n] << " " << g[n][m] << endl; 
    return 0;
}