2017-02-27, post № 161
programming, Python PIL, #amazing, #labyrinth, #Minotaur, #solver
Mazes have been a subject of human interest for thousands of years. The Greeks used them to trap a bull-man hybrid, the French built them to show how they could impose order on nature, and even nowadays people enjoy wandering around corn mazes.
The algorithmic art of using computers to solve mazes — and even to find the shortest path through a maze —, however, has only emerged in the last couple of decades.
I was inspired by a recent Computerphile video in which Michael Pound talks about implementing different path finding algorithms for use in maze solving. And as he used Python — one of my favourite languages out there —, I thought I could give it a try and came up with this maze solver.
The mazes given to the solver (through a
.png file) have to have a specific form. The maze needs to have a border all around (painted black) with two holes at the top and bottom, marking the maze’s start and exit (all path pixels are white).
Then the solver — using PIL — reads in the maze file, determines start and exit and starts at the maze’s start, labelling each maze path according to its shortest distance to the start. After it has found the exit, it stops looking at the maze and traces its origins back from the exit, marking the path it goes along as the maze’s optimal solution (highlighted in red).
The different hues of blue indicate the tile’s distance to the start, the white tiles are tiles the solver did not even look at.The different shadings also reveal information about the maze. Mazes with only one solution tend to have sharp changes as there are parts of the maze separated by only one wall, yet separated by a huge walk distance through the maze. The one maze with multiple solutions (see upper right image below) — in contrast — has only gradual changes in hue.
To solve a 𝟦 megapixel maze, the solver takes around 𝟥 seconds, for a 𝟣𝟨 megapixel maze around 𝟣𝟦 seconds and for a 𝟤𝟤𝟧 megapixel maze around 𝟩 minutes and 𝟤𝟤 seconds.
Performance was measured on an Intel Core i7 (𝟦.𝟢𝟢 GHz).
All mazes shown were downloaded from Michael Pound’s mazesolving GitHub repository, which were mostly generated using Daedalus.
The solver’s source code is listed below, though you can also download the
2017-02-25, post № 160
games, Java, programming, #2016, #clone, #exponents, #jar, #powers of two
4096 is a Java-based clone of the well-known web and mobile game 2048, which itself clones 1024 and is similiar to THREES. The naming trend is quite obvious, though note that is a power of two where the exponent is divisible by three, further connecting to the aforementioned game.
In the game, you are faced with a 𝟦 ⨉ 𝟦 matrix, containing powers of two. By swiping in the four cardinal directions (e. g. pressing the arrow keys), you shove all the non-empty cells to that side. When two equal powers of two collide, they fuse together, adding. Once you shoved, an empty tile pseudo-randomly transforms to either a two-tile (𝟫𝟢％) or a four-tile (𝟣𝟢％).
Your objective at first is to reach the tile
4096, though the real goal is to achieve the highest score. Your score is the sum of all the collisions you managed to cause.
To play 4096, you can either download the .jar file or review and compile the game for yourself, using the source code listed below.
- Up, down, left or right arrow key shoves the tiles,
- ‘Escape’ restarts the game upon a loss,
- ‘F11’ toggles fullscreen.
2017-02-11, post № 159
programming, Python, #2016, #puzzle, #shell
Slitherlink is a neat puzzle in which you are presented with a number matrix and the goal is to draw one connected, not interlooping line in between the cells. The number in each cell determines exactly how many line segments must be drawn around each cell (𝟢, 𝟣, 𝟤 or 𝟥). When a cell does not contain any number, the number of line segments adjacent to this cell are unrestricted.
* * * * * * *--* * *--*--*
2 1 3 | | 2 1| 3 |
* * * * * * * *--* *--* *
2 2 2 | | 2 2| 2|
* * * * * * *--* *--* * *
2 -> | 2 | | |
* * * * * * * *--* *--* *
1 3 1 1 1 3| 1 1|
* * * * * * *--*--* *--* *
3 2 3 | 3 2| | 3|
* * * * * * *--*--*--* *--*
A sample 5x5 Slitherlink with solution.
Slitherlink was invented by the Japanese publisher Nikoli in 1989. It has many other names than “Slitherlink”, yet I prefer this descriptive name. Imagining a snake slithering along the board, seeking to link up with itself is a bit charming.
As with most of these puzzles that have simple rules and are fairly easy to work out by hand — on small scales that is —, writing a solver for them can prove to be more difficult than one may expect.
The first solving strategy I tried out was to brute force the problem. Using the Slitherlink from above as an example, there would be different cells with line segments. With each line segment either being drawn or not, there are different boards to check. With one board being checked per nanosecond the solver would take seconds or 𝟥𝟨.𝟧𝟨 years. Brute force is definitely not a viable way to conquer Slitherlink.
After this harsh discovery, I needed a better way to approach solving a given Slitherlink puzzle. Doing some research, I even discovered that Slitherlink is an NP-complete problem (see this paper by Stefan Herting), whereby it — assuming — is not even possible to write a solving algorithm which takes polynomial time.
However, solving small Slitherlink puzzles is fortunately possible in a reasonable time frame.
The strategy I used in the solver consists of pre-programmed rules — figured out by humans — which determine parts of the board based on special arrangements and enforcing the puzzle’s rules (such as that there must only be one line). Using those clues, the solver partly solves a given Slitherlink until there are no more known rules to advance. At that point the solver guesses for a given line segment to be either crossed (marking it cannot be drawn) or drawn, building a tree.
Conflicting attempts (where the solver wrongly guessed, then later — through applying the given rules — determines the attempt as flawed) are thrown away, only leaving possible solved scenarios. Because each Slitherlink has one unique solution, this process ultimately results in one surviving attempt, which then is checked for correctness and printed out as the solution.
A list of Slitherlink rules can be found in this Wikipedia article.
Using the above described method, my solver takes roughly 𝟢.𝟢𝟧 seconds on an Intel Core i7 (𝟦.𝟢𝟢 GHz) to solve the example 𝟧 ⨉ 𝟧 Slitherlink. A 𝟣𝟢 ⨉ 𝟣𝟢 Slitherlink takes around 𝟣.𝟨 seconds whereas it takes 𝟥𝟤 seconds to solve a 𝟣𝟧 ⨉ 𝟣𝟧 Slitherlink. The non-polynomial time is clearly recognisable.
My solver best runs in a bash shell, as it uses ANSI escape sequences to give the solved line a vivid blue and is entirely written in Python as well as fully text-based. The source code is listed below.
Other people also have written solvers, including puzzle generators, such as kakuro-online or appspot. The latter even supports different polygons as the Slitherlink base.
2017-01-28, post № 158
ImageMagick, programming, Python PIL, #animation, #flight, #gif, #light, #physics
Light is a fascinating thing in our universe. We perceive it as color, warmth and vision. Yet it does things one may not expect it to do. One of the experiments that called for a better physical model of light was the double slit experiment. In this experiment, a laser is shone through two closely adjacent slits and projected on the screen behind. Using old physical models, one would expect to see one or maybe two specs of light on the screen, when in reality there appear alternating dark and bright spots.
To explain why this seemingly strange phenomenon is occurring, one can either see light as photons and comprehend that a photon presumably follows every possible path there is in the entire universe and then — through it being observed — randomly chooses one path and thus creates stripes (according to the theory of quantum mechanics) or one can see light as simply being a wave.
For more information on the double-slit experiment, I refer to this Wikipedia entry.
The animation shown below describes light as a wave. The green vectors represent the light wave’s phase at the points on the light beam, the yellow vector represents the addition of both of the slit’s light beam’s phase when hitting the screen and the red dots at the screen represent the light’s brightness at that point (defined by the yellow vector’s length).
To create the animation, Python and a Python module called PIL were used to create single frames which were then stitched together by ImageMagick to create an animated gif.
2017-01-14, post № 157
Java, mathematics, programming, #compiled, #filled Julia set, #generalized filled Julia set, #generalized Mandelbrot set, #Julia, #speed increase, #viewer
Over a year ago, I published my first Mandelbrot Set viewer, a Python program using Pygame. Since then, I have published a rather short program highlighting errors that can occur when calculating the set (Mandelbrot Set Miscalculations).
Since my first viewer was written in Python, which is an interpreted programming language, and I wanted to make my viewer faster, I decided to write one in Java. I was hoping for a speed increase since Java is compiled and thus should run at higher speeds. I was not disappointed. The new Java-based viewer runs noticeably faster and additionally I added a lot of new features, all listed below.
- Left-clicking and dragging draws a zoom frame, single left-clicking removes the frame,
- Right-clicking (and optionally dragging) moves the zoom frame,
- ‘Space’ zooms into the zoom frame,
- ‘F1’ moves one step back the zoom history,
- ‘F2’ shows the path a complex number at the cursor’s position follows when the function is iteratively applied,
- ‘F3’ shows the and axis,
- ‘F4’ displays the current cursor’s corresponding complex number,
- ‘F5’ toggles between showing and hiding the menu (text in the left upper corner describing the viewer’s functions and current states),
- ‘F6’ increments the exponent (going from to in whole-number steps),
- ‘F7’ toggles between the Mandelbrot set and the filled Julia set,
- ‘F8’ toggles between previewing a small filled Julia set at the cursor’s position based upon the cursor’s complex number,
- ‘F9’ completely resets the zoom and zoom history,
- ‘F11’ (or ‘F’) toggles between fullscreen and windowed mode,
- ‘F12’ quits the application,
- ‘L’ increases the color depth (starting at 𝟤𝟧𝟨 and increasing in steps of 𝟤𝟧𝟨),
- ‘Q’ saves the current image to disk.
To use this application, you can either trust me and download the .jar-file or view the source code listed below, verify it and compile the program yourself.
The program will start in fullscreen mode, to change to windowed mode, just press ‘F11’ (as listed above).
2016-12-31, post № 156
art, haiku, poetry, #2017, #Silvester
Aiming at the sky.
Stunningly empty vastness …
Sound of explosion.
2016-12-24, post № 155
art, haiku, poetry, programming, Python, #ascii, #ascii art, #star tree, #tree, #xmas
Christmas tree gets chopped,
Excitement fills the people.
Forming winter mood.
2016-12-18, post № 154
art, haiku, poetry, #nature, #sleep
The lakes are frozen.
All world seems barren and still.
Nature is asleep.
Jonathan Frech's blog; built 2021/10/02 17:36:09 CEST