Cod sursa(job #413151)

Utilizator johnbBaranga Ionut johnb Data 7 martie 2010 20:31:03
Problema Cel mai lung subsir comun Scor 80
Compilator c Status done
Runda Arhiva educationala Marime 2.56 kb
#include <stdio.h>
#include <stdlib.h>

#define INT_S  sizeof(int)
#define CHAR_S sizeof(char)
#define nill -1

int max(int a, int b) { return a > b ? a : b; }


/*
 * determina lungimea maxima a subsecventei / subsecventelor cumune
 */
int lcs_len(int* x, int* y, int i, int j, int **c) {
    // caz banal - unul dintre siruri este vid => lungimea este 0
    if (i == nill || j == nill) {
          return 0;
    }
    
    // daca s-a mai calculat, se intoarce rezultatul calculat, stocat
    // int tabela reprezentata prin matricea c
    if (c[i][j] == nill) {   
        if(x[i] == y[j])  
                  c[i][j] = lcs_len(x, y, i - 1, j - 1, c) + 1;
        else 
                  c[i][j] = max(lcs_len(x, y, i, j - 1, c), lcs_len(x, y, i - 1, j, c));

    }
    return c[i][j];
}

/*
 * determina o subsecventa de lungime maxima, avand tabela cu lungimile 
 * tuturor subsecventelor comune de lungime maxima ale tuturor prefixelor
 * sirurilor x si y
 */
void lcs(int* x, int* y, int i, int j, int**c, int* r, int p, int len) {
     if (c[i][j] == 0 || i == nill || j == nill)   
        return;
     
     if (x[i] == y[j]) {
              r[len - p++ - 1] = x[i]; // sau y[j]
              i--;
              j--;
     }
     else {
          if (i)
             if (c[i-1][j] > c[i][j - 1])
                i--;
          else
              if (j)
                 j--;
              else
                  i--;
     }
     lcs(x, y, i, j, c, r, p, len);
}
     
     
               
    

int main() {
    freopen("cmlsc.in", "r", stdin);
    freopen("cmlsc.out", "w", stdout);
    
    // intrari :
    int
        *x, // primul sir
        *y, // al 2 - lea sir
         m, // lungimea primului sir
         n, // lungimea celui de-al doilea sir
        
    // iesiri
        *r, // subsecventa comunade lungime maxima 
        
    // auxiliare
         i, j,
      ** c;
          
    scanf("%d %d\n", &m, &n);
    x = (int*)malloc(m * INT_S);
    y = (int*)malloc(n * INT_S);    
    for (i = 0; i < m; i++)
        scanf("%d ", x + i);
    for (i = 0; i < n; i++)
        scanf("%d ", y + i);
    c = (int**)malloc(m * sizeof(int*));
    for (i = 0; i < m; i++) {
        c[i] = (int*)malloc(n * INT_S);
        for (j = 0; j < n; j++)
            c[i][j] = nill;
    }
    int len = lcs_len(x, y, m - 1, n - 1, c);
    printf("%d\n", len);
    r = (int*)malloc(len * INT_S);
    lcs(x, y, m - 1, n - 1, c ,r , 0, len);
    for (i = 0; i < len; i++)
        printf("%d ", r[i]);
    return 0;
}