100道动态规划——9 UVA 10723 Cyborg Genes LCS的应用以及有趣的题干

好吧，我觉得这个题干挺有意思的，说我们半机械人将要毁灭你们啊，什么什么的。。没见过这么长的题干。。

回到正式上面来，题目要求的就是两个串的最短公共母串的长度以及公共母串的种类数。

我觉得首先马上就可以想到LCS吧，然后发现最短公共母串的长度就是str1.length+str2.length-lcs.length，然后就是种类数的问题了。

在这里我觉得就可以发现我对于LCS的认识还仅仅停留在会写而已，没有认真去理解里面的内涵，就像是单纯的把代码记录下来了而已，在这里求种类数其实和最开始求LCS的状态转移方程差不多。

假若str1[i]==str2[j]的话，相当于一次性放了2个，种类数依然是不变的

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

假若lcs_length[i-1][j]!=lcs_length[i][j-1]，这种情况对应着需要在较长的lcs后面再附加一个，因此种类数依然是不变的

variety[i][j]=max(variety[i-1][j],variety[i][j-1]);

否则就是lcs_length[i-1][j]==lcs_length[i][j-1]的情况了，此时就str1用了i个，str2用了j个的状态可以由variety[i-1][j]再放一个str1[i]，或者是variety[i][j-1]再放一个str2[j]得来，因此这个方程就是：  variety[i][j]=variety[i-1][j]+variety[i][j-1];

最后输出variety[len1][len2]即可，果然最基本的还是需要好好巩固一下。。。

#include <iostream>
#include <algorithm>
#include <cstring>

using namespace std;

char str1[35],str2[35];
unsigned dp[35][35],variety[35][35],t,len1,len2;

int main(){
	for(int i=0;i<=30;++i)
		variety[0][i]=1;
	for(int i=0;i<=30;++i)
		variety[i][0]=1;

	ios_base::sync_with_stdio(false);
	(cin>>t).get();
	for(int ca=1;ca<=t;++ca){
		cin.getline(str1+1,33);
		cin.getline(str2+1,33);

		len1=strlen(str1+1),len2=strlen(str2+1);
		if(!len1){
			cout<<"Case #"<<ca<<": "<<len2<<' '<<1<<endl;
			continue;
		}
		else if(!len2){
			cout<<"Case #"<<ca<<": "<<len1<<' '<<1<<endl;
			continue;
		}


		for(int i=1;str1[i];++i)
		for(int j=1;str2[j];++j){
			if(str1[i]==str2[j])
				variety[i][j]=variety[i-1][j-1],dp[i][j]=dp[i-1][j-1]+1;
			else if(dp[i-1][j]==dp[i][j-1])
				variety[i][j]=variety[i-1][j]+variety[i][j-1],dp[i][j]=dp[i-1][j];
			else if(dp[i-1][j]>=dp[i][j-1])
				variety[i][j]=variety[i-1][j],dp[i][j]=dp[i-1][j];
			else
				variety[i][j]=variety[i][j-1],dp[i][j]=dp[i][j-1];
		}


		cout<<"Case #"<<ca<<": "<<len1+len2-dp[len1][len2]<<' '<<variety[len1][len2]<<endl;
		for(int i=1;i<=len1;++i)
		for(int j=1;j<=len2;++j)
			dp[i][j]=variety[i][j]=0;
	}
}