Tuesday 11 March 2014

Backtracking Algorithms

Solution Spaces

Backtracking is a refinement of the brute force approach, which systematically searches for a solution to a problem among all available options. It does so by assuming that the solutions are represented by vectors (v1, ..., vm) of values and by traversing, in a depth first manner, the domains of the vectors until the solutions are found.
When invoked, the algorithm starts with an empty vector. At each stage it extends the partial vector with a new value. Upon reaching a partial vector (v1, ..., vi) which can’t represent a partial solution, the algorithm backtracks by removing the trailing value from the vector, and then proceeds by trying to extend the vector with alternative values.
ALGORITHM try(v1,...,vi) 
   IF (v1,...,vi) is a solution THEN RETURN (v1,...,vi) 
   FOR each v DO 
      IF (v1,...,vi,v) is acceptable vector  THEN 
        sol = try(v1,...,vi,v) 
        IF sol != () THEN RETURN sol 
      END 
   END 
   RETURN () 
If Si is the domain of vi, then S1 × ... × Sm is the solution space of the problem. The validity criteria used in checking for acceptable vectors determines what portion of that space needs to be searched, and so it also determines the resources required by the algorithm.
The traversal of the solution space can be represented by a depth-first traversal of a tree. The tree itself is rarely entirely stored by the algorithm in discourse; instead just a path toward a root is stored, to enable the backtracking.
                     •
            |------------------|
v1 :     |-•||--------S-------•|-|-
      |-----|-----------1-----|------|
v2 :|•----|•|S||--•|-  |•----|•|S||--•|-
    |----|||||2-----|  |----|||||2-----|
   ... ...... ... ... ... ... ...... ... ... ...

Traveling Salesperson

The problem assumes a set of cities, and a salesperson which needs to visit each city exactly once and return to the base city at the end. The solution should provide a route of minimal length.
The route (a, b, d, c) is the shortest one for the following one, and its length is 51.
do-----12-|co
 | ||| ||| |
 |   |||   |
15|20   9|18
ao|  -----|o
   10      b
The traveling salesperson problem is an NP-hard problem, and so no polynomial time algorithm is available for it. Given an instance = (V, E) the backtracking algorithm may search for a vector of cities (v1, ..., v||) which represents the best route.
The validity criteria may just check for number of cities in of the routes, pruning out routes longer than ||. In such a case, the algorithm needs to investigate |||| vectors from the solution space.
bo----
 |   ----
 |    ---ao
 |-----
co-
                         •
        |---------------------------------|
  |----a------|    |----b------|    |-----c-----|

|-a-| |b--| |c-| |-a-| |b--| |c--||-a-| |-b-| |c--|

a b ca b c a b c a b c a b ca b c a b c a b ca b c
On the other hand, the validity criteria may check for repetition of cities, in which case the number of vectors reduces to ||!.
    •|
 |-------|
|a| b-| c-|

b c ac|a|b|

c b ca b a

The Queens Problem

Consider a by chess board, and the problem of placing queens on the board without the queens threatening one another.
|-|-o|-|-o|-|-o|-|--|
|-|--|-|--|-|--|-|--|
|-|--|o|-o|o|--|-|--|
|o|-o|o|-o|o|-o|o|o-|
|-|--|o|-o|o|--|-|--|
|-|-o|-|-o|-|-o|-|--|
|o|--|-|-o|-|--|o|--|
| |  | | o| |  | |o
|-|--|-|-o|-|--|-|--|
--------------------
The solution space is {123, ..., n}n. The backtracking algorithm may record the columns where the different queens are positioned. Trying all vectors (p1, ..., pn) implies nn cases. Noticing that all the queens must reside in different columns reduces the number of cases to n!.
For the latter case, the root of the traversal tree has degree n, the children have degree 1, the grand children degree 2, and so forth.
             ||-||
             ||-||
     |-------------------|
    1          2         3
   •|-||     ||•||     ||-•|
   ||-||     ||-||     ||-||
  |----|    |----|    |----|
|1,2||1,3||2,1||2,3||3,1|-3,2-|
|••|||•|•||••|||-••||•|•|-|••-|
--|----|----|-- -|----|----|---
1,2,3 1,3,2  2,1,3  2,3,1 3,1,2 3,2,1
|•||||•||||-•|||-•|||-|•|-||•-|
|-••||-••||•|•||•|•||••||-••|-|
--------------- ---------------
Checking for threatening positions along the way my further reduce the number of visited configurations.
             ||-||
             ||-||
     |-------------------|
    1          2         3
   •|-||     ||•||     ||-•|
   ||-||     ||-||     ||-||
  |----|    |----|    |----|
|1,2||1,3||2,1||2,3||3,1|-3,2-|
|••|||•|•||••|||-••||•|•|-|••-|
-------|------- ------|--------
     1,3,2            3,1,2
     |•|||          |-|•|-
     |-••|          |••||-
     -----           -----

Convex Hull (Graham’s Scan)

The problem asks to construct the shortest polygon which encloses a given set of points on a plan. Intuitively, the polygon can be determined by placing a rubber around the set.
 •
  |---
 |   ----
 |•     ---
 |         ---
•|-----•     ----
       ----------•
Determine an extreme point with the largest coordinate. Sort the points in order of increasing angles, relatively to the extreme point.
  •
b  ---
     ----
   •--  ----
  c   --------
 •----d-•--------
e       ---------•
                a
Traverse the points in the sorted order, adding them to the partial solution upon making a turn of less than 180 degrees, and backtracking when making a larger turn.
  b
 •|-
  |----
  |c  ----
  •      ----
       •d   ---
•e            ---a
                •         b
 •|--
  ||||-
  |c -----
  •----  ----
 e    -•d   ----
•              -•a         b
 •|--
  | ----
  •c  |||-
    --| d|---
 e-----•  - ----
•    ||||-     -•a         b
 •|--
  | ----
  •c-  ----
    --- d ---
 e     •     ---
•              -•a        •b
  |---
  |  ---
  •c   ----
        d ----
•e     •     ----
                •a        •b
  |---
  ||||---
  •c -  ---
 ||||-  d  ---
•e     •     ----
                •a        •b-
  |---
  |  ----
  •c    ----
       •d  ---
•e            ---a
                •       •b-
   ----
   c  ----
  •      ---
       •d  ----
•e            ---a
                •         b
 •--
 ||-|--
 | c  ----
 |•      ----
 e     •d   ---
•             ---a
                •         b
 •---
 |  ---
 | c   ---
 |•     d----
 e||   •    ----
•---------------•a
The algorithm requires O(n log n) time for sorting the points, and O(n) time to select an appropriate subset.

Generating Permutations

A permutation can be obtained by selecting an element in the given set and recursively permuting the remaining elements.
             {
              ai,P(a1,...,ai-1,ai+1,...,aN)  if N > 1
P(a1,...,aN) =  aN                         if N = 1
                    --|--|--|-|
                    |a|b-|c-d-|
    a|------------b------------c-------------d
--|--|--|-|  ---|-|--|--| ---|--|-|--|  --|--|--|-|
|-|b-|c-d-|  |a-|-|c-|d-| |a-|b-|-|d-|  |a|b-|c-|-|
At each stage of the permutation process, the given set of elements consists of two parts: a subset of values that already have been processed, and a subset that still needs to be processed. This logical seperation can be physically realized by exchanging, in the i’th step, the i’th value with the value being chosen at that stage. That approaches leaves the first subset in the first i locations of the outcome.
                    --|--|--|-|
                    |a|b-|c-d-|
--|--|------------|--------------------------|--|-|
a||b |c d |  |b a |c |d | |c||b |a|d |  |d|b |c |a|
-----|--------------------------|----   -----------
--|--|--|-|  ---|-|--|--| ---|--|-|--|
b-|a-|c-d-|  |b-|c|a-|d-| |b-|d-|c|a-|
      ---|--|------------|--|-|
      |b-|c-a-|d-|  b-|c-|d-|a|
                         |
                    b-|c-|d-|a|
                    |-|--|--|-|
permute(i) 
   if i == N  output A[N] 
   else 
      for j = i to N do 
         swap(A[i], A[j]) 
         permute(i+1) 
         swap(A[i], A[j]) 

No comments:

Post a Comment

 

FACEBOOK PAGE

SKETCHES & PAINTINGS