- 对于KMP算法,都会涉及到两个串
一个是待匹配串
char T[1000]//待匹配串
待匹配串的长度定为n
一个是模板串 char P[100]//模板串 模板串的长度定为m
- 要知道KMP的F[i]数组求得的数值就是串P中的[1,i-1]的后缀与串P中的[0,i-2]前缀的最大匹配长度。
对于KMP算法一般会利用到两个函数
void find(char *T, char *P, int *f) //找到所有匹配点
void getFail(char *P, int *f)
- KMP算法经常用到的两个函数的模板
const int MAXN=1000000+100;
const int MAXM=10000+100;
int T[MAXN],P[MAXM],f[MAXM],f2[MAXM];
int n,m;
void getFail(char *P,int *f)//必须用f和f2两个数组,要不然不好计算优化指针,比如串aaaab的优化指针是000030
{
f[0]=f[1]=0;//保存优化之后的失配指针
f2[0]=f2[1]=0;//f2保存未优化的适配指针
int m =strlen(P);
for(int i=1;i<m;i++)
{
int j=f2[i];
while(j && P[i]!=P[j]) j=f2[j];
f2[i+1] = f[i+1] = (P[i]==P[j])?j+1:0;//既然i+1的失配位置指向j+1,但是P[i+1]和P[j+1]的内容是相同的
//所以就算指针从i+1跳到j+1去,还是不能匹配,所以f[i+1]直接=f[j+1]
if(f[i+1]==j+1 && P[i+1]==P[j+1]) f[i+1]=f[j+1];
//cout<<f[i]<<f2[i]<<endl;
}
}
const int MAXN=1000000+100;
const int MAXM=10000+100;
int T[MAXN],P[MAXM],f[MAXM],f2[MAXM];
int n,m;
//int cnt;
//int ex[MAXN];
void find(char *T,char *P,int *f)
{
int n=strlen(T);
int m=strlen(P);
int j=0;
for(int i=0;i<n;i++)
{
if(j && T[i]!=P[j]) j=f[j];
if(T[i]==P[j]) j++;
//ex[i]=j; 在用KMP计算两字符串前缀后缀的时候常用到(注:从此处开始下面的代码一般都以具体题目而定)
if(j==m)//j为两串相同的长度,i为在T串中出现的位置
{
//cnt++;
//printf("%d\n", i - m + 1);//输出每次P串在T串中每一次出现的位置,注意下标是从0开始的
//return i-m+1; (前面要变成 int find函数)
}
}
}
- 还有一种求多个字符串的最长公共连续子串,如果存在多个长度相同的就输出字典序最小的那个,其模板为
char s[4100][250];
char P[250];
char ans[250];//最终的结果字符串
int m,next[250];
int n;//表示有多少个主串
void getFail()
{
m=strlen(P);
next[0]=next[1]=0;
for(int i=1;i<m;i++)
{
int j =next[i];
while(j && P[i]!=P[j]) j=next[j];
next[i+1] = (P[i]==P[j])?j+1:0;
}
}
bool find(char *T)
{
int j=0;
getFail();
int n=strlen(T);
for(int i=0;i<n;i++)
{
while(j && T[i]!=P[j]) j=next[j];
if(T[i]==P[j]) j++;
if(j==m) return true;
}
return false;
}
bool solve_1(int len)//判断第0个base串的长为len的字串中是否有可行串
{
int s0_len=strlen(s[0]);
for(int i=0;i+len-1<s0_len;i++)
{
strncpy(P,s[0]+i,len);
P[len]=0;//P串末尾加'\0'
bool ok=true;
for(int j=1;j<n;j++)
if(!find(s[j]))
{
ok=false;
break;
}
if(ok) return true;
}
return false;
}
void solve_3(int len)//找出长度为len的可行字串中字典序最小的,放在ans中
{
bool first=true;
int s0_len=strlen(s[0]);
for(int i=0;i+len-1<s0_len;i++)
{
strncpy(P,s[0]+i,len);
P[len]=0;//P串末尾加'\0'
bool ok=true;
for(int j=1;j<n;j++)if(!find(s[j]))
{
ok=false;
break;
}
if(ok)
{
if(first)
{
strncpy(ans,P,len+1);
first=false;
}
else if(strcmp(ans,P)>0)//字典序小才更新
strncpy(ans,P,len+1);
}
}
}
bool solve()
{
int L=1,R=strlen(s[0]);
if(!solve_1(1)) return false;
while(R>L)
{
int m=L+(R-L+1)/2;
if(solve_1(m))
L=m;
else
R=m-1;
}
//找到了可行字串的最长长度为L,然后需要找出字典序最小的
solve_3(L);//找到长L的字典序最小的字串存在ans中
return true;
}
KMP模板题
POJ 3461Oulipo(KMP:统计一个串出现的次数):KMP模板题。解题报告!
HDU 1711Number Sequence(KMP:找模板第一次出现的位置):KMP模板题。解题报告!
HDU 3336Count the string(KMP:串前缀匹配自身+DP):简单的前缀匹配+DP问题。解题报告!
POJ3080 Blue Jeans(KMP:最长连续公共子序列):求多个字符串的最长公共连续子串,如果存在多个长度相同的就输出字典序最小的那个。解题报告!
POJ3450 Corporate Identity(KMP:最长连续公共子序列):类似于上一题。解题报告!
HDU2087剪花布条(KMP:贪心):需要贪心的从左到右有考虑。解题报告!
HDU2203亲和串(KMP:循环移位):给你串T和串P问你串P是否能和T循环移K位后的串匹配。解题报告!
HDU1867 A + B for you again(KMP:后缀与前缀):字符串合并。解题报告!
KMP前缀与后缀
HDU 2594Simpsons’ HiddenTalents(KMP:后缀与前缀):KMP的思想。解题报告!
POJ 2752Seek the Name, Seekthe Fame(KMP:后缀与前缀):找出所有前缀与后缀匹配的长度。解题报告!
求字符串最短循环节
UVA 1328Period(KMP:最短循环节):如何求字符串的循环节?解题报告!
POJ 2406Power Strings(KMP:找串循环节):求出最小循环节个数。解题报告!
HDU 3746Cyclic Nacklace(KMP:补齐循环节):要你补充完整一个串,使得该串有循环节构成。解题报告!
POJ 2185Milking Grid(KMP:循环节加强版):转换思维。解题报告!