Manacher算法???

  Mancher算法是用来找最大回文串的,因为功能太过单一,所以用处不是很多,但是在找回文串方面还是很强大的。

大致原理:

  枚举回文串的中点,用p[i]表示以i为中心的最大回文串的长度。然后对p[i]进行扩展就可以了。

    最简单的扩展就是当s[i+p[i]+1]=s[i-p[i]-1]的时候p[i]就可以加一。

    然后就是对上面进行一下优化,用一个id来记录可以延伸的最靠右的回文串的中心,mx为其长度。然后对于枚举到的中心i如果i没有超过id+mx,说明在id前面一定有一个j与i相同,并且j已经处理过了,所以就可以直接将p[i]=p[j]但是需要注意的是,如果i+p[j]超出了mx,就不能保证这是个回文串了,p[i]=min(p[id-(i-id)],mx-i)    这里的id-(i-id)就是上面的j

    同时不断地更新id和mx并且记录下最大的p[i]就行了。

一道板子题luogu3805

代码:

 1 #include<vector>
 2 #include<queue>
 3 #include<ctime>
 4 #include<algorithm>
 5 #include<deque>
 6 #include<cstdlib>
 7 #include<string>
 8 #include<cstring>
 9 #include<map>
10 #include<bitset>
11 #include<cstdio>
12 #include<iostream>
13 #define max(x,y) x>y?x:y
14 using namespace std;
15 typedef long long LL;
16 void read(int &x){
17     x=0;int f=1;char c=getchar();
18     while(c<'0'||c>'9'){if(c=='-') f=-1;c=getchar();}
19     while(c>='0'&&c<='9') { x=x*10+c-'a';c=getchar();}x*=f;
20 }
21 const int N=11000000+10;
22 char s[N],S[N*2];
23 int ans,len,id,mx,p[N*2];
24 inline void manacher()
25 {
26     for(int i=1;i<=len;++i)
27     {
28         if(id+mx>i)
29         p[i]=min(id+mx-i,p[id*2-i]);
30         while(i+p[i]+1<=len&&i-p[i]-1>=1&&S[i+p[i]+1]==S[i-p[i]-1]) p[i]++;
31         if(i+p[i]>id+mx) id=i,mx=p[i];
32         ans=max(ans,p[i]);
33     }
34     printf("%d",ans);
35 }
36 
37 int main()
38 {
39     cin>>s;
40     int l=strlen(s);
41     S[++len]='#';
42     for(int i=0;i<l;++i)
43     {
44         S[++len]=s[i];
45         S[++len]='#';
46     }
47     manacher();
48     
49     return 0;
50 }