题目链接:http://cogs.pro:8081/cogs/problem/problem.php?pid=vQzXJkgWa
【题目描述】
法国作家乔治·佩雷克(Georges Perec,1936-1982)曾经写过一本书,《敏感字母》(La disparition),全篇没有一个字母‘e’。他是乌力波小组(Oulipo Group)的一员。下面是他书中的一段话:
Tout avait Pair normal, mais tout s’affirmait faux. Tout avait Fair normal, d’abord, puis surgissait l’inhumain, l’affolant. Il aurait voulu savoir où s’articulait l’association qui l’unissait au roman : stir son tapis, assaillant à tout instant son imagination, l’intuition d’un tabou, la vision d’un mal obscur, d’un quoi vacant, d’un non-dit : la vision, l’avision d’un oubli commandant tout, où s’abolissait la raison : tout avait l’air normal mais…
佩雷克很可能在下面的比赛中得到高分(当然,也有可能是低分)。在这个比赛中,人们被要求针对一个主题写出甚至是意味深长的文章,并且让一个给定的“单词”出现次数尽量少。我们的任务是给评委会编写一个程序来数单词出现了几次,用以得出参赛者最终的排名。参赛者经常会写一长串废话,例如500000个连续的‘T’。并且他们不用空格。
因此我们想要尽快找到一个单词出现的频数,即一个给定的字符串在文章中出现了几次。更加正式地,给出字母表{'A','B','C',...,'Z'}和两个仅有字母表中字母组成的有限字符串:单词W和文章T,找到W在T中出现的次数。这里“出现”意味着W中所有的连续字符都必须对应T中的连续字符。T中出现的两个W可能会部分重叠。
【输入格式】
输入包含多组数据。
输入文件的第一行有一个整数,代表数据组数。接下来是这些数据,以如下格式给出:
第一行是单词W,一个由{'A','B','C',...,'Z'}中字母组成的字符串,保证1<=|W|<=10000(|W|代表字符串W的长度)
第二行是文章T,一个由{'A','B','C',...,'Z'}中字母组成的字符串,保证|W|<=|T|<=1000000。
【输出格式】
对每组数据输出一行一个整数,即W在T中出现的次数。
【样例输入】
3
BAPC
BAPC
AZA
AZAZAZA
VERDI
AVERDXIVYERDIAN
【样例输出】
1
3
0
思路:这是一道裸的KMP算法题,题目中模式串在第一行,主串在第二行,按顺序读入即可。如果用string存,可能会超时,因此用char数组存储,但是在存储时不要用char[0]来表示字符串长度,此时数据类型为char,范围最多只有0-255,因此要额外设int型变量来存字符串长度。用《数据结构 C语言版》清华出版的那本书上的KMP模板会超时,可能是因为 j 多余了。每次都先找到相匹配的位置,j 达到模式串长度时ans+1即可。
代码:
1 #define _CRT_SECURE_NO_WARNINGS
2 #include <iostream>
3 #include <cstdio>
4 #define MAXLENGTH 10000
5 using namespace std;
6 //typedef unsigned char SString[1000100];
7 int nextP[1000100];
8 char s[1000100], t[1000100];
9 int slen, tlen;
10
11 void getNext(char*);
12 int indexKMP(char*, char*, int);
13 int getLen(char*);
14
15 void getNext(char* T)
16 {
17 int i = 0, j = -1;
18 nextP[0] = -1;
19 while (i < tlen)
20 {
21 if (j == -1 || T[i] == T[j])
22 {
23 ++i, ++j;
24 if (T[i] != T[j])
25 nextP[i] = j;
26 else
27 nextP[i] = nextP[j];
28 }
29 else
30 j = nextP[j];
31 }
32 }
33
34 int indexKMP(char* S, char* T, int pos) //从pos处开始匹配
35 {
36 int i = pos, j = 0, res = 0;
37 getNext(T);
38 while (i < slen)
39 {
40 while (j != -1 && S[i] != T[j]) j = nextP[j];
41 i++, j++;
42 if (j == tlen)
43 res++;
44 }
45 return res;
46 }
47
48 int getLen(char* T)
49 {
50 int len = 0;
51 for (int i = 0; T[i] != '\0'; i++)
52 len++;
53 return len;
54 }
55 int main()
56 {
57 freopen("oulipo.in", "r", stdin);
58 freopen("oulipo.out", "w", stdout);
59 int x;
60 cin >> x;
61 while (x--)
62 {
63 scanf("%s%s", t, s);
64 slen = getLen(s);
65 tlen = getLen(t);
66 cout << indexKMP(s, t, 0) << endl;
67 }
68 return 0;
69 }