Cod sursa(job #396845)

Utilizator andrei.sfrentSfrent Andrei andrei.sfrent Data 15 februarie 2010 22:44:58
Problema Cel mai lung subsir comun Scor 100
Compilator c Status done
Runda Arhiva educationala Marime 1.47 kb
#include <stdio.h>

#define L1 1024
#define L2 1024
#define L 1024

#define max(a, b) (a > b ? a : b)

int sir1[L1 + 1], sir2[L2 + 1];
int l1, l2; /* lungimile sirurilor, nu imi plac literele m si n in cazul asta*/
int sc[L + 1];
int m[L1 + 1][L2 + 1]; //iniatializata deja cu 0

int lungimeCMLSC;

void citesteSiruri()
{
	FILE* fi = fopen("cmlsc.in", "r");
	fscanf(fi, "%d%d", &l1, &l2);
	int i;
	for(i = 1; i <= l1; ++i) fscanf(fi, "%d", &sir1[i]);
	for(i = 1; i <= l2; ++i) fscanf(fi, "%d", &sir2[i]);
	fclose(fi);
}

void determinaMatriceAsociata()
{
	int i, j;
	for(i = 1; i <= l1; ++i)
	{
		for(j = 1; j <= l2; ++j)
		{	
			if(sir1[i] == sir2[j]) m[i][j] = m[i - 1][j - 1] + 1;
			else m[i][j] = max(m[i - 1][j], m[i][j - 1]);
		}
	}
}

void determinaLungime()
{
	lungimeCMLSC = m[l1][l2];
}

void scrieRezultat()
{
	FILE* fo = fopen("cmlsc.out", "w");
	int i;
	fprintf(fo, "%d\n", lungimeCMLSC);
	for(i = 1; i <= lungimeCMLSC; ++i) fprintf(fo, "%d ", sc[i]);
	fprintf(fo, "\n");
	fclose(fo);
}

void determinaSubsirComun()
{
	int i = l1, j = l2; //pozitii in sir1 si sir2
	int indexSubsir = lungimeCMLSC; //prima pozitie necompletata din subsir
	
	while(i && j) /* refac drumul in matricea asociata */
	{
		if(sir1[i] == sir2[j])
		{
			sc[indexSubsir--] = sir1[i];
			i--;
			j--;
		}
		else
		{
			if(m[i][j - 1] > m[i - 1][j]) j--;
			else i--;
		}
	}
}

int main()
{
	citesteSiruri();
	determinaMatriceAsociata();
	determinaLungime();
	determinaSubsirComun();
	scrieRezultat();
	return 0;
}