hdu5371 manacher + 线段树

这题说的找出一个数组串 3等分 第一个部分和第3个部分一样，第二个部分和第一个部分回文，那么计算出这些字符串问这样的字符串最长为多少，我们先使用manacher 计算出每个位置以他为对称轴左边端点的最长回文半径 加入第i个位置 回文半径为 d[i]，那么他能影响的范围为d[i]-i至i 如果他是那个最大的那么答案就是 i-j （j为[d[i]-i,i]之间的数），我们可以再适当的点加入这个位置的影响在线段树中，然后在最后使用线段树查找出这个区间的最值是多少

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <string.h>
#include <vector>
using namespace std;
const int maxn=;
int str[maxn];
int s[maxn*];
int rad[maxn*];//以i为中心的回文串 半径
int H[maxn],to[maxn],nx[maxn],numofedg;
int d[maxn];
void addedg(int u, int v)
{
     ++numofedg; to[numofedg]=v; nx[numofedg]=H[u]; H[u]=numofedg;
}
void manacher(int len)
{

     int cur=;
     s[]=-;
     for(int i=; i<len; i++)
     {
         s[cur++]=str[i];
         s[cur++]=-;
     }
     rad[]=;
     int id=;
     for(int i=; i<cur; i++)
        {
            if(rad[id]+id>i)rad[i]=min(rad[*id-i],rad[id]+id-i);
            else rad[i]=;
            while(i-rad[i]>=&&i+rad[i]<cur&&s[ i - rad[i] ] == s[ i+rad[i] ])rad[i]++;
            if(id+rad[i]<i+rad[i])id=i;
        }
     for(int i=; i<=len; i++)
        {
               d[i]=rad[i*]/;
               if(d[i]==)continue;
               addedg(i-d[i],i);
        }

}
struct Itree
{
      int loc,cL,cR,ans;
      int val[maxn*];
      void build(int L, int R, int o)
      {
          val[o]=;
          if(L==R) return ;
          int mid=(L+R)>>;
          build(L,mid,o*);
          build(mid+,R,o*+);
      }
      void update(int L, int R, int o)
      {
           if(L==R)
            {
                 val[o]=L ; return ;
            }
           int mid=(L+R)>>;
           if(loc<=mid) update(L,mid,o*);
           else update(mid+,R,o*+);
           val[o]=max(val[o*],val[o*+]);
      }
     void query(int L, int R, int o)
     {
         if(cL<=L&&R<=cR)
         {
             ans=max(ans,val[o]); return;
         }
         int mid=(L+R)>>;
         if(cL<=mid)query(L,mid,o*);
         if(cR>mid) query(mid+,R,o*+);
     }
}T;
int main()
{
   int cas;
   scanf("%d",&cas);

   for(int cc=; cc<=cas; cc++)
   {
       int n;
       scanf("%d",&n);
       H[]=;numofedg=;
       for(int i=; i<n; i++){
            H[i+]=;
            scanf("%d",&str[i]);
       }
       manacher(n);
       int ans=;
       T.build(,n,);
       d[]=;
       for(int i=; i<n; i++)
        {
             for(int j=H[i]; j; j=nx[j])
             {
                 T.loc=to[j]; T.update(,n,);
             }
             T.cL=i+;T.cR=d[i]+i;
             if(T.cL>T.cR)continue;
             T.ans=;
             T.query(,n,);
             ans=max(ans,T.ans-i);
        }
        printf("Case #%d: %d\n",cc,ans*);
   }
    return ;
}