数据结构实验之串一:KMP简单应用
Time Limit: 1000MS Memory limit: 65536K
题目描述
给定两个字符串string1和string2,判断string2是否为string1的子串。
输入
输入包含多组数据,每组测试数据包含两行,第一行代表string1(长度小于1000000),第二行代表string2(长度小于1000000),string1和string2中保证不出现空格。
输出
对于每组输入数据,若string2是string1的子串,则输出string2在string1中的位置,若不是,输出-1。
示例输入
abc a 123456 45 abc ddd
示例输出
1 4 -1
///ACcode
#include <iostream>
#include <cstdio>
using namespace std;
const int maxn=10000010;
typedef struct
{
char *top,*base;
int Size;
int len;
} String;
void Creat(String &s) ///建串
{
s.base=new char[maxn];
s.top=s.base;
s.Size=maxn;
s.len=0;
}
void push(String &s,char c) ///入串...函数名有点眼熟...
{
s.top++;
*s.top=c;
s.len++;
}
void add(String &s,char st[]) ///为了使 两个串 下标从零开始
{
int i;
for (i=0; st[i]; i++)
{
push(s,st[i]);
}
}
void getnext(String &s,int next[]) /// 相当于KMP函数 其实就是和本身比较来寻找next[]数组;建议先看下面的KMP函数
{
int i=1,j=0; ///j从 1 开始的话 next[2]=2;显然不符合next[]数组....
next[1]=0;
while(i<=s.len)
{
if (j==0||s.base[i] == s.base[j])
{
i++;
j++;
next[i]=j; /// 动态演示中很明显看出 前 j 个字符已经匹配
}
else
{
j=next[j]; /// 失配时 j“跳”到next[j]位置(也就是说next[j]之前的位置已经匹配)
}
}
}
void kmp(String &st,String &pt,int next[]) ///kmp ,st主串 pt模式串
{
int i=1,j=1;
while(i<=st.len && j<=pt.len)
{
if(j==0||st.base[i] == pt.base[j])
{
i++;
j++;
}
else
{
j=next[j]; /// 失配时 j“跳”到next[j]位置(也就是说next[j]之前的位置已经匹配)
}
}
if (j>pt.len) ///说明模式串pt匹配完成
{
cout<<i-pt.len<<endl;
}
else ///串st遍历先完了,则pt不是st子串
{
cout<<"-1"<<endl;
}
}
void print(String &s) ///打印
{
int i;
for (i=1; i<=s.len; i++)
{
cout<<s.base[i];
}
cout<<endl;
}
int main()
{
int next[maxn];
char st[maxn],pt[maxn]; ///先数组存 主串 , 模式串
while(cin>>st>>pt)
{
String s,p; ///主串 , 模式串
Creat(s);
Creat(p);
add(s,st);
add(p,pt);
getnext(p,next);
kmp(s,p,next);
}
return 0;
}
如下图
↓ ↓ ↓ ↓当求next[4]时,前面绿色表示已经匹配,即a匹配 next[i]就等于j 也就是next[i-1]+1
↓ ↓ ↓ ↓同理 ,next[6]前面有 两个绿色表示已经匹配 next[i]=j (next[i-1]+1)
next[i]失配 模式串(下面那个串)j=next[j],继续寻找匹配的位置......
↓ ↓ ↓ ↓同理 next[8]也如此求 next[]等于1 表示紧挨着前面没有匹配成功的 2 表示1个 以此类推..