LCS

来源：互联网发布：两会期间网络保障方案编辑：程序博客网时间：2024/06/15 07:19
#include <stdio.h>#include <string.h>#include <stdlib.h>#define MAXN 9999int num[MAXN][MAXN];int step[MAXN][MAXN];char res[MAXN];void get_lcs(char s1[],char s2[],int i,int j){    int m = num[j][i]-1;    while(i>=0&&j>=0){        if(step[j][i]==0){            res[m--] = s1[i];            i--;            j--;        }        else if(step[j][i] == -1)            i--;        else            j--;    }}int calc_lcs(char s1[],char s2[],int i,int j){    if (i==-1 || j==-1)        return 0;    if (num[j][i]>-1)        return num[j][i];    int n1,n2,n;    if(s1[i] == s2[j]){        n = calc_lcs(s1,s2,i-1,j-1);        step[j][i] = 0;        num[j][i] = n+1;        return n+1;    }    else{        n1 = calc_lcs(s1,s2,i-1,j);        n2 = calc_lcs(s1,s2,i,j-1);        if(n1>n2){            step[j][i] = -1;            num[j][i] = n1;            return n1;        }        else{            step[j][i] = 1;            num[j][i] = n2;            return n2;        }    }}double calc(char s1[],char s2[]){    int len1,len2,i,j;    len1 = strlen(s1);    len2 = strlen(s2);    for(j=0;j<len2;j++)        for(i=0;i<len1;i++){num[j][i] = -1;            step[j][i] = 9;}    //int l = num[len2-1][len1-1];    int l = calc_lcs(s1,s2,len1-1,len2-1);    int shorter = (len1 < len2 ? len1 : len2);    double percent = (double)l/(double)shorter;    return percent;}int main(){/* *fgets(s1,MAXN,file);fgets(s2,MAXN,file);*/    /*    char s1[] = "asdf";    char s2[] = "aqweqwedqrqef";*/char s1[MAXN];char s2[MAXN];FILE *instream;instream = freopen("text","r",stdin);if (instream == NULL)fprintf(stderr,"freopen failed！！！！！！！！！");gets(s1);gets(s2);fclose(stdin);instream = NULL;    double p = calc(s1,s2);printf("%s\n",s1);printf("%s\n",s2);    printf("%f\n",p);    return 0;}
0 0