题意:给出两个数字串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;
}
/*
*/