Codeforces 360C Levko and Strings (dp)

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;
}

     
     
    
    
   
   

  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
区间DP是一种动态规划的方法,用于解决区间范围内的问题。在Codeforces竞赛中,区间DP经常被用于解决一些复杂的字符串或序列相关的问题。 在区间DP中,dp[i][j]表示第一个序列前i个元素和第二个序列前j个元素的最优解。具体的转移方程会根据具体的问题而变化,但是通常会涉及到比较两个序列的元素是否相等,然后根据不同的情况进行状态转移。 对于区间长度为1的情况,可以先进行初始化,然后再通过枚举区间长度和区间左端点,计算出dp[i][j]的值。 以下是一个示例代码,展示了如何使用区间DP来解决一个字符串匹配的问题: #include <cstdio> #include <cstring> #include <string> #include <iostream> #include <algorithm> using namespace std; const int maxn=510; const int inf=0x3f3f3f3f; int n,dp[maxn][maxn]; char s[maxn]; int main() { scanf("%d", &n); scanf("%s", s + 1); for(int i = 1; i <= n; i++) dp[i][i] = 1; for(int i = 1; i <= n; i++) { if(s[i] == s[i - 1]) dp[i][i - 1] = 1; else dp[i][i - 1] = 2; } for(int len = 3; len <= n; len++) { int r; for(int l = 1; l + len - 1 <= n; l++) { r = l + len - 1; dp[l][r] = inf; if(s[l] == s[r]) dp[l][r] = min(dp[l + 1][r], dp[l][r - 1]); else { for(int k = l; k <= r; k++) { dp[l][r] = min(dp[l][r], dp[l][k] + dp[k + 1][r]); } } } } printf("%d\n", dp[n]); return 0; } 希望这个例子能帮助你理解区间DP的基本思想和应用方法。如果你还有其他问题,请随时提问。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值