最长公共子序列

算法设计与分析课程实验报告

实验名称：             班级：            学号：              姓名：             指导老师：

• 实验目的

掌握使用动态规划法解决最长公共子序列问题。

• 实验关键信息记录

#include<stdio.h>

#include <iostream>   

using namespace std;

const int M = 8;

const int N = 7;

void output(char* s, int n);

void LCSLength(int m, int n, char* x, char* y, int** c);

void LCS(int i, int j, char* x, int** c);

int main()

{

    char x[] = { ' ','B','D','C','B','D','B','C','A' };

    char y[] = { ' ','D','B','A','C','B','A','B' };

    int** c = new int* [M + 1];

    for (int i = 0; i <= M; i++)

    {

        c[i] = new int[N + 1];

    }

    cout << "序列X：" << endl;

    output(x, M);

    cout << "序列Y：" << endl;

    output(y, N);

    LCSLength(M, N, x, y, c);

    cout << "序列X、Y最长公共子序列长度为：" << c[M][N] << endl;

    cout << "序列X、Y最长公共子序列为：" << endl;

    LCS(M, N, x, c);

    cout << endl;

}

void output(char* s, int n)

{

    for (int i = 1; i <= n; i++)

    {

        cout << s[i] << " ";

    }

    cout << endl;

}

void LCSLength(int m, int n, char* x, char* y, int** c)

{

    int i, j;

    for (i = 1; i <= m; i++)

        c[i][0] = 0;

    for (i = 1; i <= n; i++)

        c[0][i] = 0;

    for (i = 1; i <= m; i++)

    {

        for (j = 1; j <= n; j++)

        {

            if (x[i] == y[j])

            {

                c[i][j] = c[i - 1][j - 1] + 1;

            }

            else if (c[i - 1][j] >= c[i][j - 1])

            {

                c[i][j] = c[i - 1][j];

            }

            else

            {

                c[i][j] = c[i][j - 1];

            } } }}

void LCS(int i, int j, char* x, int** c)

{

    if (i == 0 || j == 0)

    {

        return;

    }

    if (c[i][j] == c[i - 1][j - 1] + 1)

    {

        LCS(i - 1, j - 1, x, c);

        cout << x[i] << " ";

    }

    else if (c[i - 1][j] >= c[i][j - 1])

    {

        LCS(i - 1, j, x, c);

    }

    else

    {

        LCS(i, j - 1, x, c);

    }

}

结果：

• 实验总结