码迷,mamicode.com
首页 > 其他好文 > 详细

bzoj2342

时间:2017-09-19 16:45:37      阅读:173      评论:0      收藏:0      [点我收藏+]

标签:style   数据结构   std   can   答案   max   alt   情况   logs   

manacher+set

好像这种求回文串形态的都是用manacher,然后求出中心之间奇怪的关系

我们先跑出来manacher,由于只要偶数回文串,我们就只把#位置记录下来,重新保存,然后我们考虑什么情况下区间[x+1,y]*4可以更新答案,x是整个回文串的中心。那么很明显有x+f[x]/2>=y,y-p[y]>=x

碰见这种情况,我们一般先保证一个不等式恒成立,另一个用数据结构搞一搞,那么我们让y-p[y]>=x恒成立,先对y排个序,然后将满足条件的y逐渐加入,插入set,我们枚举的是x,自然希望y越大越好,那么我们在set里查x+f[x]/2的前继,看能是否更新答案。

技术分享
#include<bits/stdc++.h>
using namespace std;
const int N = 500010;
int n, len, ans, pos, mx;
char t[N], s[N << 1];
int p[N << 1], f[N], y[N];
set<int> S;
bool cmp(int i, int j) { return i - f[i] < j - f[j]; }
int main()
{
    scanf("%d%s", &n, t + 1);
    s[0] = !;
    s[len = 1] = #;
    for(int i = 1; i <= n; ++i) s[++len] = t[i], s[++len] = #; 
    for(int i = 1; i <= len; ++i) 
    {
        if(i < mx) p[i] = min(mx - i, p[2 * pos - i]); else p[i] = 1;
        while(s[i - p[i]] == s[i + p[i]]) ++p[i];
        if(i + p[i] > mx)
        {
            pos = i;
            mx = i + p[i];
        }
    }
    for(int i = 1; i <= n; ++i) f[i] = (p[i << 1 | 1] - 1) >> 1, y[i] = i;
    sort(y + 1, y + n + 1, cmp);
    int j = 0;
    for(int i = 1; i <= n; ++i)
    {
        while(y[j + 1] - f[y[j + 1]] <= i && j < n) ++j, S.insert(y[j]);
        set<int> :: iterator it = S.upper_bound(i + (f[i] >> 1));
        if(it == S.begin()) continue;
        else --it;
        if(i + (f[i] >> 1) >= *it && *it - f[*it] <= i) ans = max(ans, (*it - i) * 4);
    }
    printf("%d\n", ans);
    return 0;
}
View Code

 

bzoj2342

标签:style   数据结构   std   can   答案   max   alt   情况   logs   

原文地址:http://www.cnblogs.com/19992147orz/p/7552497.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!