Ecosyste.ms: Awesome

An open API service indexing awesome lists of open source software.

Awesome Lists | Featured Topics | Projects

https://github.com/supermitch/chinese-postman

My Python solution to the Chinese-Postman problem.
https://github.com/supermitch/chinese-postman

circuit dijkstra graph network

Last synced: about 6 hours ago
JSON representation

My Python solution to the Chinese-Postman problem.

Awesome Lists containing this project

README

        

# Chinese-Postman Solver

I wrote this program to solve the
[Chinese Postman problem](http://en.wikipedia.org/wiki/Route_inspection_problem).

> The Chinese Postman Problem, or "route inspection problem"
> is to find a shortest closed circuit that visits _every edge_ of a
> (connected) undirected graph.

## Inspiration

I was inspired to learn about and solve this problem when I thought it would
be cool to follow every trail in
[Pacific Spirit Park](http://en.wikipedia.org/wiki/Pacific_Spirit_Regional_Park)
in one run.

Given that the park contains over 73 km of trail, I need to find the optimum
Eularian Path. Otherwise it's going to be a really, really long run!

## The Process

The solution is roughly a three-step process:

1. Determine if the graph has an
[Eularian Path](http://en.wikipedia.org/wiki/Eulerian_path)
(Very easy)
2. Make the non-Eularian graph Eularian, at the minimum expense
(Not so easy)
3. Find the fudged Eularian path
(Pretty easy)

### Solving Minimum Expense

In order to convert a non- or semi-Eularian graph to an Eularian one,
you must eliminate odd nodes (nodes having an odd number of edges.)

To eliminate an odd node, you need to add another edge to it (essentially
retracing your steps.) However, this comes as a cost! The goal then is
to find out which edges to repeat, that eliminate all the odd nodes, with
the minimum cost.

1. Find all possible combinations of pairs of odd nodes
2. Using Dijkstra's Algorithm, find the cost of the minimum path between
those pairs
3. Find which set of paths (depending on how many odd nodes you have)
that results in the least total cost
4. Modify your graph with these new parallel edges

Now you have an Eularian graph with only even nodes, for which an Eularian
Circuit can be found.

### Solving the Eularian Circuit

Solving the Eularian Circuit (now that we have one) is relatively easy. At
first, I simply walked the edges randomly until I happened to find a route
that either dead-ended, or resulted in a circuit. Then I implemented [Fleury's
Algorithm](http://en.wikipedia.org/wiki/Eulerian_path#Fleury.27s_algorithm)
which says always choose a non-bridge over a bridge (for obvious
reasons). Now it takes very few attempts to solve most circuits.

Later I will implement an alternative circuit finding method (Hierholzer's?)

## Usage

This program will probably run in Python 2.7 and definitely in Python 3.4-3.10
(Tested recently w/ Python 3.10.6.)

* Usage: `python main.py -h`
* Specify which graph to load by adding the graph name `python main.py north`
* Optionally specify the starting node, .e.g `python main.py square 4` would produce
`[4, 3, 2, 1, 4]`.

You can find all the graph names in the [data/data.py](data/data.py) file.

## Tests

There are unit tests included, in the `tests` directory.
You can run these via

```bash
python tests/run_tests.py
```

from the root project folder.

## Graph Format

* A graph is defined as a list containing tuples
* Each tuple in the list represents an edge
* Edges are defined as `(start node, end node, length)`

For example, an equilateral triangle like:
```
1
/ \
2 - 3
```
Could be represented as:
`triangle = [(1, 2, 1), (2, 3, 1), (3, 1 ,1)]`

If an edge is [directed](https://en.wikipedia.org/wiki/Directed_graph)
you can add an optional fourth argument, `True`, e.g.
`(1, 2, 5, True)` for a one-way edge from Node 1 to
Node 2 of length 5.

See [network.py](chinesepostman/network.py) for the actual implementation.