hdu 3068 最长回文子串 TLE

后缀数组+RMQ是O(nlogn)的，会TLE.....

标准解法好像是马拉车，O(n)....

 #include "algorithm"
 #include "cstdio"
 #include "cstring"
 using namespace std;
 #define maxn 220020

 int wa[maxn],wb[maxn],wv[maxn],ws[maxn];
 int rank[maxn],height[maxn];
 int r[maxn],sa[maxn],RMQ[maxn][];
 char s[];
 int n;

 int cmp(int *r,int a,int b,int l)
 {
     return r[a]==r[b]&&r[a+l]==r[b+l];
 }

 void da(int *r,int *sa,int n,int m)
 {
     int i,j,p,*x=wa,*y=wb,*t;
     for(i=; i<m; i++) ws[i]=;     //注意此处有个小问题：ws和iostream里面某地重名了..
     for(i=; i<n; i++) ws[x[i]=r[i]]++;
     for(i=; i<m; i++) ws[i]+=ws[i-];
     for(i=n-; i>=; i--) sa[--ws[x[i]]]=i;
     for(j=,p=; p<n; j*=,m=p)
     {

         for(p=,i=n-j; i<n; i++) y[p++]=i;
         for(i=; i<n; i++) if(sa[i]>=j) y[p++]=sa[i]-j;
         for(i=; i<n; i++) wv[i]=x[y[i]];
         for(i=; i<m; i++) ws[i]=;
         for(i=; i<n; i++) ws[wv[i]]++;
         for(i=; i<m; i++) ws[i]+=ws[i-];
         for(i=n-; i>=; i--) sa[--ws[wv[i]]]=y[i];
         for(t=x,x=y,y=t,p=,x[sa[]]=,i=; i<n; i++)
             x[sa[i]]=cmp(y,sa[i-],sa[i],j)?p-:p++;
     }
     return;
 }

 void calheight(int *r,int *sa,int n)
 {
     int i,j,k=;
     for(i=; i<=n; i++) rank[sa[i]]=i;
     for(i=; i<n; height[rank[i++]]=k)
         for(k?k--:,j=sa[rank[i]-]; r[i+k]==r[j+k]; k++);
     return;
 }

 bool _same(int lx,int ly,int l1,int l2)
 {
     return (((lx<=l1-) && (ly>=l1+))||((ly<=l1-) && (lx>=l1+)));
 }

 void ST()        //初始化
 {
     memset(RMQ,,sizeof(RMQ));
     for(int i=;i<=n;i++)
         RMQ[i][]=height[i];
     for(int j=;(<<j)<=n;j++)
         for(int i=;i+(<<j)-<=n;i++)
         RMQ[i][j]=min(RMQ[i][j-],RMQ[i+(<<(j-))][j-]);
 }

 int Query(int L,int R)    //求a[L..R]区间的最值
 {
     int k=;
     while((<<(k+))<=R-L+) k++;
     int tb=min(RMQ[L][k],RMQ[R-(<<k)+][k]);
     return tb;
 }

 int calc(int l,int r)       //求l开始后缀和r开始后缀的最长公共前缀
 {
     int tl=rank[l],tr=rank[r];
     if (tl>tr)      swap(tl,tr);
     tl++;
     int ans=Query(tl,tr);       //相当于RMQ问题
 //  printf("calc:  %d %d -- %d %d == %d\n",l,r,tl,tr,ans);
     return ans;
 }

 int _max(int a,int b,int c,int d)
 {
     int mx=a;
     if (b>mx) mx=b;
     if (c>mx) mx=c;
     if (d>mx) mx=d;
     return mx;
 }

 //da(r,sa,n+1,128);
 //calheight(r,sa,n);
 int main()
 {
     while(~scanf("%s",s))
     {
         n=strlen(s);

         for (int i=; i<n; i++)
             r[i]=int(s[i])-int('a')+;

         r[n]=;
         for (int i=n-; i>=; i--)
             r[*n-i]=int(s[i])-int('a')+;
         r[*n+]=;
         int n2=n;       n=*n+;

 //      for (int i=0; i<=n; i++)      printf("%d ",r[i]);
 //        printf("\n %d   %d\n",n,n2);

         da(r,sa,n+,);
         calheight(r,sa,n);

 //      for (int i=0; i<=n; i++)
 //          printf("%d  %d  %d\n",sa[i],height[i],rank[i]);
 //      printf("\n");

         ST();
         int ans=;
         for (int i=; i<=n2-; i++)
         {
             int t1=calc(i,n-i)*;
             int t2=calc(i+,n-i)*+;
             int t3=calc(i+,n-i-)*;
 //          printf("%d %d %d\n",t1,t2,t3);
             ans=_max(ans,t1,t2,t3);
         }
         printf("%d\n",ans);
     }
     return ;
 }

最近两天有点不在状态....先滚去整理模板吧