Codility峰的O（N * log（log（N）））算法？

Question

任务说明在这里：https://codility.com/demo/take-sample-test/peaks 也在这里：Codility Peaks Complexity

首先，我尝试自己解决这个问题，但只能想出我认为是蛮力解决方案的方法。但是，它的得分为 100/100：https://codility.com/demo/results/demoRNWC3P-X4U

这显然让我完全不满意。 ;) 为N 的每个因子调用外部循环（或为每个峰值，以较小者为准），为每个峰值调用内部循环（只需检查每个块中是否有峰值）。也许那是O(N^2)，也许更好一点（因为它在时间限制内通过了测试）但我几乎可以肯定它不是O(N*log(log(N)))。

然后我尝试搜索O(N*log(log(N))) 解决方案，但其他人似乎都有与我非常相似的解决方案。

那么有人对O(N*log(log(N)))（或更好的）解决方案有想法吗？

另外，如果有人能告诉我我的解决方案有多复杂，我将不胜感激。

Answer 1

您的代码是 O(n d(n))，其中 d(n) 是 n 的除数。在 [1, 100000] 上，d(n) 在 83160 = 2^3 3^3 5 7 11 处最大化，有 126 个除数。 According to Wikipedia，对于每个 epsilon>0，d(n) 为 o(n^epsilon)，因此它的增长相当缓慢。

要获得 O(n log log n) 解决方案，请构建一个部分和数组，告诉您每个点还剩下多少个峰值。然后您可以判断 O(1) 时间间隔内是否存在峰值。然后检查除数 d 需要 O(n/d) 时间。在所有除数 d 上加 n/d 与在所有除数 d 上加 d 相同，根据同一个 Wikipedia 页面，结果是 O(n log log n)。

Answer 2

我已经按照 tmyklebu（谢谢！）建议的风格实现了一个解决方案，它应该是 n.log(log(n))。 Codility 不再测试这个问题的“性能”（！），但 python 解决方案的准确度为 100%。

顺便说一句，如果您一直在学习 Codility 课程，您会记得在 Lesson 8: Prime and composite numbers 中，谐波数运算的总和将产生 O(log(n)) 复杂度。我们有一个简化的集合，因为我们只关注因子分母。 Lesson 9: Sieve of Eratosthenes 展示了素数倒数之和如何为 O(log(log(n))) 并声称“证明是不平凡的”。除数倒数的总和与素数倒数的总和不同，但我建议它也属于“非平凡”证明类别！

def solution(data):

    length = len(data)

    # array ends can't be peaks, len < 3 must return 0    
    if len < 3:
        return 0

    peaks = [0] * length

    # compute a list of 'peaks to the left' in O(n) time
    for index in range(2, length):
        peaks[index] = peaks[index - 1]

        # check if there was a peak to the left, add it to the count
        if data[index - 1] > data[index - 2] and data[index - 1] > data[index]:
            peaks[index] += 1

    # candidate is the block size we're going to test
    for candidate in range(3, length + 1):

        # skip if not a factor
        if length % candidate != 0:
            continue

        # test at each point n / block
        valid = True
        index = candidate
        while index != length:

            # if no peak in this block, break
            if peaks[index] == peaks[index - candidate]:
                valid = False
                break

            index += candidate

        # one additional check since peaks[length] is outside of array    
        if index == length and peaks[index - 1] == peaks[index - candidate]:
            valid = False

        if valid:
            return length / candidate

    return 0