回文串专题

由于博主比较水嘛，所总结的可能不全面，如果各位有什么好的思路可以在下方留言交流！欢迎各位大神前来踩。
什么是回文串？
回文串指的是正读和反读都一样的字符串，如aba，abba等。
那么怎么来判断回文串呢？
这里的思路是同时从前后两端向中间扫描如果有一对字符不能匹配到，就说明不是回文串，如果把所有的扫描完或者只剩下一个未扫描时就说明它是回文串。
代码：

int pend(char s[],int n)
{
    int i,j;
    i=0;j=n-1;
    while(i<j)
    {
        if(s[i]!=s[j])
            return 0;
        i++;
        j--;
    }
    return 1;
}

这里有一个回文串的题目，上链接：

hdu2163
题目不难，判断是否为字符串暴力即可；

#include <stdio.h>
#include <string.h>
int pendd(char s[],int n);
int main()
{
    char s[60];
    int n,t=0;
    while(~scanf("%s",s))
    {
        t++;
        n=strlen(s);
        if(n==4)
            if(strcmp(s,"STOP")==0)
            break;
        if(pendd(s,n)==1)
        printf("#%d: YES\n",t);
        else
        printf("#%d: NO\n",t);
    }
    return 0;
}
int pendd(char s[],int n)
{
    int i,j;
    i=0;j=n-1;
    while(i<j)
    {
        if(s[i]!=s[j])
            return 0;
        i++;
        j--;
    }
    return 1;
}

再举一题
问题 : 小L的逆袭

时间限制: 1 Sec 内存限制: 128 MB
提交: 5 解决: 2
[提交][状态]
题目描述

好啦好啦、不是逆袭女神神马的，想到这里的童鞋都给我面壁去。。

集训队里面的大M和大W、对小L是亦师亦队友的关系，他们之间偶尔也会弄一些小比赛来乐呵乐呵，当然小L的获胜战绩至今为 0。。 T^T + 囧~。今天他们又手痒难耐了、那比什么类？对了、据说字符串里面有种字符串叫回文字符串呀，那我们不如就来统计这个吧。


最后Lh成功获得了仅有的一次胜利、并获得两位土豪实现承诺的彩头：两位土豪正式工作后、头两年的工资的50%、无偿归Lh所有。

Lh：“(￣﹁￣)”想入非非ing。。。

两位土豪：“哪有什么彩头？又做梦了吧你。。。”

正题：

好了，擦掉口水、回归正题：

给一个n*m的矩阵，统计里面有多少回文串，我们只横向统计和纵向统计。比如对于2*3的矩阵：

bcb

aca

有9个回文串：

    六个字母各自算一个；

    纵向有一个回文串cc；

    横向有两个回文串，分别是bcb、aca

聪明的你知道答案是多少么？

回文串的定义：

一个字符串，如果从前往后看、和从后往前看、一模一样，那么这个字符串就是回文串。比如aba，就是一个回文串；而abaa就不是。

输入

输入数据有T组，第一行一个整数T（0<T<=10）。

每组数据第一行一对整数n、m（0<n,m<=50），接下来n行、每行一个长为m的字符串（只由小写字母构成）。

输出

对于每组数据，输出一个整数，表示有多少回文串。

样例输入

1
2 3
bcb
aca
样例输出

9
题意：给出一个字符串矩阵，让求出其中共有多少个回文串，由于数据不大，我暴力水过了。

#include <stdio.h>
#include <string.h>
int query(char s[],int n);
int pend(char s[],int n);
int main()
{
    char str[60][60],s[60];
    int sum=0;
    int m,n,i,j,t;
    scanf("%d",&t);
    while(t--)
    {
        sum=0;
        scanf("%d%d",&m,&n);
        for(i=0;i<m;i++)
        {
            scanf("%s",str[i]);
        }
        for(i=0;i<m;i++)
        {
            for(j=0;j<n;j++)
            {
                s[j]=str[i][j];
            }
            sum+=query(s,n);
        }
        for(i=0;i<n;i++)
        {
            for(j=0;j<m;j++)
            {
                s[j]=str[j][i];
            }
            s[j]='\0';
            sum+=query(s,m);
        }
        printf("%d\n",sum+m*n);
    }
    return 0;
}
int query(char s[],int n)
{
    int i,j,k,v=0;
    char m[60];
    for(i=0;i<n-1;i++)
    {
        k=0;
        for(j=i;j<n;j++)
        {
            m[k++]=s[j];
            m[k]='\0';
            if(j==i)
                continue;
            if(pend(m,j-i+1))
                v++;
        }
    }
    return v;
}
int pend(char s[],int n)
{
    int i,j;
    i=0;j=n-1;
    while(i<j)
    {
        if(s[i]!=s[j])
            return 0;
        i++;
        j--;
    }
    return 1;
}

接下来这个题就没那么水了，上链接。
hdu3948
题意是给你一个串，让你求出其中有多少个不重复的回文串。
思路是后缀数组什么的没看懂，这两天再看看，先把代码沾出来

#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;
#define MAXN 200100
int dp[MAXN][20];
char r[MAXN];
int sa[MAXN];
int wa[MAXN],wb[MAXN],wv[MAXN],ws[MAXN];
int height[MAXN],rk[MAXN];
int mm[MAXN];
int cas;
inline bool cmp(int *r,int a,int b,int len){
    return r[a]==r[b]&&r[a+len]==r[b+len];
}
void SA(int n,int m){
    int i,j,p,*x=wa,*y=wb,*t;
    for(i=0;i<m;i++)
        ws[i]=0;
    for(i=0;i<n;i++)
        ws[x[i]=r[i]]++;
    for(i=1;i<m;i++)
        ws[i]+=ws[i-1];
    for(i=n-1;i>=0;i--)
        sa[--ws[x[i]]]=i;
    for(j=p=1;p<n;j<<=1,m=p){
        for(p=0,i=n-j;i<n;i++)
            y[p++]=i;
        for(i=0;i<n;i++){
            if(sa[i]>=j)
                y[p++]=sa[i]-j;
        }
        for(i=0;i<m;i++)
            ws[i]=0;
        for(i=0;i<n;i++)
            ws[wv[i]=x[y[i]]]++;
        for(i=1;i<m;i++)
            ws[i]+=ws[i-1];
        for(i=n-1;i>=0;i--)
            sa[--ws[wv[i]]]=y[i];
        for(t=x,x=y,y=t,x[sa[0]]=0,p=i=1;i<n;i++)
            x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;
    }
}
void Height(int n){
    int i,j,k=0;
    for(i=0;i<=n;i++)    //这里sa[0]为‘\0’开始的子串
        rk[sa[i]]=i;
    for(i=0;i<n;height[rk[i++]]=k)
        for(k?k--:0,j=sa[rk[i]-1];r[i+k]==r[j+k];k++);
}
int init(){
    int i,j,len;
    memset(height,0,sizeof(height));
    len=strlen(r);
    r[len]='$';
    for(j=0;j<len;j++)
        r[j+len+1]=r[len-1-j];
    r[j+len+1]='\0';
    return len*2+1;
}
void st(int n){
    int i,j,p,q;
    for(i=1;i<=n;i++)
        dp[i][0]=i;
    for(j=1;j<=mm[n];j++)
        for(i=1;i+(1<<j)-1<=n;i++){
            p=height[dp[i][j-1]];
            q=height[dp[i+(1<<(j-1))][j-1]];
            if(p>q)
                dp[i][j]=dp[i+(1<<(j-1))][j-1];
            else
                dp[i][j]=dp[i][j-1];
        }

}
int RMQ_MIN(int i,int j){
    int tem;
    if(i>j){
        tem=i;
        i=j;
        j=tem;
    }
    i++;           //交换后小的要加一
    int k=mm[(j-i+1)];
    return min(height[dp[i][k]],height[dp[j-(1<<k)+1][k]]);
}
void solve(int n){
    int i,j,ans=0,pre1=0,pre2=0;
    st(n);
    for(i=2;i<=n;i++){
        pre1=min(pre1,height[i]);
        j=RMQ_MIN(i,rk[n-1-sa[i]]);
        if(j>pre1){
            ans+=j-pre1;
            pre1=j;
        }
        pre2=min(pre2,height[i]);
        j=RMQ_MIN(i,rk[n-sa[i]]);
        if(j>pre2){
            ans+=j-pre2;
            pre2=j;
        }
    }
    printf("Case #%d: %d\n",cas,ans);
}
int main(){
    int i,j,n,t,T;
    mm[0]=-1;
    for(i=1;i<MAXN;i++)
        mm[i]=((i&(i-1))==0)/*i=2的次方时i&(i-1)=0*/?mm[i-1]+1:mm[i-1];           //这里预处理一下，时间跑快了

    scanf("%d",&T);
    cas=0;
    for(t=1;t<=T;t++){
        scanf("%s",r);
        n=init();
        SA(n+1,130);
        Height(n);
        cas++;
        solve(n);
    }
}

再来一个最长回文串（manacher）
hdu3068

#include<stdio.h>
#include<string.h>
#define MAX 111000
char s[MAX*2];
int p[MAX*2];
int min(int a,int b)
{
    return a>b?b:a;
}
int main(){
    while(scanf("%s",s)!=EOF){
        int len=strlen(s),id=0,maxlen=0;
        for(int i=len;i>=0;--i){//插入'#'
            s[i+i+2]=s[i];
            s[i+i+1]='#';
        }//插入了len+1个'#',最终的s长度是1~len+len+1即2*len+1,首尾s[0]和s[2*len+2]要插入不同的字符
        s[0]='*';//s[0]='*',s[len+len+2]='\0',防止在while时p[i]越界
        for(int i=2;i<2*len+1;++i){
            if(p[id]+id>i)p[i]=min(p[2*id-i],p[id]+id-i);
            else p[i]=1;
            while(s[i-p[i]] == s[i+p[i]])++p[i];
            if(id+p[id]<i+p[i])id=i;
            if(maxlen<p[i])maxlen=p[i];
        }
        printf("%d\n",maxlen-1);
    }
    return 0;
}

后两个都用到了manachar算法，之后我会更新一篇关于它的总结，请持续关注~~~~。