1. 程式人生 > >【洛谷4287】[SHOI2011] 雙倍迴文(Manacher演算法經典題)

【洛谷4287】[SHOI2011] 雙倍迴文(Manacher演算法經典題)

點此看題面

大致題意: 求一個字串中有多少個長度為偶數的迴文串,它的一半也是迴文串。

ManacherManacher演算法

這應該是ManacherManacher演算法一道比較好的入門題,強烈建議在做這題之前先去學一學ManacherManacher演算法。

LinkLink

ManacherManacher演算法 詳見部落格 Manacher演算法學習筆記

大致思路

如果你會ManacherManacher演算法,其實這題的思路應該還是比較清晰的。

只要在更新最大右邊界的同時,判斷每一個新出現的迴文串是否為雙倍迴文串

,並更新ansans即可。

程式碼

#include<bits/stdc++.h>
#define max(x,y) ((x)>(y)?(x):(y))
#define min(x,y) ((x)<(y)?(x):(y))
#define uint unsigned int
#define LL long long
#define ull unsigned long long
#define swap(x,y) (x^=y,y^=x,x^=y)
#define abs(x) ((x)<0?-(x):(x))
#define INF 1e9
#define Inc(x,y) ((x+=(y))>=MOD&&(x-=MOD))
#define ten(x) (((x)<<3)+((x)<<1)) #define N 500000 using namespace std; int n;string st; class FIO { private: #define Fsize 100000 #define tc() (FinNow==FinEnd&&(FinEnd=(FinNow=Fin)+fread(Fin,1,Fsize,stdin),FinNow==FinEnd)?EOF:*FinNow++) #define pc(ch) (FoutSize<Fsize?Fout[FoutSize++]=ch:(fwrite(Fout,1,FoutSize,stdout),Fout[(FoutSize=0)++]=ch))
int f,FoutSize,OutputTop;char ch,Fin[Fsize],*FinNow,*FinEnd,Fout[Fsize],OutputStack[Fsize]; public: FIO() {FinNow=FinEnd=Fin;} inline void read(int &x) {x=0,f=1;while(!isdigit(ch=tc())) f=ch^'-'?1:-1;while(x=ten(x)+(ch&15),isdigit(ch=tc()));x*=f;} inline void read_char(char &x) {while(isspace(x=tc()));} inline void read_string(string &x) {x="";while(isspace(ch=tc()));while(x+=ch,!isspace(ch=tc())) if(!~ch) return;} inline void write(int x) {if(!x) return (void)pc('0');if(x<0) pc('-'),x=-x;while(x) OutputStack[++OutputTop]=x%10+48,x/=10;while(OutputTop) pc(OutputStack[OutputTop]),--OutputTop;} inline void write_char(char x) {pc(x);} inline void write_string(string x) {register int i,len=x.length();for(i=0;i<len;++i) pc(x[i]);} inline void end() {fwrite(Fout,1,FoutSize,stdout);} }F; class Class_Manacher//Manacher演算法 { private: int len,p[(N<<1)+5]; char ns[(N<<1)+5]; inline void Init(string st) { register int i,l=st.length(); for(ns[i=0]='!',ns[len=1]='%';i<l;++i) ns[++len]=st[i],ns[++len]='%'; ns[++len]='?'; } public: inline int GetAns(string st)//基本上是模板,只是在更新最大右邊界的同時多了一個更新ans的過程 { register int i,j,ans=0,id,Max=0; for(Init(st),i=1;i<len;++i) { p[i]=i<=Max?min(p[(id<<1)-i],p[id]+id-i):1; while(!(ns[i-p[i]]^ns[i+p[i]])) ++p[i]; if(i+p[i]>Max) { if(i&1) for(j=i+(Max>i+4?((Max-i+3>>2)<<2):4);j<i+p[i];j+=4) if(p[i-(j-i>>1)]>(j-i>>1)) ans=max(ans,j-i);//首先判斷長度是否為偶數,然後判斷每一個新出現的迴文串是否為雙倍迴文串,並更新ans Max=i+p[id=i]; } } return ans; } }Manacher; int main() { F.read(n),F.read_string(st),F.write(Manacher.GetAns(st)); return F.end(),0; }