实验三 动态规划算法

一、实验目的与要求

1、熟悉最长公共子序列问题的算法；

2、初步掌握动态规划算法；

3、能对设计的算法进行复杂度分析。

二、实验题目

最长公共子序列问题。给定2个序列X={x1,x2,…,xm}和Y={y1,y2,…,yn}，用动态规划算法找出X和Y的最长公共子序列，用C/C++实现该算法。

三、实验步骤

1、算法思想描述

由算法LCSLength计算得到的数组b可用于快速构造序列X和Y。首先从b[m][n]开始，依次在数组b中搜索。当在b[i][j]=1，表示Xi和Yj的最长公共子序列是由Xi-1和Yj-1的最长公共子序列在尾部加上xi所得到的子序列。当b[i][j]=2时，表示Xi和Yj的最长公共子序列与Xi-1和Yj的最长公共子序列相同。当b[i][j]=3时，表示Xi和Yj的最长公共子序列相同。

1. 实现的程序代码

#include <stdio.h>

#include <stdlib.h>

#include <string.h>

void Defind_Array(int **c,int **b,int len1,int len2);

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

void printfLCS( int i, int j,char *str, int **b);

int main()

{

    int Xm,Yn;

    printf("请分别输入X、Y两个序列的长度m和n：");

    scanf("%d%d",&Xm,&Yn);

    char x[Xm],y[Yn];

    getchar();

    printf("输入:------------------------\n    序列X的值：");

    gets(x);

    printf("    序列Y的值：");

    gets(y);

    int len1 = strlen(x),len2 = strlen(y);

    if(Xm<len1 || Yn<len2)//超出长度约束时提醒并退出;

    {

        printf("你输入的字符长度超过了约束范围,程序退出！");

        return 0;

    }//申请二维数组

    int **c = (int **)malloc(sizeof(int*)*(len1 + 1));

    int **b = (int **)malloc(sizeof(int*)*(len1 + 1));

    Defind_Array(c,b,len1,len2);

    int sum = LCSLength(strlen(x),strlen(y),x, y, c, b);//计算LCS的长度

    printf("输出:------------------------");

    printf("\n    其子序列长度为: %d\n    最长公共子序列Z为:",sum);

    printfLCS(len1,len2,x,b);

    int i;//动态内存释放

    for ( i = 0; i <= len1; i++)

    {

        free(c[i]);

        free(b[i]);

    }

    free(c);

    free(b);

    return 0;

}

void Defind_Array(int **c,int **b,int len1,int len2)

{

    int i,j;

    for( i = 0; i<= len1; i++ )  

    {

        c[i] = (int *)malloc(sizeof(int)*(len2 + 1));

        b[i] = (int *)malloc(sizeof(int)*(len2 + 1));

    }

    for ( i = 0; i<= len1; i++)

    {

        for( j = 0; j <= len2; j++)

        {

            c[i][j] = 0;

            b[i][j] = 0;

        }

    }

}

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

{

    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-1]==y[j-1])

            {

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

                b[i][j]=1;

            }

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

            {

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

                b[i][j]=2;

            }

            else

           {

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

                 b[i][j]=3;

           }

           }

            return c[m][n];

}

void printfLCS( int i, int j,char *str, int **b)

{

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

     if( b[i][j] == 1)

{

     printfLCS( i - 1, j - 1,str, b);

     printf("%c ", str[i-1]);

}

    else if ( b[i][j] == 2 )

    printfLCS(i - 1, j,str, b);

    else

    printfLCS(i , j - 1,str, b);

}

1. 运行结果图
2.