BZOJ2423 求最长公共子序列的长度和个数

题意：

给你两个字符串，求这两个字符串最长公共子序列的长度(LCS)和个数。

思路：

字符串a[max],b[max]
求长度的数组dp[2][max]
（一）求长度：
其实可以直接看白皮p56
基础的动态规划；设一个数组dp[i][j] ,表示从起始位到第一个字符串的第i位与第二个字符串的第j位的最大长度。
遍历两个字符串:
如果a[i]==b[j] 则 dp[i][j]=dp[i-1][j-1]+1;
如果a[i]!=b[j] 则 dp[i][j]=max(dp[i-1][j],dp[i][j-1]);

设求个数的数组g[2][max]
(二)求个数：
稍微绕个弯，在理解LCS的基础上稍微改一下。
（1）如果a[i]==b[j]：
此时公共子序列的长度增加1，其子序列个数至少与g[i-1][j-1]的个数相同 (只是在上次的子序列中都加了一个字母，所以个数不变)
然后去找第一个串的上一个和第二个串的上一个。
①如果dp[i][j]==dp[i][j-1] (即两者的长度相同,但内容不同或顺序不同) 则g[i][j]+=g[i][j-1];
②如果dp[i][j]==dp[i-1][j] (同上) 则g[i][j]+=g[i-1][j];
（2）如果a[i]!=b[j]：
g[i][j]=k1g[i-1][j]+k2g[i][j-1]-k3*g[i-1][j-1]
与（1）时原理相同
①如果dp[i][j]==dp[i][j-1] (即k1=1,否则为0) 则g[i][j]+=g[i][j-1];
②如果dp[i][j]==dp[i-1][j] (即k2=1,否则为0) 则g[i][j]+=g[i-1][j];
③如果dp[i][j]==dp[i-1][j-1] (即k3=1,否则为0) 则g[i][j]-=g[i-1][j-1];

③的情况是因为如果①和②同时触发则会出现加了两次g[i-1][j-1]的情况，所以要剪掉。
例：abca和abcd：
当i=4，j=4时，dp[i][j]=dp[i-1][j-1]=3，g[i][j-1]=g[i-1][j]=1，
所以表示的为同一 / 公共子序列 / 的个数，所以要减去其中一个。

滚动数组：
百度，请
因为只涉及到i与i-1和j与j-1的操作，所以可以不用将所有的i都记录下来，因为是扫描字符串，所以用1表示i，0表示i-1就行。

取模：
详见https://blog.csdn.net/qq_37025443/article/details/79339724

#include <cstdio>
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <vector>
#include <cstring>
#include <stack>
#include <cmath>
using namespace std;
const long long mod=1e8;
const long long inf=0x3f3f3f3f;
char a[5010],b[5010];
int  dp[2][5010],g[2][5010];
int main()
{
    scanf("%s%s",a+1,b+1);
    int l1=strlen(a+1)-1,l2=strlen(b+1)-1;
    for(int i=0;i<=l2;i++)
       g[0][i]=1;
    g[1][0]=1;
    for(int i=1;i<=l1;i++){
    	int nw=i&1,pr=nw^1;    //n&1,n为奇数时n&1=1；偶数时为0
    	                                    //1^1=0;  0^1 =1, 就是异或运算
    	for(int j=1;j<=l2;j++){
    		if(a[i]==b[j])
    		{
    			dp[nw][j]=dp[pr][j-1]+1;
	    		g[nw][j]=g[pr][j-1];
	    		if(dp[nw][j]==dp[nw][j-1])
	    		  g[nw][j]=(g[nw][j]+g[nw][j-1])%mod;
	    		if(dp[nw][j]==dp[pr][j])
	    		  g[nw][j]=(g[nw][j]+g[pr][j])%mod;
    		}
    		else
    		{
    			dp[nw][j]=max(dp[pr][j],dp[nw][j-1]);
    			g[nw][j]=0;
    			if(dp[nw][j]==dp[nw][j-1])
	    		  g[nw][j]=(g[nw][j]+g[nw][j-1])%mod;
	    		if(dp[nw][j]==dp[pr][j])
	    		  g[nw][j]=(g[nw][j]+g[pr][j])%mod;
	    		if(dp[nw][j]==dp[pr][j-1])
	    		  g[nw][j]=(g[nw][j]-g[pr][j-1])%mod;
    		}
    	}	
    }
    cout<<dp[l1&1][l2]%mod<<endl;
    cout<<(g[l1&1][l2]+mod)%mod<<endl;   
	return 0;
}