UVA 11475 Extend to Palindrome KMP，后缀数组

     给一个字符串，尽可能少的添加字符使他变成一个回文串，求添加之后的回文串。

    这题kmp,后缀数组都能做，大体思路就是找一个最长的回文后缀，那么要添加的字符就是字符串除掉这个后缀后剩下的前缀了。找的时候可以反向拼接，构造后缀数组然后二分长度+遍历height判断；或者直接把原串接在翻转串的后面，求一下失配函数。不管是效率还是实现，感觉都是KMP好一点....

KMP的代码:

#include <iostream>
#include <cstring>
#include <algorithm>
#include <cstdio>
#include <cmath>
using namespace std;
typedef long long ll;
const int maxn=410000;
bool ok[maxn];
int f[maxn];
char s[maxn],ss[maxn];
int cnt[maxn];
void getFail(char* P,int* f)
{
    int m=strlen(P);
    f[0]=0;
    f[1]=0;
    for (int i=1; i<m; i++)
    {
        int j=f[i];
        while(j && P[i]!=P[j]) j=f[j];
        f[i+1]=P[i]==P[j]?j+1:0;
    }
}
bool find(char* T,char* P,int* f)
{
    int n=strlen(T);
    int m=strlen(P);
//    getFail(P,f);
    int j=0;
    for (int i=0; i<n; i++)
    {
        while(j && P[j]!=T[i]) j=f[j];
        if (P[j]==T[i]) j++;
        if (j==m) return true;
    }
    return false;
}
int tt;
int n,m,k;

int main()
{
//    freopen("in.txt","r",stdin);
    while(~scanf("%s",ss))
    {
        int len=strlen(ss);
        for (int i=len-1; i>=0; i--)
        s[len-i-1]=ss[i];
        s[len]='\0';
        strcat(s,ss);
        getFail(s,f);
        int ans=min(len,f[strlen(s)]);
        printf("%s",ss);
        int ct=len-ans-1;
        for (ct;ct>=0; ct--)
        printf("%c",ss[ct]);
        puts("");
    }
    return 0;
}

后缀数组的代码：

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <string>
#include <cstring>
using namespace std;
typedef long long ll;
const int maxn=410100;
int s[maxn];
int sa[maxn],t[maxn],t2[maxn],c[maxn];
int sa2[maxn];
char s1[maxn],s2[maxn];
int rank[maxn],height[maxn];
int l1,l2;
int len,len1;
int n,m,ans;
int belong[maxn];
inline int idx(char s)
{
   if (s>='a' && s<='z') return s-'a'+2;
   if (s>='A' && s<='Z') return s-'A'+2+26;
//     return s-min('a','A')+101;

}
void getheight(int n)
{
    int i,j,k=0;
    for (i=0; i<=n; i++) rank[sa[i]]=i;

    for (i=0; i<n; i++)
    {
        if (k) k--;
        int j=sa[rank[i]-1];
        while(s[i+k]==s[j+k]) k++;
        height[rank[i]]=k;
    }
}

void build_ss(int m,int n)
{
    n++;
    int i,*x=t,*y=t2;
    for (int i=0; i<m; i++) c[i]=0;
    for (int i=0; i<n; i++) c[x[i]=s[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 k=1; k<=n; k<<=1)
    {
        int p=0;
        for (i=n-k; i<n; i++) y[p++]=i;
        for (i=0; i<n; i++) if (sa[i]>=k) y[p++]=sa[i]-k;

        for (i=0; i<m; i++) c[i]=0;
        for (i=0; i<n; i++) c[x[y[i]]]++;
        for (i=1; i<m; i++) c[i]+=c[i-1];
        for (i=n-1; i>=0; i--) sa[--c[x[y[i]]]] = y[i];
        swap(x,y);
        p=1;
        x[sa[0]]=0;
        for (i=1; i<n; i++)
        x[sa[i]]=(y[sa[i-1]]==y[sa[i]] && y[sa[i-1]+k]==y[sa[i]+k])? p-1 : p++;
        if (p>=n) break;
        m=p;
    }
}
int d[maxn][22];

void RMQ_init()
{
    for (int i=1; i<=n; i++) d[i][0]=height[i];
     for (int j=1; (1<<j)<=n; j++)
      for (int i=1; (i+(1<<j)-1)<=n; i++)
      d[i][j]=min(d[i][j-1],d[i+(1<<(j-1))][j-1]);
}
int RMQ(int L,int R)
{
    if (L>R) swap(L,R);
    L++;

    int k=0;
    while((1<<(k+1))<=R-L+1) k++;
    return min(d[L][k],d[R-(1<<k)+1][k]);
}
void init()
{
    l1=strlen(s1);
    n=0;
    for (int i=0; i<l1; i++)
    belong[n]=0,s[n++]=idx(s1[i]);
    belong[n]=-1;
    s[n++]=1;
    for (int i=l1-1; i>=0; i--)
    {
        belong[n]=1;
        s[n++]=idx(s1[i]);
    }
    belong[n]=-1;
    s[n++]=0;


    build_ss(60,n);
    getheight(n);
}

bool check(int k)
{
    int sta=0;
    int p=-1,q=-1;
    for (int i=1; i<=n; i++)
    if (height[i]>=k)
    {
        sta|=(1<<belong[sa[i]]);
        sta|=(1<<belong[sa[i-1]]);
        if (belong[sa[i]]==0 && sa[i]+height[i]==l1) p=sa[i];
        if (belong[sa[i-1]]==0 && sa[i-1]+height[i]==l1) p=sa[i-1];
        if (belong[sa[i]]==1 && sa[i]==l1+1) q=sa[i];
        if (belong[sa[i-1]]==1 && sa[i-1]==l1+1) q=sa[i-1];
        if (p>=0 && q>=0&&  sta==3) return true;
    }
    else
    {
        p=-1;
        q=-1;
        sta=0;

    }
    return false;
}
void slove()
{
    int l=1,r=l1+1;
    ans=0;
    while(l<r)
    {
        int mid=(l+r)>>1;
        if (check(mid))
        {
            ans=mid;
            l=mid+1;
        }
        else r=mid;
    }
    int ct=l1-ans-1;
    printf("%s",s1);
    for (ct;ct>=0; ct--)
    printf("%c",s1[ct]);
    puts("");
}
int tt;
int main()
{
//      freopen("in.txt","r",stdin);
      while(~scanf("%s",&s1))
      {
         init();
         slove();
      }
}