动态规划——序列匹配
问题:
设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)