第一个最长的回文子串

回文串

输入一个字符串，求出其中第一个最长的回文子串 。
子串的含义是：在原串中连续出现的字符串片段。回文的含义是：正着看和倒着看都相同，如abba(偶)和yyxyy(奇)。
在判断时，应该忽略所有标点符号和空格，且忽略大小写，但输出应保持原样(在回文串的首部和尾部不要输出多余字符)。输入字符串长度不超过5000，且占据单独一行。应该输出最长的回文串，如果有多个，输出起始位置最靠左的。

样例输入：
Hello, Madam, I'm Adam.
样例输出：
Madam, I'm Adam

[常规的方法]

这种方法的效率太低，当目标字符串有5000个a时，找到最长回文子串的时间约为109.17s。

#include <cstdio> #include <cstring> #include <cctype> #include <ctime> #define MAXN 5000+10 #define LOCAL char buf[MAXN], s[MAXN]; int p[MAXN]; int main() { #ifdef LOCAL // redirection IO freopen("input.txt","r",stdin); freopen("output.txt","w",stdout); #endif int n, m=0, max=0, x, y, i, j, k; clock_t begin,end; fgets(buf,sizeof(s),stdin);// do not use gets() // check if (!strlen(buf)) { printf("buffer is null/n"); return 0; } n=strlen(buf); for (i=0; i<n; ++i) { if (isalpha(buf[i])) { p[m]=i; s[m++]=tolower(buf[i]); } } // check printf("%s/n",s); begin=clock(); for (i=0; i<m; ++i) { for (j=i; j<m; ++j) { int ok=1; for (k=i; k<=j; ++k) { if (s[k]!=s[i+j-k]) { ok=0; } } if (ok && j-i+1>max) { max=j-i+1; x=p[i]; y=p[j]; } } } end=clock(); // print for (i=x; i<=y; ++i) { printf("%c",buf[i]); } printf("/n"); printf("Time used = %.2lfs/n",(double)(end-begin)/CLOCKS_PER_SEC); return 0; } /* input: abxba cbaabc abcd,e.fghi ?. ihg?fedc.b.a output: abxbacbaabcabcdefghiihgfedcba abcd,e.fghi ?. ihg?fedc.b.a Time used = 0.00s */

[一种比较高效的方法]

枚举回文串的“中间”位置i，然后不断向外扩展，直到有字符不同。
这种方法的效率较高，当目标字符串有5000个a时，找到最长回文子串的时间约为0.09s。
#include <cstdio> #include <cstring> #include <cctype> #include <ctime> #define MAXN 5000+10 #define LOCAL char buf[MAXN], s[MAXN]; int p[MAXN]; int main() { #ifdef LOCAL // redirection IO freopen("input.txt","r",stdin); freopen("output.txt","w",stdout); #endif int n, m=0, max=0, x, y, i, j; clock_t begin,end; fgets(buf,sizeof(s),stdin);// do not use gets() // check if (!strlen(buf)) { printf("buffer is null/n"); return 0; } n=strlen(buf); for (i=0; i<n; ++i) { if (isalpha(buf[i])) { p[m]=i; s[m++]=tolower(buf[i]); } } // check printf("%s/n",s); begin=clock(); for (i=0; i<m; ++i) { // odd string for (j=0; i-j>=0 && i+j<m; ++j) { if (s[i-j]!=s[i+j]) { break; } if (j*2+1>max) { max=j*2+1; x=p[i-j]; y=p[i+j]; } } // even string for (j=0; i-j>=0 && i+j+1<m; ++j) { if (s[i-j]!=s[i+j+1]) { break; } if (j*2+2>max) { max=j*2+2; x=p[i-j]; y=p[i+j+1]; } } } end=clock(); // print for (i=x; i<=y; ++i) { printf("%c",buf[i]); } printf("/n"); printf("Time used = %.2lfs/n",(double)(end-begin)/CLOCKS_PER_SEC); return 0; } /* input: abxba cbaabc abcd,e.fghi ?. ihg?fedc.b.a output: abxbacbaabcabcdefghiihgfedcba abcd,e.fghi ?. ihg?fedc.b.a Time used = 0.00s */