Find largest cuboid containing only 1's in an NxNxN binary array

Question

Given an NxNxN binary array (containing only 0\'s or 1\'s), how can we obtain the largest cuboid with a non-trivial solution i.e. in O(N^3) ?

--

It is the same p

Answer 1

This solution has O(N³ log² N) complexity (may be optimized to O(N³ log N)). Additional integer array of size 2*8*N³ will be needed.

1. Compute r(i,j,k): for each of the N² rows, compute cumulative sum of all non-zero elements, resetting it when a zero element is found.
2. Perform the following steps for various values of K, using Golden section search (or Fibonacci search) to find the maximum result.
3. Compute c(i,j,k): for each of the N² columns, compute cumulative sum of all elements with r(i,j,k) >= K, resetting it when an element with r(i,j,k) < K is found. For good visualization of steps 1 and 2, see this answer.
4. Perform the last step for various values of M, using Golden section search to find the maximum result.
5. Compute sum: for each of the N² values of 3rd coordinate, compute cumulative sum of all elements with c(i,j,k) >= M, resetting it when an element with c(i,j,k) < M is found. Calculate sumKM and update the best so far result if necessary.

"cross the edge" property of the array is handled in obvious way: iterate every index twice and keep all cumulative sums not larger than N.

For multidimensional case, this algorithm has O(N^D log^D-1 N) time complexity and O(D*N^D) space complexity.

---

Optimization to O(N³ log N)

Step 4 of the algorithm sets a global value for M. This step may be excluded (and complexity decreased by log N) if value for M is determined locally.

To do this, step 5 should be improved. It should maintain a double-ended queue (head of which contains local value of M) and a stack (keeping starting positions for all values of M, evicted from the queue).

While c(i,j,k) increases, it is appended to the tail of the queue.

If c(i,j,k) decreases, all larger values are removed from the queue's tail. If it decreases further (queue is empty), stack is used to restore 'sum' value and put corresponding 'M' value to the queue.

Then several elements may be removed from the head of the queue (and pushed to the stack) if this allows to increase local solution's value.

For multidimensional case, this optimization gives O(N^D log^D-2 N) complexity.