codeforces #30E Tricky and Clever Password KMP+Manacher+二分

题目大意：给定一个字符串$S$，要求分成$A+prefix+B+middle+C+suffix$6段，满足：
$|prefix|=|suffix|$
$|middle|$为奇数
$prefix+middle+suffix$为回文串
除$middle$外所有段长度都可以为$0$
要求最大化$|prefix|+|middle|+|suffix|$，输出一组方案($|prefix|=|suffix|=0$时只输出$middle$)

首先我们发现$suffix$串是顶着右端点的，因此我们可以枚举$|suffix|$
对于每个$|suffix|$我们需要求出最左侧能与$suffix$匹配的$prefix$
我们发现随着$|suffix|$的增大，$prefix$的位置是单调不减的，因此我们可以用KMP算法求出对于每个$|suffix|$最靠左的$prefix$串
现在对于一个确定的$prefix$串和$suffix$串，我们需要求出中间的区间内最长的回文子串$middle$

现在就是区间最长回文子串的问题了

一种做法是这样的：
首先利用Manacher算法求出以每个点为中心的最长回文半径，然后对于每次询问我们二分答案：
设询问区间为$[L,R]$，二分长度为$M$，那么如果$[L+M-1,R-M+1]$区间最大的最长回文半径$\geq M$则存在回文半径为$M$的回文子串，否则不存在
利用ST表即可做到单次询问$O(logn)$

我还YY出了一种SB做法：
观察这个题的所有询问，我们发现询问区间是依次包含的
因此我们可以从最小的区间开始依次拓展
每次拓展，假设当前拓展的是$R$，那么我们需要用以$R$结尾，长度不超过$R-L+1$的最长回文子串来更新答案
这个用回文自动机+树上倍增就可以了

= =所以说是傻逼做法嘛

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define M 100100
#define MOD 999911657
#define BASE 2333
using namespace std;
int n,ans;
char s[M],rev_s[M];
int next[M],f[M],pos[M];
pair<int,int> intervals[M];
int log_2[M],a[M][17],_ans[M];
void KMP()
{
    int i,fix=0;
    for(i=2;i<=n;i++)
    {
        while( fix && rev_s[fix+1]!=rev_s[i] )
            fix=next[fix];
        if( rev_s[fix+1]==rev_s[i] )
            ++fix;
        next[i]=fix;
    }
    fix=0;
    for(i=1;i<=n;i++)
    {
        while( fix && rev_s[fix+1]!=s[i] )
            fix=next[fix];
        if( rev_s[fix+1]==s[i] )
            ++fix;
        if( i+fix>=n )
            return ;
        if(!pos[fix])
        {
            pos[fix]=i;
            intervals[fix]=make_pair(i+1,n-fix);
        }
    }
}
void Manacher()
{
    int i,id=1,mx=1;
    s[0]='$';
    for(i=1;i<=n;i++)
    {
        f[i]=min(mx-i+1,f[id+id-i]);
        while(s[i+f[i]]==s[i-f[i]])
            ++f[i];
        if(i+f[i]-1>mx)
            mx=i+f[i]-1,id=i;
    }
}
int Get_Max(int x,int y)
{
    int len=log_2[y-x+1];
    return max(a[x][len],a[y-(1<<len)+1][len]);
}
int Bisection(int x,int y)
{
    int l=1,r=y-x+2>>1;
    while(l+1<r)
    {
        int mid=l+r>>1;
        if( Get_Max(x+mid-1,y-mid+1)>=mid )
            l=mid;
        else
            r=mid;
    }
    return Get_Max(x+r-1,y-r+1)>=r?r:l;
}
int main()
{
    int i,j;
    scanf("%s",s+1);n=strlen(s+1);
    for(i=1;i<=n;i++)
        rev_s[i]=s[n-i+1];
    KMP();
    Manacher();
    pos[0]=1;intervals[0]=make_pair(1,n);
    for(i=2;i<=n;i++)
        log_2[i]=log_2[i>>1]+1;
    for(i=1;i<=n;i++)
        a[i][0]=f[i];
    for(j=1;j<=log_2[n];j++)
        for(i=1;i+(1<<j)-1<=n;i++)
            a[i][j]=max(a[i][j-1],a[i+(1<<j-1)][j-1]);
    for(i=0;pos[i];i++)
    {
        _ans[i]=2*Bisection(intervals[i].first,intervals[i].second)-1;
        if(i*2+_ans[i]>ans*2+_ans[ans])
            ans=i;
    }
    if(ans==0)
    {
        cout<<1<<endl;
        for(i=1;i<=n;i++)
            if(f[i]*2-1==_ans[0])
            {
                cout<<i-f[i]+1<<' '<<_ans[0]<<endl;
                return 0;
            }
    }
    else
    {
        cout<<3<<endl;
        cout<<pos[ans]-ans+1<<' '<<ans<<endl;
        int l=intervals[ans].first,r=intervals[ans].second;
        for(i=l;i<=r;i++)
            if( min(min(i-l,r-i)+1,f[i])*2-1==_ans[ans] )
            {
                cout<<i-(_ans[ans]>>1)<<' '<<_ans[ans]<<endl;
                break;
            }
        cout<<n-ans+1<<' '<<ans<<endl;
    }
    return 0;
}