ural1297最长回文串

WA在test16过不了，我并觉得我的程序没毛问题- -

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define maxl 201000

using namespace std;

int insn,n,ansind,ansl;
int c[maxl],rnk[maxl],sa[maxl],height[maxl],X[maxl],Y[maxl];
int st[maxl][20]; 
char ch[maxl];

void prework()
{
	memset(st,0,sizeof(st));
	memset(sa,0,sizeof(sa));
	memset(rnk,0,sizeof(rnk));
	memset(sa,0,sizeof(sa)); 
	memset(height,0,sizeof(height));
	n=strlen(ch);insn=n;
	ch[insn]=1;
	int j=n-1;
	for(int i=insn+1;i<=2*insn;i++)
		ch[i]=ch[j],j--;
	n=2*insn+1;ch[2*insn+1]=0;
}

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

void da(int m,int n)
{
	int p=0,*x=X,*y=Y;
	for(int i=0;i<m;i++) c[i]=0;
	for(int i=0;i<n;i++) c[x[i]=ch[i]]++;
	for(int i=1;i<m;i++) c[i]+=c[i-1];
	for(int i=n-1;i>=0;i--) sa[--c[x[i]]]=i;
	for(int j=1;p<n;j<<=1)
	{
		p=0;
		for(int i=n-j;i<n;i++) y[p++]=i;
		for(int i=0;i<n;i++) if(sa[i]>=j) y[p++]=sa[i]-j;
		for(int i=0;i<m;i++) c[i]=0;
		for(int i=0;i<n;i++) c[x[y[i]]]++;
		for(int i=1;i<m;i++) c[i]+=c[i-1];
		for(int i=n-1;i>=0;i--) sa[--c[x[y[i]]]]=y[i];
		p=1;swap(x,y);x[sa[0]]=0;
		for(int i=1;i<n;i++) x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;
		m=p;
	}
}

void getheight(int n)
{
	int k=0;
	for(int i=1;i<=n;i++)
		rnk[sa[i]]=i;
	for(int i=0;i<n;i++)
	{
		if(k) k--;
		int j=sa[rnk[i]-1];
		while(ch[i+k]==ch[j+k])
			k++;
		height[rnk[i]]=k;
	}
}

void RMQ()//RMQ???
{
    int i,j;
    memset(st,127,sizeof(st));
    for(i=1;i<=insn*2+1;i++)
        st[i][0]=height[i];
    for(j=1;(1<<j)<=2*insn+1;j++)
        for(i=1;i+(1<<j)-1<=2*insn+1;i++)
            st[i][j]=min(st[i][j-1],st[i+(1<<(j-1))][j-1]);
}

int lcp(int l,int r)//???????
{
    int a=rnk[l],b=rnk[r];
    if(a>b)
        swap(a,b);
    a++;
    int t=(int)(log(double(b-a+1))/log(2.00));
    return min(st[a][t],st[b-(1<<t)+1][t]);
}

void mainwork()
{
	da(256,n+1);
	getheight(n);
	ansl=0;ansind=maxl-1;
	int l=log2(n);
	RMQ();
	int s,t,len,d;
/*	for(int i=0;i<insn;i++)
	{
		s=rnk[i];t=rnk[2*insn-i];
		if(s>t) swap(s,t);
		s++; 
		len=t-s+1;l=(int)(log(double(t-s+1))/log(2.00));
		d=min(st[l][s],st[l][t-(1<<l)+1]);
		if(2*d-1>ansl)
			ansl=2*d-1,ansind=i-(d-1);
		else
		if(2*d-1==ansl)
			ansind=min(ansind,i-(d-1));

		if(i>0)
		{ 
			s=rnk[i];t=rnk[2*insn-i+1];
			if(s>t) swap(s,t);
			s++;
			len=t-s+1;l=(int)(log(double(t-s+1))/log(2.00));
			d=min(st[l][s],st[l][t-(1<<l)+1]);	
			if(2*d>ansl)
				ansl=2*d,ansind=i-1-(d-1);
			else
			if(2*d==ansl)
				ansind=min(ansind,i-1-(d-1));
		}
	}*/
	int res,flag;ansl=0;
	for(int i=0;i<insn;i++)
    {
        res=lcp(i,2*insn-i)*2-1;
        if(ansl<res)//???
        {
        	ansl=res;
            flag=i;
        }
        if(i>0)//???
        {
        	res=lcp(i,2*insn-i+1)*2;
            if(ansl<res)
            {
                ansl=res;
                flag=i;
            }
        }
    }
        if(ansl%2==1)
            for(int i=flag-ansl/2;i<=flag+ansl/2;i++)
                printf("%c",ch[i]);
        else
            for(int i=flag-ansl/2;i<=flag+ansl/2-1;i++)
                printf("%c",ch[i]);
        printf("\n");
}

void print()
{
/*	for(int i=1;i<=ansl;i++)
		printf("%c",ch[ansind+i-1]);
	printf("\n");*/
}

int main()
{
	while(~scanf("%s",ch))
	{
		prework();
		mainwork();
		print();
	}
	return 0;
}