P3805 【模版】manacher演算法

題目描述

給出一個只由小寫英文字元a,b,c...y,z組成的字元串S,求S中最長迴文串的長度.

字元串長度為n

輸入輸出格式

輸入格式：

一行小寫英文字元a,b,c...y,z組成的字元串S

輸出格式：

一個整數表示答案

輸入輸出樣例

aaa

3

說明

字元串長度len <= 11000000

老呂教的manacher太low，，

寫一個T一個，

以後改寫位運算型的了。

一個點才300ms

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<cmath>
 5 #include<queue>
 6 #include<algorithm>
 7 #include<cstdlib>
 8 #define lli long long int 
 9 using namespace std;
10 const lli MAXN=31000011;
11 char s[MAXN];
12 char str[MAXN];
13 int ans[MAXN];
14 int len=0;
15 int ls=0;
16 void getstr()
17 {
18     str[0]='#';
19     str[1]='#';
20     for(int i=0;i<ls;i++)
21     str[(i<<1)+2]=s[i],str[(i<<1)+3]='#';
22     ls=(ls<<1)+2;
23     str[ls]=0;
24 }
25 void manacher()
26 {
27     getstr();
28     int mx=0,id=0;
29     
30     len=strlen(str);
31     for(int i=1;i<len;i++)
32     {
33         if(mx>i)
34             ans[i]=min(ans[2*id-i],mx-i);
35         else ans[i]=1;
36         while(str[i+ans[i]]==str[i-ans[i]])
37             ++ans[i];
38         if(i+ans[i]>mx)
39             mx=i+ans[i],id=i;
40             
41     }
42 }
43 int main()
44 {
45     scanf("%s",s);
46     ls=strlen(s);
47     manacher();
48     int out=0;
49     for(int i=0;i<len;i++)
50         out=max(out,ans[i]);
51     printf("%d",out-1);
52     return 0;
53 }