KMP --关于cogs1570 乌力波

题目链接：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 
 3 #include 
 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 }