如题,给出两个字符串s1和s2,其中s2为s1的子串,求出s2在s1中所有出现的位置。
为了减少骗分的情况,接下来还要输出子串的前缀数组next。如果你不知道这是什么意思也不要问,去百度搜[kmp算法]学习一下就知道了。
输入格式:
第一行为一个字符串,即为s1(仅包含大写字母)
第二行为一个字符串,即为s2(仅包含大写字母)
输出格式:
若干行,每行包含一个整数,表示s2在s1中出现的位置
接下来1行,包括length(s2)个整数,表示前缀数组next[i]的值。
输入样例#1:
ABABABC
ABA
输出样例#1:
1
3
0 0 1
时空限制:1000ms,128M
数据规模:
设s1长度为N,s2长度为M
对于30%的数据:N<=15,M<=5
对于70%的数据:N<=10000,M<=100
对于100%的数据:N<=1000000,M<=1000
样例说明:
所以两个匹配位置为1和3,输出1、3
裸KMP
思路来源http://www.cnblogs.com/c-cloud/p/3224788.html
1 #include<iostream>
2 #include<cstdio>
3 #include<cstring>
4 using namespace std;
5 const int MAXN=1000001;
6 char a[MAXN],b[MAXN];
7 int la,lb;
8 int p[MAXN];
9 void makep()
10 {
11 int j=0;// 记录最长公共前后缀的长度;
12 p[0]=0;// 第一个数的前缀和后缀都是空集,故为0
13 for(int i=1;i<lb;i++)// 第一个是0,所以从1开始,半闭半开区间
14 {
15 while(j>0&&b[j]!=b[i])//当出现不匹配的情况的的话,就让j从p数组中取一个能继续匹配的元素
16 j=p[j-1];// 因为对j的最后一次处理是j++,而p是从0开始,故j=p[j-1]
17 if(b[j]==b[i])//如果相同的话就让j++ 自己可以举例子试试
18 j++;
19 p[i]=j;// 将每一个元素的最长公共前后缀的长度记录下来,方便以后处理
20 }
21 }
22 void KMP()
23 {
24 int j=0;// 已经匹配的最长的长度
25 for(int i=0;i<la;i++)
26 {
27 while(j>0&&b[j]!=a[i])// 当将要匹配的元素不想同且找到相同的匹配元素的时候
28 j=p[j-1];// j==0的时候p[-1]不存在,所以j>0
29 if(b[j]==a[i])// 同理,当将要匹配的元素相同的时候,j++
30 j++;
31 if(j==lb)// 找到一个匹配
32 printf("%d\n",i+1-lb+1);
33 // 因为i,lb都是从0开始计数的,so+1
34 }
35 }
36 int main()
37 {
38 scanf("%s %s",a,b);
39 la=strlen(a);lb=strlen(b);//统计两个字符串的长度
40 makep();// 构造p数组
41 KMP();//
42 for(int i=0;i<lb;i++)printf("%d ",p[i]);
43 return 0;
44 }