HAOI2010最长公共子序列

题目描述

字符序列的子序列是指从给定字符序列中随意地（不一定连续）去掉若干个字符（可能一个也不去掉）后所形成的字符序列。令给定的字符序列X=“x₀，x₁，…，x_{m-1”，序列Y=}“y₀，y₁，…，y_k-1”是X的子序列，存在X的一个严格递增下标序列<i₀，i₁，…，i_k-1>，使得对所有的j=0，1，…，k-1，有x_ij = y_j。例如，X=“ABCBDAB”，Y=“BCDB”是X的一个子序列。

对给定的两个字符序列，求出他们最长的公共子序列长度，以及最长公共子序列个数。

题目分析

第一问非常简单吧，f[i][j]表示A序列匹配到i，B序列匹配到j，那么如果A[i]=B[j],那么f[i][j]=f[i-1][j-1];如果A[i]!=B[i].f[i][j]=max(f[i-1][j],f[i][j-1]);

但是第二问就没这么容易了，如果A[i]=B[j]。

样例： ABB和BCB

首先，num[i][j]=num[i-1][j-1]是一定的吧，但是有可能A[i]和B[j-1]相等，或者A[i-1]和B[j]相等，所以在相等时要进行加的操作。

另外的情况

样例：ABD和ABC

如果f[i-1][j]=f[i][j-1]的时候，num[i][j]=num[i-1][j]+f[i][j-1]，但是样例中的情况，显然重复计算了一次不和A[i]与B[j]产生影响的“AB”，所以要减去num[i-1][j-1]。

代码

#include<iostream>
#include<algorithm>
#include<string>
#include<vector>
#include<cstdlib>
#include<cstdio>
#include<cstring>
#include<queue>
#include<climits>
using namespace std;
int f[5005],last[5005],num[5005],ll[5005];
char s1[5005],s2[5005];
int mod=100000000;
int n1,n2;
int main()
{
    int i,j,ans,bj,kl;
    scanf("%s",s1);scanf("%s",s2);
    n1=strlen(s1);n2=strlen(s2);
    n1--;n2--;
    for(i=0;i<=n2;i++)ll[i]=1;
    for(i=0;i<n1;i++){
    	for(j=0;j<n2;j++){
    		ll[0]=1;num[0]=1;
    		if(s1[i]==s2[j]){
    			f[j+1]=last[j]+1;num[j+1]=ll[j];
    			num[j]=num[j]%mod;
    			if(last[j+1]==f[j+1])num[j+1]+=ll[j+1];
    			num[j]=num[j]%mod;
    			if(f[j]==f[j+1])num[j+1]+=num[j];
    			num[j]=num[j]%mod;
    		}
    		else {
    			if(f[j]>last[j+1]){
    				f[j+1]=f[j];num[j+1]=num[j];
    				num[j]=num[j]%mod;
    			}
    			else if(f[j]<last[j+1]){
    				f[j+1]=last[j+1];num[j+1]=ll[j+1];
    				num[j]=num[j]%mod;
    			}
    			else {
    				f[j+1]=f[j];
    				num[j+1]=num[j]+ll[j+1];
    				if(last[j]==f[j+1])num[j+1]-=ll[j];
    				num[j]=num[j]%mod;
    			}
    		}
    		num[j]=num[j]%mod;
    	}
    	for(j=0;j<n2;j++){last[j+1]=f[j+1];ll[j+1]=num[j+1];}
    }
    printf("%d\n%d",f[n2],num[n2]);
    return 0;
}