Math.Sym

Provides an efficient definition of standard permutations, StPerm, together with an abstract class, Perm, whose functionality is largely inherited from StPerm using a group action and the standardization map.

Standard permutations

data StPerm

By a standard permutation we shall mean a permutations of [0..k-1].

Instances

Eq StPerm

Ord StPerm

Show StPerm

Monoid StPerm

Perm StPerm

toList :: StPerm -> [Int]

Convert a standard permutation to a list.

fromList :: [Int] -> StPerm

Convert a list to a standard permutation. The list should a permutation of the elements [0..k-1] for some positive k. No checks for this are done.

sym :: Int -> [StPerm]

The list of standard permutations of the given size (the symmetric group). E.g.,

 sym 2 == [fromList [0,1], fromList [1,0]]

The permutation typeclass

class Perm a where

The class of permutations. Minimal complete definition: st act and idperm. The default implementations of size and neutralize can be somewhat slow, so you may want to implement them as well.

Methods

st :: a -> StPerm

The standardization map. If there is an underlying linear order on a then st is determined by the unique order preserving map from [0..] to that order. In any case, the standardization map should be equivariant with respect to the group action defined below; i.e., it should hold that

 st (u `act` v) == u `act` st v

act :: StPerm -> a -> a

A (left) group action of StPerm on a. As for any group action it should hold that

 (u `act` v) `act` w == u `act` (v `act` w)   &&   neutralize u `act` v == v

size :: a -> Int

The size of a permutation. The default implementation derived from

 size == size . st

This is not a circular definition as size on StPerm is implemented independently. If the implementation of st is slow, then it can be worth while to override the standard definiton; any implementation should, however, satisfy the identity above.

idperm :: Int -> a

The identity permutation of the given size.

inverse :: a -> a

The group theoretical inverse. It should hold that

 inverse == unst . inverse . st

and this is the default implementation.

ordiso :: StPerm -> a -> Bool

Predicate determining if two permutations are order-isomorphic. The default implementation uses

 u `ordiso` v  ==  u == st v

Equivalently, one could use

 u `ordiso` v  ==  inverse u `act` v == idperm (size u)

unstn :: Int -> StPerm -> a

The inverse of the standardization function. For efficiency reasons we make the size of the permutation an argument to this function. It should hold that

 unst n w == w `act` idperm n

and this is the default implementation. An un-standardization function without the size argument is given by unst below.

unst :: Perm a => StPerm -> a

The inverse of st. It should hold that

 unst w == unstn (size w) w

and this is the default implementation.

Instances

Perm String

Perm StPerm

Perm [Int]

Convenience functions

empty :: Perm a => a

The empty permutation.

one :: Perm a => a

The one letter permutation.

toVector :: Perm a => a -> Vector Int

Convert a permutation to a vector.

fromVector :: Perm a => Vector Int -> a

Convert a vector to a permutation. The vector should be a permutation of the elements [0..k-1] for some positive k. No checks for this are done.

bijection :: StPerm -> Int -> Int

The bijective function defined by a permutation.

generalize :: (Perm a, Perm b) => (StPerm -> StPerm) -> a -> b

Generalize a function on StPerm to a function on any permutations:

 generalize f = unst . f . st

generalize2 :: (Perm a, Perm b, Perm c) => (StPerm -> StPerm -> StPerm) -> a -> b -> c

Like generalize but for functions of two variables

normalize :: (Ord a, Perm a) => [a] -> [a]

Sort a list of permutations with respect to the standardization and remove duplicates

cast :: (Perm a, Perm b) => a -> b

Cast a permutation of one type to another

Constructions

(\+\) :: Perm a => a -> a -> a

The direct sum of two permutations.

dsum :: Perm a => [a] -> a

The direct sum of a list of permutations.

(/-/) :: Perm a => a -> a -> a

The skew sum of two permutations.

ssum :: Perm a => [a] -> a

The skew sum of a list of permutations.

inflate :: Perm a => a -> [a] -> a

inflate w vs is the inflation of w by vs. It is the permutation of length sum (map size vs) obtained by replacing each entry w!i by an interval that is order isomorphic to vs!i in such a way that the intervals are order isomorphic to w. In particular,

 u \+\ v == inflate (fromList [0,1]) [u,v]
 u /-/ v == inflate (fromList [1,0]) [u,v]

Generating permutations

unrankPerm :: Perm a => Int -> Integer -> a

unrankPerm u rank is the rank-th (Myrvold & Ruskey) permutation of size n. E.g.,

 unrankPerm 9 88888 == "561297843"

randomPerm :: Perm a => Int -> IO a

randomPerm n is a random permutation of size n.

perms :: Perm a => Int -> [a]

All permutations of a given size. E.g.,

 perms 3 == ["123","213","321","132","231","312"]

Sorting operators

stackSort :: Perm a => a -> a

One pass of stack-sort.

bubbleSort :: Perm a => a -> a

One pass of bubble-sort.

Permutation patterns

copiesOf :: Perm a => StPerm -> a -> [Set]

copiesOf p w is the list of (indices of) copies of the pattern p in the permutation w. E.g.,

 copiesOf (st "21") "2431" == [fromList [1,2],fromList [0,3],fromList [1,3],fromList [2,3]]

avoids :: Perm a => a -> [StPerm] -> Bool

avoids w ps is a predicate determining if w avoids the patterns ps.

avoiders :: Perm a => [StPerm] -> [a] -> [a]

avoiders ps vs is the list of permutations in vs avoiding the patterns ps. This is equivalent to the definition

 avoiders ps = filter (`avoids` ps)

but is usually much faster.

av :: [StPerm] -> Int -> [StPerm]

av ps n is the list of permutations of [0..n-1] avoiding the patterns ps. E.g.,

 map (length . av [st "132", st "321"]) [1..8] == [1,2,4,7,11,16,22,29]

Single point extensions/deletions, shadows and downsets

del :: Perm a => Int -> a -> a

Delete the element at a given position

shadow :: (Ord a, Perm a) => [a] -> [a]

The list of all single point deletions

downset :: (Ord a, Perm a) => [a] -> [a]

The list of permutations that are contained in at least one of the given permutaions

ext :: Perm a => Int -> a -> a

Extend a permutation by inserting a new largest element at the given position

coshadow :: (Ord a, Perm a) => [a] -> [a]

The list of all single point extensions

Left-to-right maxima and similar functions

lMaxima :: Perm a => a -> Set

The set of indices of left-to-right maxima.

lMinima :: Perm a => a -> Set

The set of indices of left-to-right minima.

rMaxima :: Perm a => a -> Set

The set of indices of right-to-left maxima.

rMinima :: Perm a => a -> Set

The set of indices of right-to-left minima.

Components and skew components

components :: Perm a => a -> Set

The set of indices of components.

skewComponents :: Perm a => a -> Set

The set of indices of skew components.

Simple permutations

simple :: Perm a => a -> Bool

A predicate determining if a given permutation is simple.

Subsets

type Set = Vector Int

A set is represented by an increasing vector of non-negative integers.

subsets :: Int -> Int -> [Set]

subsets n k is the list of subsets of [0..n-1] with k elements.