N数之和问题:给定一个对于所有自变量$x_1, x_2, \cdots, x_n$都在一个有限的整数区间中,且严格单调递增的函数$f(x_1, x_2, \cdots, x_n)$。如何找到该函数的零点?
为了解决N数之和问题,我们首先考虑两数之和问题,我们从一个无任何附加信息的两和之和问题开始,通过不断添加新的信息,利用新的信息,来降低算法的复杂度。
问题一:
给定一个关于在相同有限区间$S$内的整数变量$x$和$y$的函数$f(x, y)$。找到$f(x,y)$的零点。
通过暴力搜索所有可能的状态空间找到答案。搜索空间的大小即$|S|^2$。通过遍历所有的$(x, y)$,找到所有使得$f(x, y)=0$的解。
1 | y 0 ... N |
问题二:
给定一个关于在相同有限区间$S$内的整数变量$x$和$y$的函数$f(x, y)$。对于方程$f(x, y) = 0$,有隐函数$y = h(x)$。找到$f(x,y)$的零点。
由于函数$y = h(x)$存在,在给定$x$的时候,我们可以找到对应的$y$。因此当我们在遍历$x$的所有可能取值时,可以用$h(x)$计算出其对应的$y$,并放在一个hashset中。每次到下一个$s$的可能取值时,因为$x$和$y$的取值范围是完全相同的,所以可以先在hashset中查看是否有该值,如果有,则找到了一个解。
或者当函数$y = h(x)$计算代价不大时,我们不必在hashset
中记录$x$对应的$y$,而是直接记录$x$。
1 | x 0 ... N |
问题三:
给定一个关于在相同有限区间$S$内的整数变量$x$和$y$,且严格单调递增的函数$f(x, y)$,找到$f(x,y)$的零点。
考虑函数的单调性:
因为$f(x, y)$严格单调递增,所以我们有
- $\forall t > x, f(t, y) > f(x, y)$
- $\forall t < y, f(x, t) < f(x, y)$
为了利用上述两条性质,我们将有限区间$S$所有可能的取值进行排序。之后让$x$等于区间中的最小值$i$。让$y$等于区间中的最大值$j$。判断$f(i, j)$和0的大小关系:
- 如果$f(i, j) < 0$,则$\forall t < j, f(i, t) < f(i, j) < 0$。因此当$x = i$时,$\forall y \in S, f(i, y) < 0$。我们可以在搜索空间中划除$i$对应的这一行。
- 如果$f(i, j) > 0$,则$\forall t > i, f(t, j) > f(i, j) > 0$。因此当$y = j$时,$\forall x \in S, f(x, j) > 0$。我们可以在搜索空间中划去$j$对应的这一列。
- 如果$f(i, j) = 0$,我们找到了一个结果。但根据严格单调递增性,我们不会在$i$,$j$所在的行列再找到任何其他一对数字使得$f(x, y) = 0$。我们可以在搜索空间划去$i$对应的这一行和$j$对应的这一列。
之后我们得到了一个新的搜索空间。我们继续从这个搜索空间的右上角(即$i$取最小值,$j$取最大值)重复上面的算法,直到搜索空间为空。
1 | start start |
问题三的一个典型例题,未给出函数表达式,仅强调了单调性:LC1237 找出给定方程的正整数解
给定函数严格单调递增
$f(x, y) < f(x + 1, y)$
$f(x, y) < f(x, y + 1)$
在$0 \leq x \leq 1000$ 和 $0 \leq y \leq 1000$的条件下,找到使得$f(x, y) - z = 0$的所有可能的解。
思路:
搜索空间为$\{(x,y)| 0 \leq x \leq 1000, 0 \leq y \leq 1000, x \in \mathbb{Z}, y \in \mathbb{Z}\}$,如下图所示。1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24 (a, b)
^
|
j 0 ... 1000 | j 0 ... 1000 j 0 ... 1000
i +-+-+-+-+---+ | i +-+-+-+-+---+ i +-----------+
0 | | | | | |S+----+ 0 | | | | | |X| 0 |X|X|X|X|X|X|
+-----------+ +-----------+ +-----------+
| | | | | | | | | | | | |X| | | | | | | |
+-----------+ +-----------+ +-----------+
. | | | | | | | . | | | | | |X| . | | | | | | |
. +-----------+ . +-----------+ . +-----------+
. | | | | | | | . | | | | | |X| . | | | | | | |
+-----------+ +-----------+ +-----------+
| | | | | | | | | | | | |X| | | | | | | |
+-----------+ +-----------+ +-----------+
1000 | | | | | | | 1000 | | | | | |X| 1000 | | | | | | |
+-+---+-+-+-+ +-+-+---+---+ +-+-+-+---+-+
| ^ ^
| f(a, b) - z > 0 | |
+------------------------+ |
| |
| f(a, b) - z < 0 |
+-------------------------------------------------+
根据函数严格单调递增的性质,我们知道
$f(a, b) - z < 0$,且$b$一定时,$\Rightarrow \forall x \in [a, 1000], f(x, b) < 0$
$f(a, b) - z > 0$,且$a$一定时,$\Rightarrow \forall x \in [b, 1000], f(a, x) > 0$
$f(a, b) - z = 0$,当$a$一定时,$\Rightarrow \forall x \in (b, 1000], f(a, x) \neq 0$
$f(a, b) - z = 0$,当$b$一定时,$\Rightarrow \forall x \in (a, 1000], f(x, b) \neq 0$
从右上角$(0, 1000)$的位置开始搜索。根据$f(a, b) - z$和$0$的大小关系,可以划去一行和(或)者一列。
- 即当$f(a, b) - z > 0$时,我们可以划去当前$a$在搜索空间中对应的这一行数字,因为函数单调递增,增加$a$只会使得$f(a, b) - z$变得更大。
- 当$f(a, b) - z < 0$时,我们可以划去$b$在搜索空间中对应的这一列数字。
- 当$f(a, b) - z = 0$时,我们同时划去$(a, b)$在搜索空间对应这一行和这一列。
1 | class Solution { |
问题4:
给定一个关于在相同有限区间$S$内的整数变量$x$和$y$,且严格单调递增的函数$f(x, y)$,找到$f(x,y) < 0$的所有解。
当我们找到一个点使得$f(a, b) < 0$时,因为单调性,所有在该点同一列,且该点上方的点$(t, b), t < a$都有$f(t, b) < 0$。所有在该点同一行,且在该点左侧的点$(a, t), t < b$都有$f(a, t) < 0$
同理,当我们找到一个点使得$f(a, b) \geq 0$时,其下方和右侧所有的点$(m, n)$也都满足$f(m, n) \geq 0$。1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20 ^
|
j 0 | N j 0 N
i +-+-+--|--+-+ i +-+-+-+-+-+-+
0 | | | ||| | | 0 | | | | | | |
+------|----+ +-----------+
| | | ||| | | | | | | | | |
+------|----+ +-----------+
<---------X| | | | | | |X--------->
+-----------+ +------|----+
| | | | | | | | | | ||| | |
+-----------+ +------|----+
| | | | | | | | | | ||| | |
+-----------+ +------|----+
N | | | | | | | N | | | ||| | |
+-+-+-+-+-+-+ +-+-+--|--+-+
|
|
v
我们仍然选择从右上角(或者左下角)开始进行搜索。 我们不断地移动$x$,初始化$y$为0。
我们通过判断$f(x, y) \geq 0$来找到这一列的第一个有效解的$y$。我们记录下这个$y$的位置,在下一次循环开始后,因为$x$增加,$y$必不可能减少。我们无需从$N$开始重新搜索,只需要从上一次的$y$位置开始搜索,找到第一个有效位置即可。1
2
3
4
5
6
7
8
9
10
11int cnt = 0;
int y = N;
for(int x = 0; x <= N; ++x)
{
// update y to its first valid position in this row
while(f(x, y) >= 0)
{
--y;
}
cnt += y;
}
1 | j 0 N |
问题4的典型例题为LC719 找出第 k 小的距离对
给定一个整数数组,返回所有数对之间的第 k 个最小距离。一对$(A, B)$的距离被定义为 A 和 B 之间的绝对差值。
我们可以用二分搜索找到第k个最小距离。注意到该整数数组中距离的最大值为数组的最大值减去数组的最小值。最小值可能为0。为们在距离上进行二分搜索:
- 设$g(d)$为在该数组上,小于$d$的距离对的数目。
- 维护二分搜索区间,使得
- 对于[begin, first)所在区间的任意值$t$,$g(t) < k$
- 对于[last, end)所在区间的任意值$t$,$g(t) >= k$
- first < last
最后first == last
所停在的位置即为结果。
该题目的关键点在于函数$g$的实现。注意到如果我们将该数组排序,并总是用较大的数减去较小的数,就不用考虑绝对值的问题。对于一对点$(x, y)$,他们之间的距离与$t$的差$f(x, y)$是关于x递减,关于y递增的函数。参照问题四,如果一点$f(x, y) < 0$,那么他左侧和他下方的点都满足$f(x, y) < 0$。如果一点$f(x, y) > 0$,那么他上方和他右侧的点也都满足$f(x, y) > 0$。我们考虑从左上角开始搜索。移动j
并不断通过$f(x,y)$和0的关系判断移动i
。这里和问题4不同的地方在于,数组中的每个数只能选取一次,且距离函数是对称的。有一半的搜索空间是无效的,因此我们每次累加cnt += j - i
个。
1 | auto dist_cnt(vector<int>& nums, int k) -> size_t |
1 | auto dist_cnt(vector<int>& nums, int k) -> size_t |
两数之和
给定一个整数数组 nums 和一个目标值 target,请你在该数组中找出和为目标值的那两个整数,并返回他们的数组下标。
改变叙述方式,该题目即对于属于整数数组中的$x$和$y$,找到函数$f(x, y) = x + y - t$的一个零点。
思路一:
- 对于$y$,我们可以求得其对应的$x$为$x = t - y$。因此我们可以使用问题二的思路,通过一个hashmap来解决。
- 从左到右遍历数组,在遍历过程中,对于每个元素
a
,查询unordered_map
中是否有target - a
。如果有,就返回结果。 - 将
(a, a_i)
放入unordered_map
中。 - 即使遍历一次,仍然考虑到了所有的
(i, j)
的可能性。
- 从左到右遍历数组,在遍历过程中,对于每个元素
1 | class Solution { |
我们可以发现,函数表达式$f(x, y) = x + y - t$对于$x$和$y$都是严格单调递增的。因此我们可以采用问题三的思路。先将数组排序,之后从$x$取最小值,$y$取最大值的位置开始搜索。即双指针法。然而LC1 两数之和这个题目,要求返回原来的元素下标。因此我们必须使用一个map
来建立有序的数值-下标对应关系。或者使用一个unordered_map
来记录对应关系后再进行排序。这两种方法都不如思路一简单,且效率较低。最适合利用单调性的是输入的数组已经排序好了的167 两数之和 II - 输入有序数组
当数组有序时
我们设置两个指针,一个位于数组第一个元素front
,一个位于数组最后一个元素back
。
- 保证
front < back
- 如果
target < *front + *back
,说明*first
小了,将first
往前移动一位。 - 如果
target > *front + *back
,说明*back
大了,将back
往后移动一位。 - 如果
target == *front + *back
,得到结果返回。
- 如果
1 | j 0 1 2 3 4 5 j 0 1 2 3 4 5 j 0 1 2 3 4 5 |
算法开始时,我们搜索位置$(0,5)$。因为数组是有序的,所以第一行中位置$(0,0)$到$(0,4)$的值都要比$(0,5)$小。最后一列中$(1,5)$到$(5,5)$都要比$(0,5)$大。
- 如果我们发现target $ > (0,5)$,那么第一行中所有的值因为必然小于$(0, 5)$,也就小于target,就不必再搜索了。因此我们向前移动
i
- 同理,如果$(0,5)$比target大,那么最后一列的内容就不必搜索了。因此我们向前移动
j
- 之后重复上述过程,即可找到
target
。
解释:一张图告诉你 O(n) 的双指针解法的本质原理(C++/Java)- nettee, Leetcode
1 | class Solution { |
接下来是一个在动态的数组上寻找两数之和的问题:170. 两数之和 III - 数据结构设计
不同于之前的两数之和问题,该题目中可能出现重复的数字,也可能出现多个可能解。
选择一:排序数组。每次add
操作必须要保证数组有序,因此会有$O(N)$的时间复杂度。find
操作可以使用双指针。
选择二:排序链表。add
不能进行二分查找,$O(N)$的时间复杂度。find
可以使用双指针。内存局部性差。
选择三:hashmap。add
期望$O(1)$,查找$O(N)$。
选择四:数组不排序。每次find
时再排序。之后双指针查找。
可以看到,如果想要利用单调性来解决该问题,就必须得维持数组有序,代价较高。在已知表达式的情况下,显然是hashmap
法要更简洁高效。
注意:因为存在重复元素,我们需要统计元素的个数。当hashmap
中存在一个元素a
,2 * a
和所查询元素相等时,如果a
只出现1次,要返回false
。
1 | class TwoSum { |
在BST上的两数之和问题,可以尝试设计BST的迭代器来解决。
653. 两数之和 IV - 输入 BST
思路一:双指针。为了实现在BST上的双指针,我们需要一个BST的正向迭代器,和一个反向迭代器。迭代其的实现见173. 二叉搜索树迭代器
1 | /** |
思路二:花费很大经历去写BST的迭代器比较麻烦。直接中序便利整个BST,将结果放入一个数组里再用双指针。
思路三:无视BST的有序性质,直接采用LC1 两数之和的unordered_map
来做。
按照问题二的思路:
我们不能再使用无序的hashset来记录已经遍历的结果。因为对于我们当前的数$x$,我们要找到以前遍历过的数字中,小于$K - x$且最最大的。因此我们需要使用一个有序容器set
来记录之前遍历的数字。我们以greater<int>
作为set
的比较函数,这样我们在使用upper_bound
时,就可以找到首个小于$K - x$的元素了。
1 | class Solution { |
按照问题三的思路,我们将数组排序后利用单调性:
当*l + *r >= K
时,因为数组有序,在[l, r)之间的任意的位置p都会使得*p + *r >= K
,因此--r
。当*l + *r < K
时,记录下此时的*l + *r
,[begin, l)之间任意的位置p,都会使得*p + *r < *l + *r
,绝不比*l + *r
更接近K
。因此++l
。
1 | class Solution { |
利用数组的前缀和,将问题转化为两数之和问题:
- 寻找两个前缀和
a
,b
,使得b - a == k
或寻找一个数字a
,使得a == k
- 因为有重复元素,因此按照[LC1 两数之和]的去重
hashmap
法来做。
1 | class Solution { |
使用前缀和。因为要求总和为零,所以表达式已知。使用一个hashmap从前向后遍历。要求直到不存在这样的序列为止。因此在遍历过程中递归。
1 | class Solution { |
类似LC560 和为K的子数组,利用前缀和,将问题转化为两数之和。
- 因为不是寻找相等的,而是寻找在一个区间内的,我们使用multiset来代替hashmap
- 遍历数组,计算每个位置对应的前缀和
- 如果该前缀和
sum
本身位于[l, u]
之间,那么该位置到数组头一个元素构成的区间满足条件,++cnt
(也可以在multiset
中放置一个guard元素0
来省略这一步) - 之后在
multiset
中寻找能够使得sum - b
满足位于[l, u]
之间的b
。b
应该满足sum - b >= l && sum - b <= u + 1
。即寻找lower_bound(sum - l + 1)
到lower_bound(sum - u)
之间元素的个数。
- 如果该前缀和
1 | class Solution { |
1 | class Solution { |
N数之和问题
固定a
的位置,移动b
和c
来当作两数之和来做。
- 找到所有可能的三元组:每当找到一个可能的三元组时,因为数组有序,我们可以确定地认为:
- 当
b
一定时,改变c
到(b, c)
之间的任何一个元素t
,都不会有a + b + t == 0
,因此++b
- 同理,当
c
一定时,改变b
到(b, c)
之间的任何一个元素,都不会有a + b + t == 0
,因此--c
- 当
- 去除重复元素:
- 每当
a
,b
或c
移动后,都和移动前的位置比较一下。 b
和c
仅需要在确定有结果后,才需要比较。因为如果a + b + c != 0
,那么即使有重复也无所谓。
- 每当
1 | class Solution { |
排序后,固定i
指针,移动j
和k
指针。
- 如果
*i + *j + *k >= target
,那么无论j
移动到[j, k)
之间的任何一个位置,也都有*i + *j + *k >= target
。 - 如果
*i + *j + *k < target
,那么无论k
移动到[j, k)
之间的任何一个位置,都有*i + *j + *k < target
,因此计数器应增加k - j
个。
1 | class Solution { |
仍然按照三数之和的代码进行,但是当sum != target
时要差最小的sum
。当sum == target
时,直接返回sum
。
1 | class Solution { |
在三数之和外面再套一层循环即可。
1 | class Solution { |
- 用一个hashmap去记录AB数组中的所有可能的和
- 遍历CD数组,寻找CD数组的和的相反数在hashmap中的个数
1 | class Solution { |