(advanced) 动态规划 UVA 1437 - String painter

There are two strings A and B with equal length. Both strings are made up of lower case letters. Now you have a powerful string painter. With the help of the painter, you can change a segment of characters of a string to any other character you want. That is, after using the painter, the segment is made up of only one kind of character. Now your task is to change A to B using string painter. What's the minimum number of operations?

Input 

Input contains multiple cases. Each case consists of two lines:

• The first line contains string A.
• The second line contains string B.

The length of both strings will not be greater than 100.

Output 

A single line contains one integer representing the answer.

Sample Input 

zzzzzfzzzzz 
abcdefedcba 
abababababab 
cdcdcdcdcdcd

Sample Output 

6 
7

题意:给出两个同样长度的串A,B，你要经过一些操作将A串变成B串，你的操作是在A串中任意选一段连续的片段，把他们变成同一个字母。问最少要多少步才能完成目标。

思路：注意到假设对于一个片段[i,j],我们可以由[i,k] [k+1,j] 合并过来，如果B[i]==B[j] [i,j]可以由[i+1,j],[i,j-1]中的更小的转移过来。这个不需要考虑A和B的关系，因为如果我们把[i,j]操作了，那么其实这个时候和原来的A就没有任何关系了，你的目标只是用尽可能少的步数将[i,j]中的字母变成和B中对应的[i,j]的字母,而这个我们用一个g[i][j]表示更改了[i,j]区域后至少需要啊多少步才能达到目标。在处理这个g[i][j]的时候需要三重循环，分别是头、尾和拆开的地方。为了然后状态时从“小”往“大”的推，我们可以每次处理的片段是先处理完长度为3的，然后由长度是3的推出长度是4的。。。以此类推。长度为1和长度为2的可以一开始预处理。

那么还有条件没用上，A和B的相等的关系怎么处理呢？我们用一个f[n]表示前n个字母最少要多少次操作才能变成一样。那么对于一个字母，如果A[i]==B[i]，那么可以有f[i]=min(f[i-1],f[k]+g[k+1][i]) 0<=k<i ,如果不等的话，f[i]=min(f[k]+g[k+1][i])

那么最后f[n]就是要求的

启示：有时候dp只用一个数组表示状态可能是不够的，这时候不要把自己的思维只停留在一个状态表示，可能需要多个状态并用。

代码：

#include<iostream>
#include<cstdio>
#include<string.h>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=100+5;

int f[maxn];
int g[maxn][maxn];
int n;
char A[maxn],B[maxn];

void solve()
{
    memset(f,0x3f,sizeof(f));
    memset(g,0x3f,sizeof(g));
    f[0]=0;
    for(int i=1;i<=n;++i)
    {
        g[i][i]=1;
        if(i+1<=n)
        {
            if(B[i]==B[i+1]) g[i][i+1]=1;
            else g[i][i+1]=2;
        }
    }
    for(int k=3;k<=n;++k)
    {
        for(int i=1;i+k-1<=n;++i)
        {
            int j=i+k-1;
            for(int q=i;q<i+k-1;++q)
                g[i][j]=min(g[i][j],g[i][q]+g[q+1][j]);
            if(B[i]==B[j]) g[i][j]=min(g[i][j],min(g[i+1][j],g[i][j-1]));
           // printf("%d %d: %d\n",i,j,g[i][j]);
        }
    }
    for(int i=1;i<=n;++i)
    for(int j=0;j<i;++j)
    {
        if(A[i]==B[i]) f[i]=f[i-1];
        f[i]=min(f[i],f[j]+g[j+1][i]);
    }
    printf("%d\n",f[n]);
}

int main()
{
    while(scanf("%s%s",A+1,B+1)==2)
    {
        n=strlen(A+1);
        solve();
    }
}