POJ-3974 Palindrome Manacher算法

　　题目链接：http://poj.org/problem?id=3974

　　Manacher算法效率真不错，用后缀数组A的都沙茶了。。

 1 //STATUS:C++_AC_235MS_10904KB

 2 #include<stdio.h>

 3 #include<stdlib.h>

 4 #include<string.h>

 5 #include<math.h>

 6 #include<iostream>

 7 #include<string>

 8 #include<algorithm>

 9 #include<vector>

10 #include<queue>

11 #include<stack>

12 #include<map>

13 using namespace std;

14 #define LL __int64

15 #define pii pair<int,int>

16 #define Max(a,b) ((a)>(b)?(a):(b))

17 #define Min(a,b) ((a)<(b)?(a):(b))

18 #define mem(a,b) memset(a,b,sizeof(a))

19 #define lson l,mid,rt<<1

20 #define rson mid+1,r,rt<<1|1

21 const int N=1000010,INF=0x3f3f3f3f,MOD=10000,STA=8000010;

22 const double DNF=1e13;

23 

24 char str[N<<1],s[N];

25 int p[N<<1];

26 int n,len;

27 

28 void Manacher(char *str,int *p)

29 {

30     int i,j,id,mx;

31     id=1,mx=1;

32     p[0]=p[1]=1;

33     for(i=2;i<n;i++){

34         p[i]=1;

35         if(mx>i){

36             p[i]=Min(p[(id<<1)-i],mx-i);

37         }

38         while(str[i+p[i]]==str[i-p[i]])p[i]++;

39         if(i+p[i]>mx){

40             id=i;

41             mx=i+p[i];

42         }

43     }

44 }

45 

46 void getstr(char *s)

47 {

48     int i;

49     str[0]='$';str[1]='#';

50     for(i=0;i<len;i++){

51         str[(i<<1)+2]=s[i];

52         str[(i<<1)+3]='#';

53     }

54     str[n]=0;

55 }

56 

57 int main()

58 {

59  //   freopen("in.txt","r",stdin);

60     int i,j,ans,ca=1;

61     while(~scanf("%s",s))

62     {

63         len=strlen(s);

64         if(len==3 && s[0]=='E' && s[1]=='N' && s[2]=='D')break;

65         n=len*2+4;

66         getstr(s);

67         Manacher(str,p);

68 

69         ans=1;

70         for(i=2;i<n;i++){

71             ans=Max(ans,p[i]-1);

72         }

73 

74         printf("Case %d: %d\n",ca++,ans);

75     }

76     return 0;

77 }