动态规划 （五）

听说周末就要开始学新的知识了，感觉最近十分的紧张，要做的题还有很多，还需要再看一些例题。上周学的是区间dp问题，基本概念就是把适合通过转化为区间dp的题目，将这一个大区间分别变为几个小区间，分别去求小区间的最优解，合成为大区间的最优解。

首先是第一种模型，将n堆石子合成一堆的最小花费。
我们首先先简单看到初始的几个值
sum[1]=0.sum[2]=2，sum[3]=min(a[1]+a[2],a[2]+a[3])+sum[3].
动态规划的方程式：
dp[i][j] = min(dp[i][k] + dp[k+1][j]) + sum[j] - sum[i-1]
我觉得动态规划有两个难点，一个是动态转移方程的建立，还有就是循环的手段。
这道题循环的手段是
for (int len = 2; len <= n; len++)
{
for (int i = 1, j = len; j <= n; i++, j++)
{
for (int k = i; k < j; k++)
{
if(dp[i][j] > dp[i][k] + dp[k+1][j] + sum[j] - sum[i-1])
dp[i][j] = dp[i][k] + dp[k+1][j] + sum[j] - sum[i-1];
}
}
len可以变相理解为j的起点，i是全方程的起点，k是方程中循环的序。
这道题就可以解出了。

第二种模型是括号匹配问题，
规定dp[i][j]为合并第i到第j个字符的最大匹配度
长度为n时，我们可以先检测a[i]和a[j]是否匹配， 匹配dp[i][j] = dp[i+1][j-1] + 2
不匹配，那我们可以按第一种模型处理，从任意位置分成两个区间
这道题的循环方式换第一种模型类似。

我们可以把[i,j]区间的字符当成由[i+1,j]在前面加个字符或[i,j-1]在后面加一个字符得来的.
这里我们只考虑[i,j]由[i+1,j]在前面加一个字符的情况
如果a[i+1]到a[j]没有和a[i]匹配的，
那么dp[i][j] = dp[i+1][j]
如果a[k]和a[i]匹配（i < k <= j），那么
dp[i][j] = max(dp[i][j], dp[i+1][k-1] + dp[k+1][j] + 2);
主函数代码如下：

while (gets(a+1))
{
    if(a[1] == 'e') break;
    memset(dp, 0, sizeof(dp));
    int n = strlen(a+1);
    for (int len = 2; len <= n; len++)
    {
        for(int i = 1, j = len; j <= n; i++, j++)
        {
            dp[i][j] = dp[i+1][j];
            for (int k = i; k <= j; k++)
                if((a[i]=='('&&a[k]==')') || (a[i]=='['&&a[k]==']'))
                    dp[i][j] = max(dp[i][j], dp[i+1][k-1] + dp[k+1][j] + 2);
        }
    }
    printf("%d\n",dp[1][n]);
}

第三种模型是只与左右边界有关的相关问题。（直接上例题了）
n个字符组成长度为m的字符串，给出增删字符的花费，可在字符串任意位置增删字符，求把字符串修改成回文串的最小花费。
规定dp[i][j]为将[i,j]区间改成回文串的最小花费，可以看成有回文串
[i+1,j]在前面加一个字符 =>前面删或后面增
[i,j-1]在后面加一个字符 =>前面增或后面删
共四种情况，当a[i] == a[j]时，加个dp[i+1][j-1]的情况就好了
大体代码如下：

const int MAXN = 2005;
char a[MAXN], ch;
int dp[MAXN][MAXN];
int add[30],sub[30];
int main()
{
    int n, m;
    scanf("%d %d", &n, &m);
    scanf("%s", a+1);
    for (int i = 1; i <= n; i++)
    {
        scanf(" %c", &ch);
        scanf("%d %d", &add[ch-'a'], &sub[ch-'a']);
    }
    memset(dp, 0x3f, sizeof(dp));
    for (int i = 1; i <= m; i++) 
        dp[i][i] = 0;
    for (int len = 2; len <= m; len++)
        for(int i = 1, j = len; j <= m; i++, j++)
        {
            dp[i][j] = min(dp[i][j], min(add[a[i]-'a'],sub[a[i]-'a']) + dp[i+1][j]);
            dp[i][j] = min(dp[i][j], dp[i][j-1] + min(add[a[j]-'a'],sub[a[j]-'a']));
            if (a[i] == a[j])
            {
                if (len==2) 
                    dp[i][j] = 0;
                else
                    dp[i][j] = min(dp[i][j], dp[i+1][j-1]);
            }
        }
    printf("%d\n", dp[1][m]);
    return 0;
}

补充：所有的dp[i][i]由于代表是两个数之间，相同之间没有数字，直接归0