动态规划——序列匹配

动态规划——序列匹配

问题：

设S1 [1,m]和S2 [1,n]为两个字符序列，对这两个字符序列进行顺序扫描比较，比较时可以在某个序列的两个字符之间插入空格，以使他们局部相同区域对应。在比较中，如果两个对应字符相等，赋权值2，如果对应字符不等，赋权值-2，如果一个是字符一个是空格，赋权值-1。

问题：寻找权值最大的匹配度（匹配权值和）与对准方式。

分析：

c[i][j]代表两序列从1分别到i,j的最大匹配值。对于每个序列的末端而言，有三种情况：1.字符对字符，二者相等，则两序列末端均-1；2.字符对字符，二者不相等，两序列末端同样-1；3.字符对空格，字符序列末端-1

因此递推公式如下：

c[i][j] = max{c±2,c[i-1][j]-1,c[i][j-1]-1}

初始化第0行第0列为0

代码：

#include <bits/stdc++.h>

using namespace std;

#define MAXN 105

// 寻找权值最大的匹配度
// 输入：两待匹配序列s1，s2，长度分别为m，n，记录最大匹配度矩阵c，过程记录矩阵a
void match(char s1[],char s2[],int m,int n,int c[][MAXN],int a[][MAXN]){
    // 初始化
    for(int i = 0;i <= m;i++)
        c[i][0] = 0;
    for(int i = 0;i <= n;i++)
        c[0][i] = 0;
    
    // 计算
    for(int i = 1;i <= m;i++){
        for(int j = 1;j <= n;j++){
            int equals = 0;
            if(s1[i] == s2[j]){
                equals = 2;
            }else{
                equals = -2;
            }
            int tmp1 = c[i-1][j-1] + equals;
            int tmp2 = c[i-1][j] - 1;
            int tmp3 = c[i][j-1] - 1;
            if(tmp1 >= tmp2 && tmp1 >= tmp3){
                c[i][j] = tmp1;
                a[i][j] = 0;
            }else if(tmp2 >= tmp1 && tmp2 >= tmp3){
                c[i][j] = tmp2;
                a[i][j] = 1;
            }else if(tmp3 >= tmp1 && tmp3 >= tmp2){
                c[i][j] = tmp3;
                a[i][j] = 2;
            }
        }
    }        
}

// 根据结果记录矩阵a构造最优解
void traceback(int a[][MAXN],int i,int j,char s1[],char s2[]){
    if(i==0){
        for(int k = 1;k <= j;k++){
            cout<<"s1: "<<"-"<<" s2: "<<s2[k]<<" ";
        }
        return ;
    }
    
    if(j==0){
        for(int k = 1;k <= i;k++){
            cout<<"s1: "<<s1[k]<<" s2: "<<"-"<<" ";
        }
        return ;
    }
    
    if(a[i][j]==0){
        traceback(a, i-1, j-1, s1, s2);
        cout<<"s1: "<<s1[i]<<" s2: "<<s2[j]<<" ";
    }else if(a[i][j]==1){
        traceback(a, i-1, j, s1, s2);
        cout<<"s1: "<<s1[i]<<" s2: "<<"-"<<" ";
    }else if(a[i][j]==2){
        traceback(a, i, j-1, s1, s2);
        cout<<"s1: "<<"-"<<" s2: "<<s2[j]<<" ";
    }
}

int main(){
    int m = 8,n = 8;
    char s1[9] = {'*','a','x','a','b','c','d','e','s'} ,s2[9] = {'*','a','x','b','a','c','f','e','s'};
    int c[m+1][MAXN];
    memset(c, 0, sizeof(c));
    int a[m+1][MAXN];
    memset(a, 0, sizeof(a));
    
    match(s1, s2, m, n, c, a);
    
    cout<<c[m][n]<<endl;
    
    traceback(a, m, n, s1, s2);
}

复杂度分析：

时间复杂度：O(mn)

空间复杂度：O(mn)