2019牛客暑期多校训练营（第五场）G：subsequence 1（递推）

【题解】

题意：输出S中所有大于T的数字子序列的个数。

思路：首先考虑到如果选取的数字子序列的长度大于m并去除前导0的干扰，那么肯定成立。所以比较复杂的是长度相等的情况，这里我用了递推。

dp[i][j] 表示当S推导到位置 i 时，长度为 j 的数字子序列的个数，我们只要在长度为1的时候特殊处理0的情况就可以排除前导0的干扰了。递推式：dp[i][j]=dp[i-1][j]+dp[i-1][j-1] (当 j >1 时)

dp1[i][j]表示当S推导到位置 i 时，长度为 j 的已经大于长度为 j 的T的数字子序列的个数，当然光是考虑这个是不够的，因为会有前面相等，后面大于，结合一下也是大于的情况，所以我们需要定义一个数组存储前面相等的情况个数。

pre[i][j]表示当S推导到位置 i 时，长度为 j 的等于长度为 j 的T的数字子序列的个数。

那么我们在什么时候增加大于的情况个数呢？

首先，当dp1[i-1][j-1]>0时，就是前面j-1个长度已经大于了，那么后面不管是什么都还是大于，所以这是每次都要加上的。如果当前位置不被选中，那么dp[i][j]=dp[i-1][j]，这个应该好理解的吧，直接顺延继承就好了。其次，当S[i]==T[j]时，更新相同的情况；当S[i]>T[j]时，前面的相同的部分就派上用场了，dp[i][j]+=pre[i-1][j-1]。

【代码】

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const ll mod=998244353;
char s[3005],t[3005];
ll dp[3005][3005];
ll dp1[3005][3005];
ll pre[3005][3005];

int main()
{
    int T; scanf("%d",&T);
    while(T--){
        int n,m; scanf("%d%d",&n,&m);
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                dp[i][j]=0,dp1[i][j]=0,pre[i][j]=0;
        scanf("%s%s",s+1,t+1);
        for(int i=1;i<=n;i++){
            for(int j=1;j<=i;j++){
                if(j==1){
                    if(s[i]=='0') dp[i][j]=dp[i-1][j];
                    else dp[i][j]=dp[i-1][j]+1;
                }
                else dp[i][j]=dp[i-1][j]+dp[i-1][j-1];
                dp[i][j]=dp[i][j]%mod;
            }
        }
        for(int i=1;i<=n;i++){
            for(int j=1;j<=m;j++){
                if(j==1){
                    pre[i][j]=pre[i-1][j],dp1[i][j]=dp1[i-1][j];
                    if(s[i]==t[j])
                        pre[i][j]+=1;
                    else if(s[i]>t[j])
                        dp1[i][j]+=1;
                }
                else{
                    pre[i][j]=pre[i-1][j],dp1[i][j]=dp1[i-1][j]+dp1[i-1][j-1];
                    if(s[i]==t[j])
                        pre[i][j]+=pre[i-1][j-1];
                    else if(s[i]>t[j])
                        dp1[i][j]+=pre[i-1][j-1];
                }
                 dp1[i][j]=dp1[i][j]%mod;
                 pre[i][j]=pre[i][j]%mod;
            }
        }
        ll sum=dp1[n][m];
        for(int i=m+1;i<=n;i++)
            sum=(sum+dp[n][i])%mod;
        printf("%lld\n",sum);
    }
}