【精】POJ 3267 The Cow Lexicon（区间动归）

poj 3267

• 题目大意

给你一个字母序列L(2 ≤ L ≤ 300)，以及W个单词 (1 ≤ W ≤ 600) ，问最少去掉多少字母使得剩下的序列都是由单词组成的。

• 分析
  这道题是在做拓扑排序的时候出现的，看了之后怎么都感觉和拓扑没太大关系，于是就没管拓扑了。这种题看着就像动归，算是动态规划的基本题型吧。

• 动态规划的几个特点：

  • 相同子问题
  • 最优子结构
  • 无后效性

这里就借这道题写一下分析动态规划的一般思路吧。分析问题时的一个好的习惯是从简单到复杂，不要一开始就希望想出多么高效的算法（不过熟之后确实可以这样，大神请无视）。拿动态规划来说，动态规划的高效体现在具有相同子问题上，动态规划通过把相同的子问题记录下来从而避免了重复的对相同问题的求解。我觉得很多问题都可以从简单粗暴的算法入手，比如枚举暴力搜索，在用枚举暴力分析是慢慢就会发现有些地方是可以通过贪心动归分治等思想来优化的，从而提升了算法的效率。
比如这道题中，枚举单词的序列肯定是不行的，但我们可能会想到搜索的方法，遍历一遍单词，对每个单词，找出第一次能匹配上的位置或不能匹配，然后在剩余的子串中再进行所有单词的匹配。仔细思考之后就会发现有很多相同的问题我们是重复求解了的，比如在第i个和第j个单词作为第一个单词匹配上之后剩下的子串长度相同，这就是相同子问题，我们可以只求解一次这样的问题并且将结果保存下来以后遇到就直接返回结果不用重复求解一次。这也叫做记忆化搜索，用空间来换取时间。
现在给出问题的更严格的定义：
令dp[i]表示从下标为i的字母开始到字母序列L结尾需要去掉的字母数，状态确定之后就是状态转移方程的事的，我认为动态规划最重要的就是状态的确定，状态一旦确定了，转移方程是相对比较容易的。转移方程：
—-dp[i]=min{dp[i], dp[j]+j-i-len_word[k] | 第k个单词能在i开头的字符串中匹配上}
—-dp[i]=L-i; 所有单词都不能匹配上，也就是需要把所有字母都去掉（字符串下标从0开始）
可以用递归搜索（记忆化搜索）的写法写也可以用递推的方式写。个人喜欢用递归，但这样有时候效率会比较低。

• 代码

#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstring>
#include<cstdlib>
#include<queue>
using namespace std;
const int INF=99999999;
int W,L;
char  s[400];
char  word[700][30];
int dp[700];//dp[i]表示以字符串中编号为i的字母开始的字符串需要去掉多少字母
int Dp(int k)
{
      int ans=INF;
      if(dp[k]>=0)return dp[k];
      if(k>L-1)return dp[k]=0;//边界条件
      if(k==L-1)return dp[k]=1;
      int i,j;
      int len_word;
      for(int t=1;t<=W;t++)//分别判断第t个单词能否和k开头的字符串匹配上
      {
            len_word=strlen(word[t]);
            if(len_word>(L-k))continue;
            i=0;
            j=k;
            while(j<L)
            {
                  if(word[t][i]==s[j]){j++;i++;}
                  else j++;
                  if(i==len_word){ans=min(ans,Dp(j)+j-k-len_word);break;}
            }
      }
      if(ans==INF)return dp[k]=L-k;
      else return dp[k]=ans;
}
int main()
{
   while(scanf("%d%d",&W,&L)!=EOF)
   {
        scanf("%s",&s);
        for(int i=1;i<=W;i++)scanf("%s",&word[i]);
        memset(dp,-1,sizeof(dp));
        cout<<Dp(0)<<endl;
   }
}