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 ;
 }