codeforces 10 D. LCIS LCIS O(n^2)算法

题目链接

给出两个序列， 求出他们的最长公共上升子序列。

两层循环， 内层循环j， 外层i。 如果a[i] == b[j]， 那么dp[j] = max(dp[j], dp[best]+1)， best是一个指针， 指向小于j的元素中dp值最大并且b[best]的值小于a[i]的元素。

如果a[i]>b[j]， 那么就看dp[j]是否大于dp[best], 如果大于就更新best。

具体看代码。

#include <iostream>
#include <vector>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <map>
#include <set>
#include <string>
#include <queue>
#include <stack>
#include <bitset>
using namespace std;
#define pb(x) push_back(x)
#define ll long long
#define mk(x, y) make_pair(x, y)
#define lson l, m, rt<<1
#define mem(a) memset(a, 0, sizeof(a))
#define rson m+1, r, rt<<1|1
#define mem1(a) memset(a, -1, sizeof(a))
#define mem2(a) memset(a, 0x3f, sizeof(a))
#define rep(i, n, a) for(int i = a; i<n; i++)
#define fi first
#define se second
typedef pair<int, int> pll;
const double PI = acos(-1.0);
const double eps = 1e-8;
const int mod = 1e9+7;
const int inf = 1061109567;
const int dir[][2] = { {-1, 0}, {1, 0}, {0, -1}, {0, 1} };
int dp[505], a[505], b[505], pre[505], ans[505];
int main()
{
    int n, m;
    cin.tie(0);
    ios::sync_with_stdio(0);
    cout.tie(0);
    cin>>n;
    for(int i = 1; i<=n; i++)
        cin>>a[i];
    cin>>m;
    for(int i = 1; i<=m; i++)
        cin>>b[i];
    for(int i = 1; i<=n; i++) {
        int best = 0;
        for(int j = 1; j<=m; j++) {
            if(a[i] == b[j] && dp[j]<dp[best]+1) {
                dp[j] = dp[best]+1;
                pre[j] = best;
            } else if(a[i]>b[j]&&dp[j]>dp[best]) {
                best = j;
            }
        }
    }
    int maxx = 0, cnt = 0;
    for(int i = 1; i<=m; i++)
        if(dp[i]>dp[maxx])
            maxx = i;
    cout<<dp[maxx]<<endl;
    for(int i = maxx; i!=0; i = pre[i])
        ans[cnt++] = b[i];
    for(int i = cnt-1; i>=0; i--)
        cout<<ans[i]<<" ";
    return 0;
}