前段时间我有一次有趣的面试经历。问题一开始很简单:
Q1:我们有一个袋子,里面有数字1,2,3,…,100。每个数字恰好出现一次,所以有100个数字。现在从袋子里随机抽取一个数字。找到丢失的号码。
当然,我以前听过这个面试问题,所以我很快就回答了这个问题:
A1:嗯,1 + 2 + 3 +…+ N的和是(N+1)(N/2)(参见维基百科:等差级数的和)。当N = 100时,和是5050。
因此,如果所有的数字都在袋子里,总和将恰好是5050。因为少了一个数,总和就会小于这个数,差的就是这个数。所以我们可以在O(N)时间和O(1)空间中找到这个缺失的数。
在这一点上,我认为我做得很好,但突然间,问题发生了意想不到的转变:
这是正确的,但是如果少了两个数字,你会怎么做?
我以前从未见过/听过/考虑过这种变化,所以我很恐慌,无法回答这个问题。面试官坚持要知道我的思考过程,所以我提到,也许我们可以通过与预期产品进行比较来获得更多信息,或者在从第一次传递中收集到一些信息后再进行第二次传递,等等,但我真的只是在黑暗中拍摄,而不是真正有一个明确的解决方案的路径。
面试官试图鼓励我说,有第二个方程确实是解决问题的一种方法。在这一点上,我有点不安(因为事先不知道答案),并问这是一种通用的(阅读:“有用的”)编程技术,还是只是一个技巧/答案。
面试官的回答让我惊讶:你可以把这个技巧概括为3个缺失的数字。事实上,你可以推广它来找到k个缺失的数。
Qk:如果袋子里少了k个数字,你如何有效地找到它?
这是几个月前的事了,我还不明白这个技巧是什么。显然有一个Ω(N)的时间下限,因为我们必须扫描所有的数字至少一次,但面试官坚持认为,解决技术的时间和空间复杂度(减去O(N)次输入扫描)定义为k而不是N。
所以问题很简单:
如何解决Q2?
你会如何解决Q3?
如何求解Qk?
澄清
Generally there are N numbers from 1..N, not just 1..100.
I'm not looking for the obvious set-based solution, e.g. using a bit set, encoding the presence/absence each number by the value of a designated bit, therefore using O(N) bits in additional space. We can't afford any additional space proportional to N.
I'm also not looking for the obvious sort-first approach. This and the set-based approach are worth mentioning in an interview (they are easy to implement, and depending on N, can be very practical). I'm looking for the Holy Grail solution (which may or may not be practical to implement, but has the desired asymptotic characteristics nevertheless).
当然,你必须以O(N)为单位扫描输入,但你只能捕获少量的信息(用k而不是N定义),然后必须以某种方式找到k个缺失的数字。
基于数字和的解决方案的问题是,它们没有考虑到存储和处理具有大指数的数字的成本……在实践中,为了使它适用于非常大的n,将使用大数库。我们可以分析这些算法的空间利用率。
我们可以分析sdcvvc和Dimitris Andreou算法的时间和空间复杂度。
储存:
l_j = ceil (log_2 (sum_{i=1}^n i^j))
l_j > log_2 n^j (assuming n >= 0, k >= 0)
l_j > j log_2 n \in \Omega(j log n)
l_j < log_2 ((sum_{i=1}^n i)^j) + 1
l_j < j log_2 (n) + j log_2 (n + 1) - j log_2 (2) + 1
l_j < j log_2 n + j + c \in O(j log n)`
所以l_j \in \ (j log n)
所使用的总存储空间:\sum_{j=1}^k l_j \in \Theta(k^2 log n)
使用的空间:假设计算a^j需要ceil(log_2 j)时间,总时间:
t = k ceil(\sum_i=1^n log_2 (i)) = k ceil(log_2 (\prod_i=1^n (i)))
t > k log_2 (n^n + O(n^(n-1)))
t > k log_2 (n^n) = kn log_2 (n) \in \Omega(kn log n)
t < k log_2 (\prod_i=1^n i^i) + 1
t < kn log_2 (n) + 1 \in O(kn log n)
总使用时间:\Theta(kn log n)
如果这个时间和空间是令人满意的,您可以使用一个简单的递归
算法。让b !I是袋子里的第I个元素,n个之前的数字
移除次数,k是移除次数。在Haskell语法中…
let
-- O(1)
isInRange low high v = (v >= low) && (v <= high)
-- O(n - k)
countInRange low high = sum $ map (fromEnum . isInRange low high . (!)b) [1..(n-k)]
findMissing l low high krange
-- O(1) if there is nothing to find.
| krange=0 = l
-- O(1) if there is only one possibility.
| low=high = low:l
-- Otherwise total of O(knlog(n)) time
| otherwise =
let
mid = (low + high) `div` 2
klow = countInRange low mid
khigh = krange - klow
in
findMissing (findMissing low mid klow) (mid + 1) high khigh
in
findMising 1 (n - k) k
使用的存储:O(k)用于列表,O(log(n))用于堆栈:O(k + log(n))
该算法更直观,具有相同的时间复杂度,占用的空间更少。
我们可以用O(log n)来做Q1和Q2。
假设我们的存储芯片由n个试管阵列组成。试管中的数字x用x毫升化学液体表示。
假设我们的处理器是一束激光。当我们点燃激光时,它垂直穿过所有的管子。每次它通过化学液体,光度就降低1。在某毫升处通过光是O(1)的运算。
现在如果我们在试管的中间点上激光就会得到光度的输出
等于预先计算的值(当没有数字缺失时计算),则缺失的数字大于n/2。
如果我们的输出更小,那么至少有一个小于n/2的数字缺失。我们也可以检查光度是否降低了1或2。如果它减少1,那么一个缺失数小于n/2,另一个大于n/2。如果它减2,那么两个数都小于n/2。
我们可以一次又一次地重复上述过程,缩小我们的问题域。在每一步中,我们将定义域缩小一半。最后我们可以得到结果。
值得一提的是并行算法(因为它们很有趣),
sorting by some parallel algorithm, for example, parallel merge can be done in O(log^3 n) time. And then the missing number can be found by binary search in O(log n) time.
Theoretically, if we have n processors then each process can check one of the inputs and set some flag that identifies the number(conveniently in an array). And in the next step each process can check each flag and finally output the number that is not flagged. The whole process will take O(1) time. It has additional O(n) space/memory requirement.
请注意,如上所述,上面提供的两个并行算法可能需要额外的空间。