最小表示法:
int getMin(int *str,int len){
int p1=0,p2=1,k=0;
while(p1<len && p2<len && k<len){
int dif=str[p1+k]-str[p2+k];
if(dif){
if(dif>0) p1= (p1+k+1<=p2)? p2:p1+k+1;
else p2= (p2+k+1<=p1)? p1:p2+k+1;
k=0;
if(p1==p2) p2++;
}else k++;
}
return min(p1,p2);
}
KMP:
#define type char
int jump[M];
type str[M],pat[M];
void getNext(type *p,int *jump,int pLen){
jump[0]=-1;
for(int i=1,j=-1;i<pLen;i++){
while(j!=-1&&p[j+1]!=p[i]) j=jump[j];
jump[i]= p[j+1]==p[i]? ++j:j;
}
}
int kmp(type *s,type *p,int sLen,int pLen){
getNext(p,jump,pLen);
int ret=-1;
for(int i=0,j=-1;i<sLen;i++){
while(j!=-1&&p[j+1]!=s[i]) j=jump[j];
if(p[j+1]==s[i]) j++;
if(j+1==pLen){
ret=i-pLen;break;//求最早匹配位置
j=jump[j],ret++;//求匹配数
}
}
return ret+1;
}
struct manacher{
char init[M],now[M<<1];
int p[M<<1];//p[i]-1表示当前位置为中心的最大回文串长度
int read(){return scanf("%s",init);}
void rebuild(){
int i,j;
now[0]='*',now[1]='#';//添加的字符必须是串中没有出现过的字符
for(i=0,j=2;init[i];i++,j+=2)
now[j]=init[i],now[j+1]='#';
now[j]=0;
}
int cal(){
rebuild();
int i,mx=0,id,res=1;
for(i=0;now[i];i++){
p[i]=(mx>i? _min(p[2*id-i],mx-i):1);
for(;now[i-p[i]]==now[i+p[i]];p[i]++);
if(mx<i+p[i]) mx=i+p[i],id=i;
res=_max(res,p[i]);
}
return res-1;
}
}p;