【BZOJ】2342 [Shoi2011]双倍回文 manacher+set

来源:互联网 发布:three.js 案例 编辑:程序博客网 时间:2024/06/06 07:19

题目传送门

“马拉车”什么的好暴力啊……总感觉就是暴力做法,但是好像时间复杂度是O(len)的,比较强啊……

对于这道题,我还是too young了,没想到w和wR是可以重叠的……

首先我们要明白manacher算法,这个看看大佬的博客就行了:大佬传送门

然后这题就变的比较水了,首先求出每两个字符间的回文半径p。

然后考虑答案中的w wR w wR的第二个空格处,对于第三个字符串w,设它的左端点为x,右端点为y。

若当前的双倍回文满足成为答案的条件,y-p[y]<=x且x+a[x]/2>=y。

因此我们可以对所有的位置按i-p[i]排序,然后按顺序把所有位置插入set中,每次求出最大的x满足i+a[i]/2>=x并更新答案即可。

附上AC代码:

#include <cstdio>#include <algorithm>#include <set>using namespace std;const int N=500010;struct note{int wz,w;bool operator < (const note lyf) const {return w<lyf.w;}}f[N];set <int> st;int n,pos,len,a[N<<1],w[N],h,ans;char t[N],s[N<<1];int main(void){scanf("%d%s",&n,t+1);for (int i=1; i<=n; ++i) s[i*2-1]='%',s[i*2]=t[i];s[len=n*2+1]='%',pos=1;for (int i=1; i<=len; ++i){if (i<pos+a[pos]) a[i]=min(a[2*pos-i],pos+a[pos]-i); else a[i]=1;while (i-a[i]>0&&i+a[i]<len+1&&s[i-a[i]]==s[i+a[i]]) ++a[i];if (i+a[i]>pos+a[pos]) pos=i;}for (int i=1; i<=n; ++i) w[i]=(a[i<<1|1]-1)>>1;for (int i=1; i<=n; ++i) f[i]=(note){i,i-w[i]};sort(f+1,f+1+n),h=0;for (int i=1; i<=n; ++i){while (h<n&&f[h+1].w<=i) st.insert(f[++h].wz);if ((*st.lower_bound(i+w[i]/2))==*st.begin()) continue;ans=max(ans,(*st.lower_bound(i+w[i]/2)-i));}printf("%d",ans<<2);return 0;}