【JZOJ5237】【GDOI2018模拟8.7】最长公共子序列

Description

Data Constraint

对于20%的数据，n,m<=10
对于40%的数据，n,m<=20
对于60%的数据，n,m<=100
对于80%的数据，n,m<=1000
对应100%的数据，n,m<=5000，保证序列只包含小写字母。

Solution

这道题和之前做的一道题好像啊，然后就往那个方向去想，搞了1个半小时，最后才发现这道题是弱化版，根本不需要之前的做法。
我们设f[i][j]表示当前A序列1~i，B序列1~j的最长不下降子序列，g[i][j]表示在f[i][j]长度下的A序列1~i，B序列1~j的好的二元组的数量。那么我们分别考虑i,j，判断若i选的f[i][j-1]==f[i][j],则g[i][j]+=g[i][j-1]，j同理，同时要判断若s[i]==s1[j]，还要判断i，j都选的情况，但这样有可能会算重，说以要判断若i，j都不选f[i-1][j-1]==f[i][j],那么g[i][j]-=g[i-1][j-1]。

Code

#include<iostream>
#include<cmath>
#include<cstring>
#include<cstdio>
#include<algorithm>
#define ll long long
using namespace std;
const ll maxn=5e3+5,mo=1e9+7;
ll f[maxn][maxn],g[maxn][maxn],a[maxn];
ll n,m,i,t,j,k,l,x,y,z,p;
char s[maxn],s1[maxn];
int main(){
    freopen("lcs.in","r",stdin);freopen("lcs.out","w",stdout);
    scanf("%s\n",s+1);n=strlen(s+1);
    scanf("%s\n",s1+1);m=strlen(s1+1);
    for (i=1;i<=m;i++)g[0][i]=1;
    for (i=1;i<=n;i++)g[i][0]=1;
    g[0][0]=1;
    for (i=1;i<=n;i++){p=0;
        for (j=1;j<=m;j++){
            f[i][j]=f[i-1][j-1]+(s[i]==s1[j]);
            f[i][j]=max(max(f[i][j],f[i][j-1]),f[i-1][j]);
            if (f[i-1][j]==f[i][j])g[i][j]=g[i-1][j];
            if (f[i-1][j-1]+1==f[i][j] && s[i]==s1[j]) g[i][j]=(g[i][j]+g[i-1][j-1])%mo;
            if (f[i][j-1]==f[i][j]) g[i][j]=(g[i][j]+g[i][j-1])%mo;
            if (f[i-1][j-1]==f[i][j]) g[i][j]=(g[i][j]-g[i-1][j-1]+mo)%mo;
        }
    }
    printf("%lld\n",f[n][m]);
    printf("%lld\n",g[n][m]);
}