Cod sursa(job #427588)

Utilizator hendrikHendrik Lai hendrik Data 28 martie 2010 01:44:45
Problema Cel mai lung subsir comun Scor 80
Compilator cpp Status done
Runda Arhiva educationala Marime 1.47 kb
#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
using namespace std;

void open(){
    freopen("cmlsc.in","r",stdin);
    freopen("cmlsc.out","w",stdout);
}

short n,m,a[1030],b[1030],maxi,idxx,idxy;
short dp[1030][1030],par[1030][1030][2];

void find_path(short x,short y,short c){
    if (c==0) return;
    find_path(par[x][y][0],par[x][y][1],dp[par[x][y][0]][par[x][y][1]]);
    if (c>dp[par[x][y][0]][par[x][y][1]]){
        if (c>1) printf(" ");
        printf("%d",a[x-1]);
    }
}

int main(){
    open();
    scanf("%d%d",&n,&m);
    for (short i=0;i<n;i++) scanf("%d",&a[i]);
    for (short i=0;i<m;i++) scanf("%d",&b[i]);
    for (short i=1;i<=n;i++){
        for (short j=1;j<=m;j++){
            if (a[i-1]==b[j-1]){
                dp[i][j]=dp[i-1][j-1]+1;
                par[i][j][0]=i-1;
                par[i][j][1]=j-1;
            }
            else {
                if (dp[i][j-1]>dp[i-1][j]){
                    dp[i][j]=dp[i][j-1];
                    par[i][j][0]=i;
                    par[i][j][1]=j-1;
                }
                else {
                    dp[i][j]=dp[i-1][j];
                    par[i][j][0]=i-1;
                    par[i][j][1]=j;
                }
            }
            if (maxi<dp[i][j]){
                maxi=dp[i][j];idxx=i;idxy=j;
            }
        }
    }
    printf("%d\n",maxi);
    find_path(idxx,idxy,maxi);
    printf("\n");
    return 0;
}