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;
}