问题描述

给你两个字符串 haystack 和 needle ,请你在 haystack 字符串中找出 needle 字符串的第一个匹配项的下标(下标从 0 开始)。如果 needle 不是 haystack 的一部分,则返回  -1 。

示例 1:

输入:haystack = “sadbutsad”, needle = “sad”
输出:0
解释:“sad” 在下标 0 和 6 处匹配。
第一个匹配项的下标是 0 ,所以返回 0 。

示例 2:

输入:haystack = “leetcode”, needle = “leeto”
输出:-1
解释:“leeto” 没有在 “leetcode” 中出现,所以返回 -1 。

提示:

  • 1 <= haystack.length, needle.length <= 104
  • haystack 和 needle 仅由小写英文字符组成

核心思路

记字符串 haystack 的长度为 n,字符串 needle 的长度为 m。

我们记字符串 str=needle+#+haystack,即将字符串 needle 和 haystack 进行拼接,并用不存在于两串中的特殊字符 # 将两串隔开,然后我们对字符串 str 求前缀函数。

因为特殊字符 # 的存在,字符串 str 中 haystack 部分的前缀函数所对应的真前缀必定落在字符串 needle 部分,真后缀必定落在字符串 haystack 部分。当 haystack 部分的前缀函数值为 m 时,我们就找到了一次字符串 needle 在字符串 haystack 中的出现(因为此时真前缀恰为字符串 needle)。

实现时,我们可以进行一定的优化,包括:

我们无需显式地创建字符串 str。
为了节约空间,我们只需要顺次遍历字符串 needle、特殊字符 # 和字符串 haystack 即可。
也无需显式地保存所有前缀函数的结果,而只需要保存字符串 needle 部分的前缀函数即可。
特殊字符 # 的前缀函数必定为 0,且易知 π(i)m\pi(i)≤m(真前缀不可能包含特殊字符 #)。
这样我们计算 π(i) 时,j=π(π(π()1)1)j=\pi(\pi(\pi(…)-1)-1) 的所有的取值中仅有 π(i1)\pi(i-1) 的下标可能大于等于 m。我们只需要保存前一个位置的前缀函数,其它的 j 的取值将全部为字符串 needle 部分的前缀函数。
我们也无需特别处理特殊字符 #,只需要注意处理字符串 haystack 的第一个位置对应的前缀函数时,直接设定 j 的初值为 0 即可。
这样我们可以将代码实现分为两部分:

第一部分是求 needle 部分的前缀函数,我们需要保留这部分的前缀函数值。
第二部分是求 haystack 部分的前缀函数,我们无需保留这部分的前缀函数值,只需要用一个变量记录上一个位置的前缀函数值即可。当某个位置的前缀函数值等于 m 时,说明我们就找到了一次字符串 needle在字符串 haystack 中的出现(因为此时真前缀恰为字符串 needle,真后缀为以当前位置为结束位置的字符串 haystack 的子串),我们计算出起始位置,将其返回即可。

实现要点

对于长度为 m 的字符串 s,其前缀函数 π(i)(0i<m)\pi(i)(0≤i<m) 表示 s 的子串 s[0:i]s[0:i] 的最长的相等的真前缀与真后缀的长度。特别地,如果不存在符合条件的前后缀,那么 π(i)=0。其中真前缀与真后缀的定义为不等于自身的的前缀与后缀。

如何求解前缀函数

长度为 m 的字符串 s 的所有前缀函数的求解算法的总时间复杂度是严格 O(m) 的,且该求解算法是增量算法,即我们可以一边读入字符串,一边求解当前读入位的前缀函数。

为了叙述方便,我们接下来将说明几个前缀函数的性质:

  1. π(i)π(i1)+1\pi(i)≤\pi(i-1)+1

    • 依据 π(i)\pi(i) 定义得:s[0:π(i)1]=s[iπ(i)+1:i]s[0:\pi(i)-1]=s[i-\pi(i)+1:i]
    • 将两区间的右端点同时左移,可得:s[0:π(i)2]=s[iπ(i)+1:i1]s[0:\pi(i)-2]=s[i-\pi(i)+1:i-1]
    • 依据 π(i1)\pi(i-1) 定义得:π(i1)π(i)1\pi(i-1)≥\pi(i)-1,即 π(i)π(i1)+1\pi(i)≤\pi(i-1)+1
  2. 如果 s[i]=s[π(i1)]s[i]=s[\pi(i-1)],那么 π(i)=π(i1)+1\pi(i)=\pi(i-1)+1

    • 依据 π(i1)\pi(i-1) 定义得:s[0:π(i1)1]=s[iπ(i1):i1]s[0:\pi(i-1)-1]=s[i-\pi(i-1):i-1]
    • 因为 s[π(i1)]=s[i]s[\pi(i-1)]=s[i],可得 s[0:π(i1)]=s[iπ(i1):i]s[0:\pi(i-1)]=s[i-\pi(i-1):i]
    • 依据 π(i)\pi(i) 定义得:π(i)π(i1)+1\pi(i)≥\pi(i-1)+1,结合第一个性质可得 π(i)=π(i1)+1\pi(i)=\pi(i-1)+1

这样我们可以依据这两个性质提出求解 π(i)\pi(i) 的方案:找到最大的 j,满足 s[0:j1]=s[ij:i1]s[0:j-1]=s[i-j:i-1],且 s[i]=s[j]s[i]=s[j](这样就有 s[0:j]=s[ij:i]s[0:j]=s[i-j:i],即 π(i)=j+1\pi(i)=j+1)。

注意这里提出了两个要求:

  1. j 要求尽可能大,且满足 s[0:j1]=s[ij:i1]s[0:j-1]=s[i-j:i-1]

  2. j 要求满足 s[i]=s[j]s[i]=s[j]

π(i1)\pi(i-1) 定义可知:

s[0:π(i1)1]=s[iπ(i1):i1](1)s[0:\pi(i-1)-1]=s[i-\pi(i-1):i-1] (1)

那么 j=π(i1)j=\pi(i-1) 符合第一个要求。如果 s[i]=s[π(i1)]s[i]=s[\pi(i-1)],我们就可以确定 π(i)\pi(i)

否则如果 s[i]s[π(i1)]s[i]≠s[\pi(i-1)],那么 π(i)π(i1)\pi(i)≤\pi(i-1),因为 j=π(i)1j=\pi(i)-1,所以 j < π(i1)j < \pi(i-1),于是可以取 (1) 式两子串的长度为 j 的后缀,它们依然是相等的:s[π(i1)j:π(i1)1]=s[ij:i1]s[\pi(i-1)-j:\pi(i-1)-1]=s[i-j:i-1]

s[i]s[π(i1)]s[i]≠s[\pi(i-1)] 时,我们可以修改我们的方案为:找到最大的 j,满足 s[0:j1]=s[π(i1)j:π(i1)1]s[0:j-1]=s[\pi(i-1)-j:\pi(i-1)-1],且 s[i]=s[j]s[i]=s[j](这样就有 s[0:j]=s[π(i1)j:π(i1)]s[0:j]=s[\pi(i-1)-j:\pi(i-1)],即 π(i)=π(i1)+1\pi(i)=\pi(i-1)+1)。

注意这里提出了两个要求:

  1. j 要求尽可能大,且满足 s[0:j1]=s[π(i1)j:π(i1)1]s[0:j-1]=s[\pi(i-1)-j:\pi(i-1)-1]

  2. j 要求满足 s[i]=s[j]s[i]=s[j]

π(π(i1)1)\pi(\pi(i-1)-1) 定义可知 j=π(π(i1)1)j=\pi(\pi(i-1)-1) 符合第一个要求。如果 s[i]=s[π(π(i1)1)]s[i]=s[\pi(\pi(i-1)-1)],我们就可以确定 π(i)\pi(i)

此时,我们可以发现 j 的取值总是被描述为 π(π(π(1)1)\pi(\pi(\pi(…-1)-1) 的结构(初始为 π(i1)\pi(i-1))。于是我们可以描述我们的算法:设定 π(i)=j+1\pi(i)=j+1,j 的初始值为 π(i1)\pi(i-1)。我们只需要不断迭代 j(令 j 变为 π(j1)\pi(j-1))直到 s[i]=s[j]s[i]=s[j]j=0j=0 即可,如果最终匹配成功(找到了 jj 使得 s[i]=s[j]s[i]=s[j]),那么 π(i)=j+1\pi(i)=j+1,否则 \pi(i)=0。

code

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
class Solution {
public:
int strStr(string haystack, string needle) {
int n = haystack.size(), m = needle.size();
if (m == 0) {
return 0;
}
vector<int> pi(m);
for (int i = 1, j = 0; i < m; i++) {
while (j > 0 && needle[i] != needle[j]) {
j = pi[j - 1];
}
if (needle[i] == needle[j]) {
j++;
}
pi[i] = j;
}
for (int i = 0, j = 0; i < n; i++) {
while (j > 0 && haystack[i] != needle[j]) {
j = pi[j - 1];
}
if (haystack[i] == needle[j]) {
j++;
}
if (j == m) {
return i - m + 1;
}
}
return -1;
}
};
  • 时间复杂度:O(n+m)O(n+m),其中 n 是字符串 haystack 的长度,m 是字符串 needle 的长度。我们至多需要遍历两字符串一次。

  • 空间复杂度:O(m)O(m),其中 m 是字符串 needle 的长度。我们只需要保存字符串 needle 的前缀函数。