The Art Gallery Guardian

Arrays and permutation

1 Permutation, functional form

What exactly is an array. In functional form, an array just have to support the following two operations in constant time.

  1. get(A,i)get(A,i): returns A[i]A[i].
  2. set(A,i,x)set(A,i,x): update the array AA such that A[i]A[i] returns xx.

So in some sense, array is just encoding a function f:[n]Xf:[n]\to X. A permutation would be a bijective function π:[n][n]\pi:[n]\to [n]. If we are interested in applying a permutation π\pi, then to program it is easy, we need to output a new function gg such that g(i)=f(π(i))g(i) = f(\pi(i)).

This allows us to apply permutations pretty easily by composing functions and cache outputs. In the purely functional view, the layout of the array in memory and the indexing can be different.

2 Permutation, physically

Sometimes one might ask to physically apply the permutation to an array. That is, the iith position in the array contains the element in π(i)\pi(i). This is helpful because it helps with cache locality: accessing consecutive elements would be in the same location. Although if the ordering of loops does not matter, there is no harm considering the functional view.

Often, one is tasked to apply permutation to an array physically. It usually ask for O(n)O(n) running time and O(1)O(1) space. Unfortunately, there is no way to obtain this running time for all permutations. There are some permutations where this is impossible. A discussion can be found in cstheory.

3 Mix and match

Lingyu Xu asked me about Wiggle Sort, where the actual physical layout has to be changed, but it takes both the functional and physical view of the array.

Let’s consider the simple case where every element is distinct. One simple solution is the following. Partition the numbers into the median, elements smaller than median, and elements larger than the median. We map the smallest n/2n/2 elements into even positions, and remaining elements into the odd position. The algorithm has O(n)O(n) time. The problem is how to get constant extra place. Because it is unclear how to apply the following permutation in place. The permutation σ(i)=(1+2i)%(n1)\sigma(i) = (1+2i) \% (n|1).

However, one we take the functional view, the problem can be solved, physically too. The answer by Stefan Pochmann shows one can simply do the mapping in place.

What is happening in the physical location. If we applied index transform π\pi, then apply physical permutation σ\sigma on the index, what happens for the physical array? It is the same as applying π1σπ\pi^{-1}\sigma\pi to the physical array. Note this means we can only use this trick to obtain conjugates of σ\sigma.

Posted by Chao Xu on .
Tags: Permutation.