BZOJ2342 双倍回文 [manacher][set]

来源:互联网 发布:angularjs repeat源码 编辑:程序博客网 时间:2024/05/21 07:30

手写字符串读优莫名RE,改成scanf就AC了,莫(m)队(d)指(z)针(z)。。。

枚举x为对称轴,即枚举第一个wR的最后一位。
找规律发现len(x+1,y)4能更新答案,仅当yp[y]<=xy<=x+p[x]/2

按照yp[y]排序一下,递推x的时候将符合1式的y插入set,在set中查找x+p[x]/2的前驱更新答案即可

#include<set>#include<cstdio>#include<cstring>#include<iostream>#include<algorithm>using namespace std;const int N=1001000;set<int> mp;set<int> :: iterator it;int n,p[N<<1],q[N<<1],flag=1;char s[N<<1],now[N<<1];void readin(int &res){    static char ch;register int flag=1;    while((ch=getchar())<'0'||ch>'9')if(ch=='-')flag=-1;res=ch-48;    while((ch=getchar())>='0'&&ch<='9')res=res*10+ch-48;res*=flag;}bool cmp(int a,int b){    return (a-p[a])<(b-p[b]);}int Manacher(){    int ans=0,mx=0,id;    for(int i=1;i<=n;i++){        if(mx>=i)p[i]=min(mx-i,p[2*id-i]);        else p[i]=0;        for(;s[i+p[i]+1]==s[i-p[i]];p[i]++);        if(p[i]+i>mx)id=i,mx=p[i]+i;        q[i]=i;    }    sort(q+1,q+1+n,cmp);    for(register int i=1;i<=n;i++){        while(flag<=n&&q[flag]-p[q[flag]]<=i)mp.insert(q[flag++]);        it=mp.upper_bound(i+p[i]/2);        if(it!=mp.begin())ans=max(ans,(*--it-i)*4);    }    return ans;}int main(){    readin(n),scanf("%s",s+1),cout<<Manacher()<<endl;    return 0;}

这里写图片描述

原创粉丝点击