2019牛客暑期多校训练营（第五场）G subsequence 1（DP+组合数学）

链接：https://ac.nowcoder.com/acm/contest/885/G

题意：T组样例，每组样例第一行给出n、m，表示字符串s的长度和字符串t的长度。接下来两行分别给出s和t。s和t都由阿拉伯数字组成。保证n>=m,并且他们的第一个字符都不是0。问s中有多少子序列大于t（当然不能有前导零）。

思路：显然，如果子序列第一位不为0，并且长度大于t，那么子序列组成的数肯定大于t，这一部分的方案可以用组合数算。那么现在就要考虑，子序列的长度等于t的长度的方案数。

状态：dp[i][j]表示从s的前i位中，选出j位数并且与t的前j位相同的子序列数目。

转移方程：

(1)s[i]==t[j],dp[i][j]=dp[i-1][j]+dp[i-1][j-1]。

(2)否则，dp[i][j]=dp[i-1][j]。

当s[i]>t[j]时，在前i-1位中选出j-1位，后面n-i位中选出m-j位即可。即贡献为。

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N = 3e3+10;
const ll mod = 998244353;
ll dp[N][N],ans,C[N][N];
char s[N],t[N];
int n,m; 
//dp[i][j]:从s串的前i个为数里选出j位数(长度为j的子序列)
//与t的前j位数相等的子序列数 
int main(void)
{
	int T;
	scanf("%d",&T);
	C[0][0]=1;
	for(int i=1;i<=N-10;i++)
	{
		C[i][0]=1;
		for(int j=1;j<=i;j++)
			C[i][j]=(C[i-1][j]+C[i-1][j-1])%mod;
	} 
	while(T--)
	{
		scanf("%d%d%s%s",&n,&m,s+1,t+1);
		ans=0;
		
		//选出子序列长度等于t的长度，要保证前面某几位相同
		//中间至少一位大于t 
		
		//初始化，选出0位数的子序列数肯定为1 
		//即一个也不选 
		for(int i=0;i<=n;i++)	
			dp[i][0]=1;
	
		for(int i=1;i<=n;i++)
			for(int j=1;j<=min(i,m);j++)
			{
				//首先不管s[i]与t[j]有怎样的大小关系，总会有dp[i-1][j]个子序列满足条件 
				//也就是不选s串中的第i位数		 
				dp[i][j]=dp[i-1][j];
				//只有s[i]==t[j]时，才能选s[i]，那要保证从s串的前i-1为中
				//选出来的j-1位数要相等 
				if(s[i]==t[j]) dp[i][j]=(dp[i][j]+dp[i-1][j-1])%mod;	
				//既然从s串中选出的子序列组成的数第j位已经大于t的第j位了，
				//那么后面无论怎么选都大于t表示的数 
				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;
}