hdu2087 KMP+简单动归

裸KMP+简单动态规划

先用KMP求出所有匹配位置的尾节点，标记为1，然后对标记数组进行动态规划，每段成功匹配的位置有两种选择，要么取要么不取，dp[i]=max(dp[i-1],dp[i-m]+1)，m为匹配串的长度。

代码：

#include <iostream>
#include <cstring>
#include <cstdio>
#include <string>

using namespace std;

int dp[1001];
char s1[1001];
char s2[1001];
int p[1001];
int mark[1001];
int n,m;

void getnext()
{
    for (int i=2;i<=m;i++)
    {
        int j=p[i-1];
        while (j>0&&s2[i]!=s2[j+1])
            j=p[j];
        if (s2[i]==s2[j+1])
            j++;
        p[i]=j;
    }
}

int max(int a,int b)
{
    return a>b?a:b;
}

int main()
{
    while (scanf("%s",s1+1))
    {
        if (s1[1]=='#')
            break;
        scanf("%s",s2+1);
        memset(mark,0,sizeof(mark));
        memset(dp,0,sizeof(dp));
        p[0]=p[1]=0;
        getnext();
        n=strlen(s1+1);
        m=strlen(s2+1);
        int j=1;
        for (int i=1;i<=n;i++,j++)
        {
            if (j==m&&s2[j]==s1[i])
            {
                mark[i]=1;
                i=i-j+1;
                j=0;
               //cout<<i<<" "<<j<<endl;
                continue;
            }
            if (s1[i]==s2[j])
                continue;
            j--;
            while (j>0&&s2[j+1]!=s1[i])
                j=p[j];
            if (j==0)
                j--;
            j++;
        }
        for (int i=1;i<=n;i++)
        {
            dp[i]=dp[i-1];
            if (mark[i]==1)
                dp[i]=max(dp[i],dp[i-m]+1);
            //cout<<i<<" "<<dp[i]<<endl;
        }
        cout<<dp[n]<<endl;
    }
}