tyoi_c(lcs)

参考这位大佬
一、最长公共子序列

int dp[1000][1000]
int Lcs(string x,string y)
{
	for(int i=0;i<=x.length();i++)
	  for(int j=0;j<=y.length();j++)
	    if(i==0||j==0)
	      dp[i][j]=0;
	    else if(x[i-1]==y[j-1])
	      dp[i][j]=dp[i-1][j-1]+1;
	    else
	      dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
	return dp[x.length()][y.length()];
}

二、最长公共子串

int dp[1000][1000]
int Lcs(string x,string y)
{
	int ans=0;
	for(int i=0;i<=x.length();i++)
	  for(int j=0;j<=y.length();j++)
	    if(i==0||j==0||x[i-1]!=y[j-1])
	      dp[i][j]=0;
	    else{
	       dp[i][j]=dp[i-1][j-1]+1;
	       ans=max(ans,dp[i][j]);
		}
	return ans;
}

回到题目

代码：

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
using namespace std;
const int mod=1e9+7;
char a[1005],b[1005];
int dp[1005][1005],f[1005][1005];
int pre[1005][27],pos[27];
int main()
{
	scanf("%s%s",a+1,b+1);
	int n=strlen(a+1),m=strlen(b+1);
	for(int i=1;i<=n;i++)
	for(int j=1;j<=m;j++)
	{
		dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
		if(a[i]==b[j]) dp[i][j]=max(dp[i][j],dp[i-1][j-1]+1);
	}
	memset(pos,-1,sizeof(pos));
	for(int i=1;i<=m;i++)
	{
		pos[b[i]-'a']=i;
		for(int j=0;j<26;j++) pre[i][j]=pos[j];
	}
	for(int i=0;i<=n;i++) f[i][0]=1;
	for(int i=0;i<=m;i++) f[0][i]=1;
	for(int i=1;i<=n;i++)
	for(int j=1;j<=m;j++)
	{
		f[i][j]=0;
		if(dp[i][j]==dp[i-1][j]) 
		{
			f[i][j]+=f[i-1][j];
			if(f[i][j]>=mod) f[i][j]-=mod;
		}
		int p=pre[j][a[i]-'a'];
		if(p!=-1 && dp[i-1][p-1]+1==dp[i][j])
		{
			f[i][j]+=f[i-1][p-1];
			if(f[i][j]>=mod) f[i][j]-=mod;
		}
	}
	printf("%d\n",f[n][m]);
}