Cod sursa(job #897409)

Utilizator harababurelPuscas Sergiu harababurel Data 27 februarie 2013 20:24:35
Problema Cel mai lung subsir comun Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.12 kb
#include <iostream>
#include <vector>
#include <algorithm>
#include <fstream>
#define nmax 1025
using namespace std;

int a[nmax], b[nmax], dp[nmax][nmax];
int pre[nmax];
vector <int> sol;

int cmlsc(int i, int j) {
	if(i==0 || j==0) return 0;
	
	if(a[i] == b[j]) {
		pre[i] = i-1;
		return cmlsc(i-1, j-1) + 1;
	}

	return max(cmlsc(i,j-1), cmlsc(i-1,j));
}

int backtrace(int i, int j) {
	if(i==0 || j==0) return 0;
	
	if(a[i] == b[j]) {
		sol.push_back(a[i]);
		backtrace(i-1, j-1);
		return 0;
	}
	
	if(dp[i-1][j] > dp[i][j-1]) backtrace(i-1, j);
	else backtrace(i, j-1);
	
	return 0;
}


int main() {
	ifstream f("cmlsc.in");
	ofstream g("cmlsc.out");
	
	int n, m, i, j;
	
	f>>n>>m;
	for(i=1; i<=n; i++) f>>a[i];
	for(j=1; j<=m; j++) f>>b[j];
	
	for(i=0; i<=n; i++) dp[i][0] = 0;
	for(j=0; j<=m; j++) dp[0][j] = 0;

	for(i=1; i<=n; i++) {
		pre[i] = 0;
		for(j=1; j<=m; j++) {
			if(a[i] == b[j]) dp[i][j] = dp[i-1][j-1] + 1, pre[i] = true;
			else dp[i][j] = max(dp[i-1][j], dp[i][j-1]);
		}
	}
	
	g<<dp[n][m]<<"\n";
	
	backtrace(n, m);
	
	reverse(sol.begin(), sol.end());
	for(i=0; i<sol.size(); i++) g<<sol[i]<<" "; g<<"\n";
	
	return 0;
}