字符串总结

回文串
Codeforces 835D
心累，博客又被吞了，本来不想写的，但这道题包含的知识点太多了，还是再写一遍吧。
该题题意:定义k阶回文串的左边右边为k-1阶回文串，而1阶回文串就是我们之前理解的回文串。

思路:该题的做法必然是先找到回文串，再判断该回文串是第几阶。
那么如何找一个字符串的所有回文串呢？
博主提供了3种方法来解决这个问题。

• 区间dp
  估计这是最常见的思路了。。代码也非常好理解。直接上代码啦。

#include <iostream>
#include <algorithm>
using namespace std;
int dp[5005][5005],ans[5005];
int main(){
    char s[5005];
    int i,j,l,len;
    while(scanf("%s",s+1)!=EOF)
    {
        len=strlen(s+1);
        memset(dp,0,sizeof(dp));
        memset(ans,0,sizeof(ans));
        for(l=1;l<=len;l++)
        {
            for(i=1;i<=(len-l+1);i++)
            {
                j=i+l-1;
                if(s[i]!=s[j])continue;
                if(l==1)
                    dp[i][j]=1;
                else if(l==2)
                    dp[i][j]=2;
                else if(dp[i+1][j-1]>0)
                    dp[i][j]=dp[i][i+l/2-1]+1;
                else
                    continue;
                ans[dp[i][j]]++;
            }
        }
        for(i=len,ans[len+1]=0;i>=1;i--)
            ans[i]+=ans[i+1];
        for(i=1;i<=len;i++)
        {
            if(i>1)
                printf(" ");
            printf("%d",ans[i]);
        }
        printf("\n");
    }
    return 0;
}

• 马拉车算法
  不了解马拉车算法的同学请点这里
  马拉车算法的最初用途是寻找字符串的最长回文串，那么如何用它来找所有的回文串呢，其实只需改一点就可以为我们所用了。

#include<iostream>
#include<stdio.h>
#include<string>
#include<queue>
#include<cstring>
using namespace std;
string str;
const int maxn=5005;
char tmp[maxn<<1];
int Len[maxn<<1];
queue<int>Q;
void init()
{
    int i,len=str.size();
    tmp[0]='@';
    for(i=1;i<=2*len;i+=2)
    {
        tmp[i]='#';
        tmp[i+1]=str[i/2];
    }
    tmp[2*len+1]='#';
    tmp[2*len+2]='$';
    tmp[2*len+3]=0;
}

void MANACHER()
{
    int len=strlen(tmp);
    int mx=0,ans=0,po=0;
    for(int i=1;i<=len;i++)
    {
        if(mx>i)Len[i]=min(mx-i,Len[2*po-i]);
        else Len[i]=1;
        while(tmp[i-Len[i]]==tmp[i+Len[i]])
            Len[i]++;
        if(Len[i]+i-1>mx)
        {
            mx=Len[i]+i-1;
            po=i;
        }
        if(Len[i]==i)
            Q.push(i);
    }
}

int main()
{
    cin>>str;
    int len=str.size();
    int ans[maxn<<1];
    memset(ans,0,sizeof(ans));
    int sum[maxn<<1];
    for(int i=len;i>=1;i--)
    {
        //len=i;
        init();
        MANACHER();
        Q.pop();
        memset(sum,0,sizeof(sum));
        while(!Q.empty())
        {
            int now=Q.front();Q.pop();
            int l=(now-1)/2+1,r=now+(now)/2;
            sum[now]=sum[l]+1;
            ans[sum[now]]++;
        }
        str=str.substr(1,len);
    }
    for(int i=len-1;i>=1;i--)
        ans[i]+=ans[i+1];
    for(int i=1;i<=len;i++)
    {
        printf("%d",ans[i]);
        if(i!=len)printf(" ");
    }
    printf("\n");
    return 0;
}

• Hash判断回文
  将其左右串转换成一个数字，判断是否相等，如果相等就是回文串了，再通过dfs来判断字串是第几阶就行了。

#include<iostream>
#include<stdio.h>
#include<string>
using namespace std;
const unsigned long long hashsize=233333;
const int maxn=5005;
const int MAXN=1e6+5;
int len;
string a;
unsigned long long p[maxn+10];
void init()
{
    p[0]=1;
    for(int i=1;i<=maxn;i++)
        p[i]=p[i-1]*hashsize;
}
unsigned long long S[maxn+10];
unsigned long long RS[maxn+10];
bool vis[maxn+10][maxn+10];
int dp[maxn+10][maxn+10];

unsigned long long HS(int l,int r)
{
    if(l==r)return a[l];
    l++,r++;
    return S[r]-S[l-1]*p[r-l+1];
}

unsigned long long RH(int l,int r)
{
    if(l==r)return a[l];
    l=len-l-1,r=len-r-1;swap(l,r);
    l++;r++;
    return RS[r]-RS[l-1]*p[r-l+1];

}

void DFS(int l,int r)
{
    if(vis[l][r])return ;
    vis[l][r]=1;
    if(l==r){dp[l][r]=1;return ;}
    int rr=(l+r+1)/2-1,ll=(l+r)/2+1;
    if(HS(ll,r)==RH(l,rr))dp[l][r]=1;
    DFS(l,rr);DFS(ll,r);
    if(dp[l][r]&&HS(ll,r)==HS(l,rr))
        dp[l][r]=dp[ll][r]+1;
}
int ans[maxn];
int main()
{
    init();

    cin>>a;
    len=a.size();
    for(int i=1;i<=len;i++)
        S[i]=S[i-1]*hashsize+a[i-1];
    for(int i=1;i<=len;i++)
        RS[i]=RS[i-1]*hashsize+a[len-i];
    for(int i=0;i<len;i++)
        for(int j=i;j<len;j++)
        DFS(i,j);
    for(int i=0;i<len;i++)
        for(int j=i;j<len;j++)
        ans[dp[i][j]]++;
    for(int i=len-1;i>=1;i--)
        ans[i]+=ans[i+1];
    for(int i=1;i<=len;i++)
    {
        printf("%d",ans[i]);
        if(i!=len)printf(" ");
    }
    return 0;
}