牛客多校第五场G

subsequence 1

只要处理长度等于t的.

转移方程没想出来QAQ

$dp(i,j,0)$代表到$s[i]$为止有多少个前缀序列与$t[0\cdots j]$相同

所以有$dp(i,j,0)=dp(i-1,j,0)+s[i]==t[j]?dp(i-1,j-1,0):0;$

$dp(i,j,1)$代表到$s[i]$为止有多少个子序列大于$t[0\cdots j]$

故有$dp(i,j,1)=dp(i-1,j,1)+(s[i]>t[j])*dp(i,j-1,0)+(j>0)*dp(i-1,j-1,0)$

#include<bits/stdc++.h>
using namespace std;
char s[],t[];
typedef long long ll;
ll dp[][][];
ll C[][];
ll mod=;
void init()
{
    for(int i=;i<;i++){
        C[i][i]=;
        C[i][]=;
    }
    for(int i=;i<;i++){
        for(int j=i+;j<;j++){
            C[j][i]=(C[j-][i-]+C[j-][i])%mod;
        }
    }
}
int main()
{
    int T;
    int n,m;
    ll ans;
    init();
    scanf("%d",&T);
    while(T--){
        ans=;
        scanf("%d%d",&m,&n);
        scanf("%s",s);
        scanf("%s",t);
       for(int i=;i<=max(m,n);i++)
        dp[i][][]=dp[i][][]=;
       // for(int i=0;i<=n;i++)dp[0][i]=1;
        for(int i=;i<m;i++){
            for(int j=;j<n;j++){

                dp[i+][j+][]=((dp[i][j][])*(s[i]==t[j])+dp[i][j+][])%mod;
                dp[i+][j+][]=(dp[i][j+][]+(s[i]>t[j])*dp[i][j][]+(j>)*dp[i][j][])%mod;
             //   cout<<dp[i+1][j+1][1]<<' ';
            }
           // cout<<'\n';
           // cout<<i<<" "<<ans<<'\n';
        }
        ans=dp[m][n][];
        //cout<<ans<<'\n';
        for(int k=n+;k<=m;k++)
        for(int i=;i+k-<=m;i++){
            if(s[i-]!='')
            ans=(ans+C[m-i][k-])%mod;
        }
        cout<<ans<<'\n';

    }
}