Find Black Hand__HDU_4271__DP

题目地址：http://acm.hdu.edu.cn/showproblem.php?pid=4271

/*
Author:Bob Lee
2012.9.11
=======================================
感谢大牛：http://www.acforfun.com/?p=84
本程序核心思想来自于此博文
代码为自己所写
首先每个子串独立处理；对于每个子串，忽略母串是个环这个事实，
先考虑是个链的情况。
显然有个方程dp[i][j]表示母串前i位于子串前j位匹配且母串结尾没有多余字符，
最少需要的处理次数。
1.s1[i]==s2[j]，dp[i][j]=min{dp[i-1][j-1],dp[i-1][j]+1,dp[i][j-1]+1};
2.s1[i]!=s2[j],dp[i][j]=min{dp[i-1][j-1]+1,dp[i-1][j]+1,dp[i][j-1]+1};

下面考虑环的问题。
首先想到将母串s复制一遍粘贴到后面，然后dp，这是对的，但是不完整。
考虑这样一个母串abcd,子串abcda，最少需要操作次数是1，
但是将abcd复制后变成abcdabcd，发现没？abcda成了子串，答案变成了0。
也就是一个位置的字符被用了两次。分析一下会发现当这种情况发生时，
相当于跨越了长度大于len(母串长)的一个范围，而由于答案不大于子串长度L[i],
即不可能由大于L[i]*2的母串的substring生成，亦即当len>L[i]*2时，
此情况不会发生，此时只需要复制前L[i]*2个字符粘贴到母串后即可。

详细细节见代码：
*/


#include<iostream>
#include<cstring>
#include<cstdio>
#include<cstdlib>
using namespace std;
#define MAXN 100005
#define INF 99999999

int dp[MAXN*2][11];
char goal_s[MAXN];
char tmp_s[MAXN*2];
char small_s[15][15];

int min(int a,int b)
{
    return a<b?a:b;
}

int slove(char *s1,int len1,char *s2,int len2)
{
    //这里不用memset，很费时间，自己写的快一些
    for(int i=0;i<=len1;i++)
        for(int j=9;j<=len2;j++)
            dp[i][j] = 0;

    for(int j=1;j<=len2;j++)
    {
        dp[0][j] = j;//要十分注意这里的初始化
        //这个初始化的意思就是表示母串从空与目标串S2前j位匹配要花的最少步数位j
        //若不加上这个初始化就是个错误的程序
        for(int i=1;i<=len1;i++)
        {
            if(s1[i-1] == s2[j-1])
            {
                dp[i][j] = min(dp[i-1][j-1],min(dp[i][j-1]+1,dp[i-1][j]+1));
            }
            else
            {
                dp[i][j] = min(dp[i-1][j-1]+1,min(dp[i][j-1]+1,dp[i-1][j]+1));
            }
        }
    }

    int ans = INF;

    //找到最小的那一个返回
    for(int i=1;i<=len1;i++)
    {
        if(ans>dp[i][len2])
            ans = dp[i][len2];
    }

    return ans;
}

int main()
{
    while(scanf("%s",&goal_s) != EOF)
    {
        int n;
        scanf("%d",&n);
        int len1 = strlen(goal_s);
        //先把整个大的字符串，复制一遍到这个字符串的后面，留着在后面用
        for(int i=len1;i<len1+len1;i++)
            goal_s[i] = goal_s[i-len1];
        int ans = INF;
        int tmp;
        int anss;
        for(int i=1;i<=n;i++)
        {
            scanf("%s",&small_s[i]);
            int len2 = strlen(small_s[i]);

            /*
            如果说目标的子串的长度的2倍已经大于了母串的长度的话，
            那么我们就枚举在2倍母串中的长度小于等于len1的串，用来求目标子串
            最后求出这个里面的最小值，就是我们要求的
            */
            if(len2+len2 > len1)
            {
                int ii,jj,kk;
                for(ii=0;ii<len1;ii++)
                {
                    for(jj=1;jj<=len1;jj++)
                    {
                        for(kk=0;kk<jj;kk++)
                        {
                            tmp_s[kk] = goal_s[ii+kk];
                        }
                        tmp = slove(tmp_s,jj,small_s[i],len2);
                        if(ans > tmp || (ans == tmp && strcmp(small_s[anss],small_s[i])>0))
                        {
                            ans = tmp;
                            anss = i;
                        }
                    }
                }
            }
            else
            {
                tmp = slove(goal_s,len1+2*len2,small_s[i],len2);
                if(ans > tmp || (ans == tmp && strcmp(small_s[anss],small_s[i])>0))
                {
                    ans = tmp;
                    anss = i;
                }
            }
        }
        printf("%s %d\n",small_s[anss],ans);
    }
}