POJ3461 Oulipo

嘟嘟嘟

今天又学了一遍kmp，总算是理解了。虽然这就是一道板子题，但我还想叨叨一会儿。

kmp匹配的大致过程就是对于主串A，我们总想让A[i - j + 1 ~ i]和模式串B的最长前缀B[j]相等。匹配的时候如果A[i]和B[j + 1]不等，就把j = f[j]，因为f[j]是保证了A[j - k + 1 ~ j]和B[k]相等的最大的k，那么可以得出对于A[i - 1]，B和A[i - 1]的后缀相等的前缀中，最长的是j, 第二长的是f[j]，第三长的是f[f[j]]……所以到了A[i]，我们就从最长的开始试，如果B[j + 1] != A[i]的话，再试第二长的，看B[f[j] + 1]是否等于A[i]，不行的话再看第三长的……所以说，我感觉kmp就是一种贪心……

构造f数组和匹配几乎一样，这里说一点，就是f[1] = 0，然后我们从f[2]开始构造。为啥咧？因为一个串本身也是自己的前缀，如果从1开始，那么f[1] = 1,然后f[2]必定等于f[1] + 1，结果f数组就变成了f[1] = 1,f[2] = 2, f[3] = 3……很明显gg了。

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<cstdlib>
#include<cctype>
#include<stack>
#include<queue>
#include<vector>
using namespace std;
#define enter puts("")
#define space putchar(' ')
#define Mem(a, x) memset(a, x, sizeof(a))
#define rg register
typedef long long ll;
typedef double db;
const int INF = 0x3f3f3f3f;
const db eps = 1e-8;
const int maxn = 1e4 + 5;
inline ll read()
{
  ll ans = 0;
  char ch = getchar(), las = ' ';
  while(!isdigit(ch)) las = ch, ch = getchar();
  while(isdigit(ch)) ans = ans * 10 + ch - '0', ch = getchar();
  if(las == '-') ans = -ans;
  return ans;
}
inline void write(ll x)
{
  if(x < 0) putchar('-'), x = -x;
  if(x >= 10) write(x / 10);
  putchar(x % 10 + '0');
}

char s[maxn], t[maxn * 100];
int f[maxn], cnt = 0;

void init()
{
  int m = strlen(s + 1);
  for(int i = 2, j = 0; i <= m; ++i)
    {
      while(j && s[j + 1] != s[i]) j = f[j];
      if(s[j + 1] == s[i]) j++;
      f[i] = j;
    }
}
void kmp()
{
  cnt = 0;
  int m = strlen(s + 1), n = strlen(t + 1);
  for(int i = 1, j = 0; i <= n; ++i)
    {
      while(j && s[j + 1] != t[i]) j = f[j];
      if(s[j + 1] == t[i]) j++;
      if(j == m) cnt++, j = f[j];
    }
}

int main()
{
  int T = read();
  while(T--)
    {
      scanf("%s
%s", s + 1, t + 1);
      init();
      kmp();
      write(cnt); enter;
    }
  return 0;
}