BZOJ 2342 [Shoi2011]双倍回文

---

【题目分析】
qy神犇讲的是用二维区间查询的方式来更新答案，需要用到扫描线。

p[i]表示i和i+1为中心的最长回文子串长度/2（str[i-k]=str[i+1+k]）。。。
用manacher On计算p数组
题目要求计算w wR w wR的最大长度
枚举x为对称轴。。。实际上对称轴在x到x+1之间，即x是第一个wR的最后一位
举些例子推一推发现len(x+1,y)*4能更新答案，仅当y-p[y]<=x且y<=x+p[x]/2
按照y-p[y]排序一下，递推x的时候将符合1式的y插入set，在set中查找x+p[x]/2的前驱更新答案即可
复杂度On+nlogn
usedtobe提出如果把题目要求改成w wR w。。。这样的话限制应该是y<=x+p[x]且y-p[y]<=x

——以上内容来自hzw
感觉这种方法更有技巧，剪掉了许多不必要的更新，确实好。但是我也确实没有看懂。。。

---

【代码】

#include<set>#include<cstdio>#include<cstdlib>#include<algorithm>using namespace std;int n,ans;char ch[500005];int p[500005],q[500005];set<int> t;void manacher(){    int 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(;ch[i+p[i]+1]==ch[i-p[i]];p[i]++);        if(p[i]+i>mx) id=i,mx=p[i]+i;    }}bool cmp(int a,int b){return (a-p[a])<(b-p[b]);}int main(){    scanf("%d",&n);    scanf("%s",ch+1);    ch[0]='#';    manacher();    for(int i=1;i<=n;i++)q[i]=i;    sort(q+1,q+n+1,cmp);    int now=1;    for(int i=1;i<=n;i++)    {        while(now<=n&&q[now]-p[q[now]]<=i)  t.insert(q[now]),now++;        set<int>::iterator tmp=t.upper_bound(i+p[i]/2);        if(tmp!=t.begin())  ans=max(ans,(*--tmp-i)*4);    }    printf("%d\n",ans);    return 0;}