http://codeforces.com/problemset/problem/360/C
令串下标从1到n,分析 t[i] 与s[i] 的三种情况:大于,等于,小于。仅当t[i]>s[i]时会产生beauty relative,且产生的数量为 (n-i+1)*(i-j),j为i之前最后一个t与s不同的位置。
设状态 f[i][j][k] 为在第i位,有j对beauty relative,且i之前最后一个与s不同的位置在k 的串的数量,然而,第三维是可以化简的,将状态变成这样:
f[i][j] :在 i 位置且 s[i] != t[i] 且有 j 对 beauty relative 的串的数量,则:
f[i][j] = sigma( f[ k ][ j-(i-k)*(n-i+1) ] ) * ( 'z'-s[i] ) + ( sum[j] * (s[i]-'a' ) )
sum[j] 为长度为i有j对beauty relative的串的数量
每次转移复杂度 i / (n-i+1),需枚举i,j,总复杂度 n*n*log(n)
其实ans = sum[K],稍微写繁了
/* Created Time: Wednesday, November 13, 2013 PM07:17:06 CST */
#include
#include
#include
using namespace std;
typedef long long lld;
const int MOD = 1e9+7;
const int N = 2222;
int n,K;
// sum[j] : 长度为i有j对beauty relative的串的数量
// tail[i] : 从i~n所有位置都小于等于s串的t串总数
lld sum[N],tail[N],f[N][N];
char s[N];
lld work() {
lld ret = 0;
tail[n+1] = 1;
for (int i = n; i >= 1; i --)
tail[i] = tail[i+1]*(s[i]-'a'+1)%MOD;
if (K==0) return tail[1];
f[0][0] = 1;
sum[0] = 1;
for (int i = 1; i <= n; i ++) {
int dt = n-i+1;
for (int j = 0; j <= K; j ++) {
for (int k = i-1; k>=0 && (i-k)*dt<=j; k--) {
f[i][j] = (f[i][j]+f[k][j-(i-k)*dt])%MOD;
}
f[i][j] = f[i][j]*('z'-s[i])%MOD;
}
ret = (ret+f[i][K]*tail[i+1]%MOD)%MOD;
for (int j = 0; j <= K; j ++) {
f[i][j] = (f[i][j]+sum[j]*(s[i]-'a')%MOD)%MOD;
sum[j] = (sum[j]+f[i][j])%MOD;
}
}
return ret;
}
int main() {
scanf("%d%d%s",&n,&K,s+1);
printf("%lld\n",work());
return 0;
}