最長公共子序列

代碼一是讓你輸入兩個序列,然後輸出最長公共子序列和長度。

代碼二是讓你輸入三個序列,然後輸出最長公共子序列的長度。

#include <stdio.h>
#include <string.h>
#include <stdlib.h>
int LCSLength(char* str1, char* str2, int **b)
{
    int i,j,length1,length2,len;
    length1 = strlen(str1);
    length2 = strlen(str2);
 
    //雙指針的方法申請動態二維數組
    int **c = new int*[length1+1]; //共有length1+1行
    for(i = 0; i < length1+1; i++)
        c[i] = new int[length2+1];//共有length2+1列
 
    for(i = 0; i < length1+1; i++)
        c[i][0]=0;        //第0列都初始化爲0
    for(j = 0; j < length2+1; j++)
        c[0][j]=0;        //第0行都初始化爲0
 
    for(i = 1; i < length1+1; i++)
    {
        for(j = 1; j < length2+1; j++)
        {
            if(str1[i-1]==str2[j-1])//由於c[][]的0行0列沒有使用,c[][]的第i行元素對應str1的第i-1個元素
            {
                c[i][j]=c[i-1][j-1]+1;
                b[i][j]=0;          //輸出公共子串時的搜索方向
            }
            else if(c[i-1][j]>c[i][j-1])
            {
                c[i][j]=c[i-1][j];
                b[i][j]=1;
            }
            else
            {
                c[i][j]=c[i][j-1];
                b[i][j]=-1;
            }
        }
    }
    /*
    for(i= 0; i < length1+1; i++)
    {
    for(j = 0; j < length2+1; j++)
    printf("%d ",c[i][j]);
    printf("\n");
    }
    */
    len=c[length1][length2];
    for(i = 0; i < length1+1; i++)    //釋放動態申請的二維數組
        delete[] c[i];
    delete[] c;
    return len;
}
void PrintLCS(int **b, char *str1, int i, int j)
{
    if(i==0 || j==0)
        return ;
    if(b[i][j]==0)
    {
        PrintLCS(b, str1, i-1, j-1);//從後面開始遞歸,所以要先遞歸到子串的前面,然後從前往後開始輸出子串
        printf("%c",str1[i-1]);//c[][]的第i行元素對應str1的第i-1個元素
    }
    else if(b[i][j]==1)
        PrintLCS(b, str1, i-1, j);
    else
        PrintLCS(b, str1, i, j-1);
}
 
int main(void)
{
    char str1[100],str2[100];
    int i,length1,length2,len;
    printf("請輸入第一個字符串:");
    gets(str1);
    printf("請輸入第二個字符串:");
    gets(str2);
    length1 = strlen(str1);
    length2 = strlen(str2);
    //雙指針的方法申請動態二維數組
    int **b = new int*[length1+1];
    for(i= 0; i < length1+1; i++)
        b[i] = new int[length2+1];
    len=LCSLength(str1,str2,b);
    printf("最長公共子序列的長度爲:%d\n",len);
    printf("最長公共子序列爲:");
    PrintLCS(b,str1,length1,length2);
    printf("\n");
    for(i = 0; i < length1+1; i++)//釋放動態申請的二維數組
        delete[] b[i];
    delete[] b;
    system("pause");
    return 0;
}

 

 

#include <stdio.h>
#include <string.h>
#include <stdlib.h>
int max1(int m,int n)
{
    if(m>n)
        return m;
    else
        return n;
}
int max2(int x,int y,int z,int k,int m,int n)
{
    int max=-1;
    if(x>max)
        max=x;
    if(y>max)
        max=y;
    if(z>max)
        max=z;
    if(k>max)
        max=k;
    if(m>max)
        max=m;
    if(n>max)
        max=n;
    return max;
}
int LCSLength(char* str1, char* str2, char* str3) //求得三個字符串的最大公共子序列長度並輸出公共子序列
{
    int i,j,k,length1,length2,length3,len;
    length1 = strlen(str1);
    length2 = strlen(str2);
    length3 = strlen(str3);
 
    //申請動態三維數組
    int ***c = new int**[length1+1];      //共有length1+1行
    for(i = 0; i < length1+1; i++)
    {
        c[i] = new int*[length2+1];      //共有length2+1列
        for(j = 0; j<length2+1; j++)
            c[i][j] = new int[length3+1];
    }
 
    for(i = 0; i < length1+1; i++)
    {
        for(j = 0; j < length2+1; j++)
            c[i][j][0]=0;
    }
    for(i = 0; i < length2+1; i++)
    {
        for(j = 0; j < length3+1; j++)
            c[0][i][j]=0;
    }
    for(i = 0; i < length1+1; i++)
    {
        for(j = 0; j < length3+1; j++)
            c[i][0][j]=0;
    }
 
    for(i = 1; i < length1+1; i++)
    {
        for(j = 1; j < length2+1; j++)
        {
            for(k = 1; k < length3+1; k++)
            {
                if(str1[i-1]==str2[j-1] && str2[j-1]==str3[k-1])
                    c[i][j][k]=c[i-1][j-1][k-1]+1;
                else if(str1[i-1]==str2[j-1] && str1[i-1]!=str3[k-1])
                    c[i][j][k]=max1(c[i][j][k-1],c[i-1][j-1][k]);
                else if(str1[i-1]==str3[k-1] && str1[i-1]!=str2[j-1])
                    c[i][j][k]=max1(c[i][j-1][k],c[i-1][j][k-1]);
                else if(str2[j-1]==str3[k-1] && str1[i-1]!=str2[j-1])
                    c[i][j][k]=max1(c[i-1][j][k],c[i][j-1][k-1]);
                else
                {
                    c[i][j][k]=max2(c[i-1][j][k],c[i][j-1][k],c[i][j][k-1],c[i-1][j-1][k],c[i-1][j][k-1],c[i][j-1][k-1]);
                }
            }
        }
    }
    len=c[length1][length2][length3];
    for(i = 1; i < length1+1; i++) //釋放動態申請的三維數組
    {
        for(j = 1; j < length2+1; j++)
            delete[] c[i][j];
        delete[] c[i];
    }
    delete[] c;
    return len;
}
 
int main(void)
{
    char str1[100],str2[100],str3[100];
    int len;
 
    printf("請輸入第一個字符串:");
    gets(str1);
    printf("請輸入第二個字符串:");
    gets(str2);
    printf("請輸入第三個字符串:");
    gets(str3);
    len=LCSLength(str1,str2,str3);
    printf("最長公共子序列的長度爲:%d\n",len);
    system("pause");
    return 0;
}
發佈了29 篇原創文章 · 獲贊 14 · 訪問量 4147
發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章