G—subsequence 1 2019牛客第五场 （组合数+dp）

题目链接：

题意：给出两个数字串s，t，找出s的子序列，满足大于t的方案数？

题解：

参考博客：

组合数学+dp。

对于长度大于t的s的子序列，都满足。

现在考虑长度等于t的s子序列。

dp[i][j]表示s的前i为，t的前j位满足的条件数。

转移方程见标程。

#include<bits/stdc++.h>
using namespace std;
 
typedef long long LL;
const int N=3010;
const LL mod=998244353;
 
LL C[N][N],dp[N][N];
 
void init()
{
    C[1][1]=C[1][0]=C[0][0]=1LL;
 
    for(int i=2;i<=3000;i++){
        C[i][0]=C[i][i]=1;
        for(int j=1;j<i;j++)
            C[i][j]=(C[i-1][j-1]+C[i-1][j])%mod;
    }
}
 
char s[N],t[N];
int main()
{
    init();
 
    int ncase;
    scanf("%d",&ncase);
 
    while(ncase--)
    {
        int n,m;
        scanf("%d%d",&n,&m);
        scanf("%s%s",s+1,t+1);
        for(int i=1;i<=n;i++)
            for(int j=1;j<=m;j++)
            dp[i][j]=0;
 
        for(int i=0;i<=n;i++)
            dp[i][0]=1;
 
        LL ans=0;
 
        for(int i=1;i<=n;i++)
        for(int j=1;j<=min(m,i);j++){
            dp[i][j]=dp[i-1][j];
            ///相同
            if(s[i]==t[j]) dp[i][j]=(dp[i][j]+dp[i-1][j-1])%mod;
            ///此时后面的n-i位可以随便弄都满足条件
            if(s[i]>t[j]) ans=(ans+dp[i-1][j-1]*C[n-i][m-j]%mod)%mod;
 
        }
 
        ///计算长度大于t的方案数
        for(int i=1;i<=n;i++){
            if(s[i]=='0') continue;
            for(int j=m;j<=n-i;j++)
                ans=(ans+C[n-i][j])%mod;
        }
 
        printf("%lld\n",ans);
 
    }
 
 
    return 0;
}
/*
*/