- iOS/Objective-C 元类和类别
- objective-c - -1001 错误,当 NSURLSession 通过 httpproxy 和/etc/hosts
- java - 使用网络类获取 url 地址
- ios - 推送通知中不播放声音
最大子数组和问题有一个非常简单的线性时间解 https://en.m.wikipedia.org/wiki/Maximum_subarray_problem .
如果相反,我们想要最大化 sum(subarray)/sqrt(subarray length),是否存在次二次时间解?
输入数组的元素将是 -infinity 到 +infinity 范围内的浮点值。
最佳答案
我在下面的测试中添加了一个版本的 estabroo 基于 Kadane 的代码。在我的测试中似乎显示出高达 10% 的差异(运行代码段进行随机测试)。
(结束更新)
就近似值而言,我能想到的最好办法是在搜索过程中使用窗口大小的随机样本对目标进行二分搜索 O(log m * n * num_samples_constant)
,其中m
是范围。在测试中,我看到暴力破解(限于 5000 个元素阵列,范围为 ±1000000000)和后者之间的变化在 0% 到 30% 之间变化,样本大小为 200 个窗口长度。 (也许另一个例程可以进一步完善?)
下面的 JavaScript 代码运行 10 次测试并报告最小和最大的差异,然后仅在更长的数组上进行二进制搜索。
其他想法包括使用 FFT 生成总和,但我不知道是否有一种有效的方法可以将每个总和与生成它的子数组长度相关联;以及尝试找到问题的另一种表示形式:
f = sqrt(i - j) * (si - sj), for j < i
f^2 = sqrt(i - j) * (si - sj) * sqrt(i - j) * (si - sj)
= (i - j) * (si^2 - 2si*sj + sj^2)
= i*si^2 - 2i*si*sj + i*sj^2
-j*si^2 + 2j*si*sj - j*sj^2
= i*si^2 +
(-2sj, sj^2, -j, 2j*sj, -j*sj^2) // known before i
dot (i*si, 1, si^2, si, 1)
(因此,如果我们在对数时间内解决了一个 5 维凸包更新,即 5 维极值点问题,并弄清楚我们的候选对象是正数还是负数,我们就可以开始了:)
function prefix_sums(A){
let ps = new Array(A.length + 1).fill(0)
for (let i=0; i<A.length; i++)
ps[i + 1] = A[i] + ps[i]
return ps
}
function brute_force(ps){
let best = -Infinity
let best_idxs = [-1, -1]
for (let i=1; i<ps.length; i++){
for (let j=0; j<i; j++){
let s = (ps[i] - ps[j]) / Math.sqrt(i - j)
if (s > best){
best = s
best_idxs = [j, i - 1]
}
}
}
return [best, best_idxs]
}
function get_high(A){
return A.reduce((acc, x) => x > 0 ? acc + x : acc, 0)
}
function get_low(A){
return Math.min.apply(null, A)
}
function f(A){
let n = A.length
let ps = prefix_sums(A)
let high = get_high(A)
let low = get_low(A)
let best = [-1, -1]
let T = low + (high - low) / 2
let found = false
while (low + EPSILON < high){
T = low + (high - low) / 2
// Search for T
found = false
for (let l=0; l<NUM_SAMPLES; l++){
let w = Math.max(1, ~~(Math.random() * (n + 1)))
for (let i=w; i<ps.length; i++){
let s = (ps[i] - ps[i - w]) / Math.sqrt(w)
if (s >= T){
found = true
best = [i - w, i - 1]
break
}
}
if (found)
break
}
// Binary search
if (found)
low = T
else
high = T - EPSILON
}
return [low, best]
}
function max_subarray(A){
var max_so_far = max_ending_here = A[0]
var startOld = start = end = 0
var divb = divbo = 1
//for i, x in enumerate(A[1:], 1):
for (let i=1; i<A.length; i++){
var x = A[i]
divb = i - start + 1
divbo = divb - 1
if (divb <= 1){
divb = 1
divbo = 1
}
undo = max_ending_here * Math.sqrt(divbo)
max_ending_here = Math.max(x, (undo + x)/Math.sqrt(divb))
if (max_ending_here == x)
start = i
max_so_far = Math.max(max_so_far, max_ending_here)
if (max_ending_here < 0)
start = i + 1
else if (max_ending_here == max_so_far){
startOld = start
end = i
}
}
if (end == A.length-1){
start = startOld + 1
var new_max = max_so_far
divbo = end - startOld + 1
divb = divbo - 1
while (start < end){
new_max = (new_max * Math.sqrt(divbo) - A[start-1])/Math.sqrt(divb)
if (new_max > max_so_far){
max_so_far = new_max
startOld = start
}
start += 1
}
}
return [max_so_far , startOld, end]
}
const EPSILON = 1
const NUM_SAMPLES = 200
let m = 1000000000
let n = 5000
let A
let max_diff = 0
let min_diff = Infinity
let max_diff2 = 0
let min_diff2 = Infinity
let num_tests = 10
for (let i=0; i<num_tests; i++){
A = []
for (let i=0; i<n; i++)
A.push([-1, 1][~~(2 * Math.random())] * Math.random() * m + Math.random())
let f_A = f(A)
let g_A = brute_force(prefix_sums(A))
let m_A = max_subarray(A)
let diff = (g_A[0] - f_A[0]) / g_A[0]
max_diff = Math.max(max_diff, diff)
min_diff = Math.min(min_diff, diff)
let diff2 = (g_A[0] - m_A[0]) / g_A[0]
max_diff2 = Math.max(max_diff2, diff2)
min_diff2 = Math.min(min_diff2, diff2)
}
console.log(`${ n } element array`)
console.log(`${ num_tests } tests`)
console.log(`min_diff: ${ min_diff * 100 }%`)
console.log(`max_diff: ${ max_diff * 100 }%`)
console.log(`min_diff (Kadane): ${ min_diff2 * 100 }%`)
console.log(`max_diff (Kadane): ${ max_diff2 * 100 }%`)
n = 100000
A = []
for (let i=0; i<n; i++)
A.push([-1, 1][~~(2 * Math.random())] * Math.random() * m + Math.random())
var start = +new Date()
console.log(`\n${ n } element array`)
console.log(JSON.stringify(f(A)))
console.log(`${ (new Date() - start) / 1000 } seconds`)
关于algorithm - 最大归一化子数组和的快速算法?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/55781154/
如何更改循环中变量的名称?比如 number1 、 number2 、 number3 、 number4 ? var array = [2,4,6,8] func ap ( number1: Int
我想设置 View 的背景颜色并在一定延迟后将其更改为另一种颜色。这是我的尝试方式: print("setting color 1") self.view.backgroundColor = UICo
我在使用 express-session 时遇到问题。 session 数据不会在请求之间持续存在。 正如您在下面的代码中看到的那样,/join 路由设置了一些 session 属性,但是当 /sur
我试图从叶渲染器获得一个非常简单的结果,用于快速 Steam 的 for 循环。 我正在上传叶文件 HTML,因为它不接受此处格式正确的代码 - 下面的pizza.swift代码- import
你们中有人有什么好的链接可以与我分享吗?我正在寻找一个 FAST 程序员编辑器,它可以非常快速地打开包含超过 100, 000 行代码的文件?我目前正在使用记事本自动取款机,打开一个 29000 行长
我现在正在处理眼动追踪数据,因此拥有一个巨大的数据集(想想数百万行),因此希望有一种快速的方法来完成此任务。这是它的简化版本。 数据告诉您眼睛在每个时间点正在查看的位置以及我们正在查看的每个文件。 X
我是新手,想为计时器或其他设备选择提示音。 如何打开此列表,以选择其中一种声音? Alert sound list 最佳答案 您将无法在应用中使用系统声音。 但是,您可以包括自己的声音文件,并将其显示
我编写了以下代码来构建具有顺序字符串的数组。 它的工作方式与我预期的一样,但我希望它能更快地运行。有没有更有效的方法在PowerShell中产生我想要的结果? 我是PowerShell的新手,非常感谢
我有一个包含一些非唯一行的矩阵,例如: x 尝试 y <- rle(apply(x, 1, paste, collapse = " ")) # y$lengths is the vector con
我的函数“keyboardWillShown”有问题。所以我想要的是菜单打开时,菜单正好出现在键盘上方。它可以在Iphone 8 plus,8、7、6上完美运行。但是,当我在模拟器上运行Iphone
我正在尝试通过Swift 5中的HTTP get方法从API提取数据。它在启动时成功加载了数据,但是当我刷新页面时,它说“索引超出范围”,这是因为数据是不再会在我的日志中读取,因此索引中没有任何内容。
我想做什么: 从我的数据库中获取时间戳并将其转换为用户的时区。 我的代码: let tryItNow = "\(model.timestampName)" let format = D
给定字体名称和字体大小,如何查找字符串的宽度(CGFloat)? (目标是将UIView的宽度设置为足以容纳字符串的宽度。) 我有两个字符串:一个重复“1”,重复36次,另一个重复“M”,重复36次。
我正在尝试解析此JSON ["Items": ( { AccountBalance = 0; AlphabetType = 3; Description = "\U0631\U
我在UINavigationBar内放置了一个UILabel。 我想根据navigationBar的高度增加该标签的字体大小。当navigationBar很大时,我希望字体大小更大;当滚动并缩小nav
我想将用户输入限制为仅有效数字并使用以下内容: func textView(_ textView: UITextView, shouldChangeTextIn range: NSRange, rep
目前我有一个包含超过 100.000 张图像的数据库,它们大小不一或类似,但我想为我的公司制作以下内容: 我插入/上传一张图片,系统返回最有可能相同的图片。我不知道使用什么算法,但它需要快速。我可以预
在我的 swift 项目中,我有一个按钮,我想在标签上打印按下该按钮的时间。 如何解决这个问题? 最佳答案 添加到DHEERAJ的答案中,您只需在func press(sender: UIButton
我必须发表评论,尝试在解析中导入数组。然而,有一个问题。 当我尝试从 Parse 加载数组时,我的输出是 ("Blah","Blah","Blah")这是一个元组...而不是一个数组 TT... 如何
我的应用程序有一个名为 MyDevice 的类,我用它来与硬件通信。该硬件是可选的,实例变量也是可选的: var theDevice:MyDevice = nil 然后,在应用程序中,我必须初始化设备
我是一名优秀的程序员,十分优秀!