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中出现的次数。

【样例输入】

BAPC

AZA

AZAZAZA

VERDI

AVERDXIVYERDIAN

【样例输出】

思路：这是一道裸的KMP算法题，题目中模式串在第一行，主串在第二行，按顺序读入即可。如果用string存，可能会超时，因此用char数组存储，但是在存储时不要用char[0]来表示字符串长度，此时数据类型为char，范围最多只有0-255，因此要额外设int型变量来存字符串长度。用《数据结构 C语言版》清华出版的那本书上的KMP模板会超时，可能是因为 j 多余了。每次都先找到相匹配的位置，j 达到模式串长度时ans+1即可。

代码：

 #define _CRT_SECURE_NO_WARNINGS

 #include <iostream>

 #include <cstdio>

 #define MAXLENGTH 10000

 using namespace std;

 //typedef unsigned char SString[1000100];

 int nextP[];

 char s[], t[];

 int slen, tlen;

 void getNext(char*);

 int indexKMP(char*, char*, int);

 int getLen(char*);

 void getNext(char* T)

 {

     int i = , j = -;

     nextP[] = -;

     while (i < tlen)

     {

         if (j == - || T[i] == T[j])

         {

             ++i, ++j;

             if (T[i] != T[j])

                 nextP[i] = j;

             else

                 nextP[i] = nextP[j];

         }

         else

             j = nextP[j];

     }

 }

 int indexKMP(char* S, char* T, int pos) //从pos处开始匹配

 {

     int i = pos, j = , res = ;

     getNext(T);

     while (i < slen)

     {

         while (j != - && S[i] != T[j]) j = nextP[j];

         i++, j++;

         if (j == tlen)

             res++;

     }

     return res;

 }

 int getLen(char* T)

 {

     int len = ;

     for (int i = ; T[i] != '\0'; i++)

         len++;

     return len;

 }

 int main()

 {

     freopen("oulipo.in", "r", stdin);

     freopen("oulipo.out", "w", stdout);

     int x;

     cin >> x;

     while (x--)

     {

         scanf("%s%s", t, s);

         slen = getLen(s);

         tlen = getLen(t);

         cout << indexKMP(s, t, ) << endl;

     }

     return ;

 }