近似回文串、最长回文子串

输入一行文本，输出最长近似回文词连续子串。所谓近似回文词是指满足以下条件的字符串：

1. S以字母开头，字母结尾

2. a(S)和b(S)最多有2k个位置不同，其中a(S)是S删除所有非字母字符并且把所有字母转化成小写之后得到的串，b(S)是a(S)的逆序串。

比如当k=1时，Race cat是一个近似回文词，因为a(S)=racecat和b(S)=tacecar只有2个位置不同。

Input
输入包含不超过25组数据，每组数据包含两行。第一行是整数k（0<=k<=200），第二行为字符串S，包含至少一个字母但不超过1000个字符（换行符不算）。S只包含字符、空格和其他可打印字符（比如逗号，句号），并且不会以空白字符开头。

Output
对于每组测试数据，输出最长近似回文子串的长度和起始位置（S的第一个字符是位置1）。如果有多个最长近似回文子串解，起始位置应尽量小。

Sample Input
1
Wow, it is a Race cat!
0
abcdefg
0
Kitty: Madam, I'm adam.
Sample Output
Case 1: 8 3
Case 2: 1 1

Case 3: 15 8

该题解题思路类似于求一个字符串的最长回文子串，只是他允许2*k个地方不同。

求最长回文子串的代码：

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

using namespace std;
const int N=1005;

int main()
{
    char str[N];
    while(~scanf("%s",str))
    {
        int len=strlen(str);
        int maxlen=0,start=0;
        for(int i=0;i<len;i++)
        {
            int j;
            for(j=0;j<=i&&j+i<len;j++)//判断ababa类型的回文串
            {
                if(str[i-j]!=str[i+j])
                    break;
            }
            --j;
            if(2*j+1>maxlen)
            {
                maxlen=2*j+1;
                start=i-j;
            }

            for(j=0;j<=i&&j+i+1<len;j++)//判断abaaba类型的回文串
            {
                if(str[i-j]!=str[i+j+1])
                    break;
            }
            --j;
            if(j<0)
                continue;
            if(2*j+2>maxlen)
            {
                maxlen=2*j+2;
                start=i-j;
            }

        }
        for(int i=start;i<start+maxlen;i++)
            printf("%c",str[i]);
        printf("\n");
        printf("%d\n",maxlen);
    }
    return 0;
}

求近似回文串的代码：

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

using namespace std;
const int N=1005;
int main()
{
    int pos[N];
    char str[N];
    char as[N],bs[N];
    int k,cas=1;
    while(~scanf("%d",&k))
    {
        getchar();
        memset(as,0,sizeof(as));
        memset(bs,0,sizeof(bs));
        memset(pos,0,sizeof(pos));
        gets(str);
        int len=strlen(str);
        int j=0;
        for(int i=0; i<len; i++)//得到纯字母的字符串
        {
            if(isalpha(str[i]))
            {
                pos[j]=i;//记录在原字符串中的位置
                as[j++]=tolower(str[i]);
            }

        }
        int maxlen=0,start=0,l;
        for(int i=0; i<j; i++)
        {
            int cnt=0;
            for(l=0; l<=i&&i+l<j; l++)//判断ababa类型
            {
                if(as[i-l]!=as[i+l])//记录不同的个数。
                    cnt++;
                if(cnt>k)
                {
                    break;
                }
            }
            --l;//退回到上一次相同的地方
            if(pos[i+l]-pos[i-l]+1>maxlen)//在原字符串中的距离大于maxlen
            {
                maxlen=pos[i+l]-pos[i-l]+1;
                start=pos[i-l];
            }
            cnt=0;
            for(l=0; l<=i&&i+l+1<j; l++)//判断abaaba类型
            {
                if(as[i-l]!=as[i+l+1])
                    cnt++;
                if(cnt>k)
                {
                    break;
                }
            }
            --l;
            if(l<0)
                continue;
            if(pos[i+l+1]-pos[i-l]+1>maxlen)
            {
                maxlen=pos[i+l+1]-pos[i-l]+1;
                start=pos[i-l];
            }

        }
        printf("Case %d: %d %d\n",cas++,maxlen,start+1);
    }
    return 0;
}