用dp[i][j]表示,考虑前i个字符,且t[i]!=s[i],且对于所有的k>i,都有s[k]=t[k],在这些所有的串中,满足恰有j个子串的字典序比原串大的串的个数。
那么考虑第i个串怎么填,只有两种填法,要么比s[i]大,要么比s[i]小,考虑第一种填法,那么枚举在i之前一段长为k的串,即t[i-k,i-1],满足t[i-k,i-1]==s[i-k,i-1],那么在这k个中任选一个,再在i后面n-i中任选一个,由这两个端点组成的串字典序一定大于原字串。注意可以选空,就有(k+1)*(n-i+1)种选法,又因为这个数不等大于dp[i][j]的j,所有可以用j控制对k的枚举。
填比s[i]小的,对答案没贡献,就是求出dp[1..i-1][j]的和。
#include <cstdio>
#include <cstring>
#include <algorithm>
#define maxn 2009
#define MOD 1000000007
using namespace std;
__int64 dp[maxn][maxn],sum[maxn];
char s[maxn];
int n,m;
int main()
{
scanf("%d%d",&n,&m);
scanf("%s",s+1);
dp[0][0]=1;sum[0]++;
for(int i=1;i<=n;i++)
{
for(int j=0;j<=m;j++)
{
dp[i][j]=sum[j]*(s[i]-'a')%MOD;
for(int k=0;k<i&&(k+1)*(n-i+1)<=j;k++)
{
dp[i][j]=(dp[i][j]+dp[i-k-1][j-(k+1)*(n-i+1)]*('z'-s[i]))%MOD;
}
sum[j]+=dp[i][j];
sum[j]%=MOD;
}
}
__int64 ans=0;
for(int i=0;i<=n;i++)
{
ans+=dp[i][m];
ans%=MOD;
}
printf("%I64d\n",ans);
//system("pause");
return 0;
}