实现最长公共子序列问题

问题描述：
  最长公共子序列（LCS）是一个在一个序列集合中（通常为两个序列）用
  来查找所有序列中最长子序列的问题。这与查找最长公共子串的问题同
  的地方是：子序列不需要在原序列中占用连续的位置。而最长公共子串
 （要求连续）和最长公共子序列是不同的。

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
int longestSequence[50][50] = {0};//最长子序列表
int main(){
    puts("--------------动态规划-----最长公共子序列问题----------------\n");
     int *sequence1,*sequence2;
     int length1,length2;
     puts("输入序列1的长度:");
     scanf("%d",&length1);
     puts("输入序列1每个数字:");
     sequence1 = (int *)malloc(sizeof(int)*(length1+1));
     int i = 1;
     for(i; i <= length1;i++){
        scanf("%d",&sequence1[i]);
     }
     puts("输入序列2的长度:");
     scanf("%d",&length2);
     puts("输入序列1每个数字:");
     sequence2 = (int *)malloc(sizeof(int)*(length2+1));
     i = 1;
     for(i; i <= length2;i++){
        scanf("%d",&sequence2[i]);
     }
     i = 0;
     dynamicPrograming(sequence1,length1,sequence2,length2);
     for(i;i<=length1;i++){
        int j=0;
        for(j;j<=length2;j++){
            printf("f[%d][%d]=%d\n",i,j,longestSequence[i][j]);
        }
     }
    printf("序列1和序列2的最长公共子序列长度为:%d\n",longestSequence[length1][length2]);
    return 0;
}

void dynamicPrograming( int *sequence1,int length1,int *sequence2,int length2){
    //思想主要是两个序列的最后一个字符相等的话，那么问题规模从f(m,n)降为f(m-1,n-1);f(m,n)=f(m-1,n-1)+1，f(m,n)代表两个子序列的最长长度
    //否者不相等那么问题f(m,n) = max{f(m,n-1),f(m-1,n)} 

    //先处理好所有只有f[m][n]中有一个长度为1的序列的的问题，这时候需要比较长度为1的那个序列是否在长度>1的序列存在
    int sequence2Count = 1;
    int sequence1Count = 1;
    for(sequence2Count;sequence2Count<=length2;sequence2Count++){
        int i = 1;
        //遍历寻找是sequence1[1]否存在当前的长度sequence2中
        for(i;i<=sequence2Count;i++){
            if(sequence1[1]==sequence2[i]){
                longestSequence[1][sequence2Count] = 1;
                break;
            }
        }
    }

    for(sequence1Count;sequence1Count<=length1;sequence1Count++){
        int i = 1;
        //遍历寻找是sequence1[1]否存在当前的长度sequence2中
        for(i;i<=sequence1Count;i++){
            if(sequence2[1]==sequence1[i]){
                longestSequence[sequence1Count][1] = 1;
                break;
            }
        }
    }

    sequence1Count = 2;
    for(sequence1Count;sequence1Count <= length1;sequence1Count++){
        sequence2Count = 2;
        for(sequence2Count; sequence2Count <= length2;sequence2Count++){
            if(sequence1[sequence1Count] == sequence2[sequence2Count]){//如果当前的两个序列最后两个数字相等

                longestSequence[sequence1Count][sequence2Count] = longestSequence[sequence1Count-1][sequence2Count-1] + 1;
            } else {
                //有可能存在f[1][2]但是最后两个数字都不相等那么就会让f[1][2] = f[0][1];这时候就等于0了 这是一个bug
                int f1 = longestSequence[sequence1Count-1][sequence2Count];
                int f2 = longestSequence[sequence1Count][sequence2Count-1];
                longestSequence[sequence1Count][sequence2Count] = f1>f2?f1:f2;
            }
        }
    }
}