jfrech.com
jblog
table of contents

TImg

2017-10-07, post № 181

BASIC, PIL, programming, Python, TI-84 Plus, #bitmap, #image, #TI, #TI-BASIC

Texas Instrument’s TI-84 Plus is a graphing calculator with a variety of features. It has built-in support for both fractions and complex numbers, can differentiate and integrate given functions and supports programming capabilities. The latter allows to directly manipulate the calculator’s monochrome display’s 𝟧𝟫𝟪𝟧 pixels (the screen has dimensions 𝟫𝟧 ⨉ 𝟨𝟥). TImg is a Python program (source code is listed below and can also be downloaded) which takes in an image and outputs TI-BASIC source code which, when run on the graphing calculator, will produce the given image — in potentially lower quality.

timg_dim_photo.jpg
TI-84 Plus’ screen dimensions (bitmap).

PIL — the Python Imaging Library — is used to read in the image and further for processing. The supplied image may be rotated and resized to better fit the TI-84’s screen and any color or even grayscale information is reduced to an actual bitmap — every pixel only has two distinct values.
Direct pixel manipulation on the TI-84 is done via the Graph screen. To get remove any pixels the system draws on its own, the first three commands are ClrDraw, GridOff and AxesOff which should result in a completely blank screen — assuming that no functions are currently being drawn. All subsequent commands are in charge of drawing the previously computed bitmap. To turn certain pixels on, Pxl-On(Y,X is used where 𝑌 and 𝑋 are the pixel’s coordinates.

timg_fractal_photo.jpg
A fractal (bitmap).

Since the TI-84 Plus only has 𝟤𝟦 kilobytes of available RAM, the source code for a program which would turn on every single pixel individually does not fit. Luckily, though, a program which only individually turns on half of the screen’s pixels fits. To ensure that TImg’s output fits on the hardware it is designed to be used with, an image’s bitmap is inverted when the required code would otherwise exceed 𝟥𝟧𝟢𝟢 lines — a value slightly above the required code to draw half of the pixels.

timg_tiception_photo.jpg
A fractal (input image).
Source code: timg.py

Rainbowify

2017-09-23, post № 180

PIL, programming, Python, #color model, #hsl, #hsv, #image filter, #images, #nature, #photography, #rainbow

To digitally represent colors, one most often uses the RGB color system. By combining three fundamental light colors in certain ways, one can define a variety of different wavelengths of light. The human eye has three distinct photoreceptors for the aforementioned three colors, nearly all screens use pixels consisting of three parts in those colors and most image formats store the image data in the RGB color system.

rainbowify-137_rainbowified.jpg
Honey bee (original)

However, there are other color systems than RGB with other strengths. Cycling through the colors of the rainbow, for example, is a lot easier using the HSL (or HSV) color model, as it is simply controlled by the hue.

rainbowify-154_rainbowified.jpg
Fruit (original)

Rainbowify uses the HSL color model to rainbowify a given image. To do so, the image is first converted into a grayscale image (averaging all three color channels). A pixel’s brightness is then interpreted as its hue with its saturation and lightness set to the maximum. As a final touch, the hue gets offset by a pixel-position dependent amount to create the overall appearance of a rainbow.
Source code is listed below and can also be downloaded.

rainbowify-124_rainbowified.jpg
Sunflower (original)
rainbowify-123_rainbowified.jpg
Thistle (original)
Source code: rainbowify.py
Extra assets: rainbowify-120.jpg, rainbowify-120_rainbowified.jpg, rainbowify-121.jpg, rainbowify-121_rainbowified.jpg, rainbowify-122.jpg, rainbowify-122_rainbowified.jpg, rainbowify-125.jpg, rainbowify-125_rainbowified.jpg, rainbowify-126.jpg, rainbowify-126_rainbowified.jpg, rainbowify-127.jpg, rainbowify-127_rainbowified.jpg, rainbowify-128.jpg, rainbowify-128_rainbowified.jpg, rainbowify-129.jpg, rainbowify-129_rainbowified.jpg, rainbowify-130.jpg, rainbowify-130_rainbowified.jpg, rainbowify-131.jpg, rainbowify-131_rainbowified.jpg, rainbowify-132.jpg, rainbowify-132_rainbowified.jpg, rainbowify-133.jpg, rainbowify-133_rainbowified.jpg, rainbowify-134.jpg, rainbowify-134_rainbowified.jpg, rainbowify-135.jpg, rainbowify-135_rainbowified.jpg, rainbowify-136.jpg, rainbowify-136_rainbowified.jpg, rainbowify-138.jpg, rainbowify-138_rainbowified.jpg, rainbowify-139.jpg, rainbowify-139_rainbowified.jpg, rainbowify-140.jpg, rainbowify-140_rainbowified.jpg, rainbowify-141.jpg, rainbowify-141_rainbowified.jpg, rainbowify-142.jpg, rainbowify-142_rainbowified.jpg, rainbowify-143.jpg, rainbowify-143_rainbowified.jpg, rainbowify-144.jpg, rainbowify-144_rainbowified.jpg, rainbowify-145.jpg, rainbowify-145_rainbowified.jpg, rainbowify-146.jpg, rainbowify-146_rainbowified.jpg, rainbowify-147.jpg, rainbowify-147_rainbowified.jpg, rainbowify-148.jpg, rainbowify-148_rainbowified.jpg, rainbowify-149.jpg, rainbowify-149_rainbowified.jpg, rainbowify-150.jpg, rainbowify-150_rainbowified.jpg, rainbowify-151.jpg, rainbowify-151_rainbowified.jpg, rainbowify-152.jpg, rainbowify-152_rainbowified.jpg, rainbowify-153.jpg, rainbowify-153_rainbowified.jpg, rainbowify-155.jpg, rainbowify-155_rainbowified.jpg, rainbowify-156.jpg, rainbowify-156_rainbowified.jpg, rainbowify-157.jpg, rainbowify-157_rainbowified.jpg, rainbowify-158.jpg, rainbowify-158_rainbowified.jpg, rainbowify-159.jpg, rainbowify-159_rainbowified.jpg

Arithmetic Golfing

2017-09-09, post № 179

code golf, mathematics, programming, Python, #byte count, #code golfing, #golfing, #short code

A recent PCG golfing question When do I get my sandwich? asked to find a mapping between seven input strings (sandwich names) and the seven days of the week (indexed by number).

The first answer was made by a user named i cri everytim and utilized a string of characters which uniquely appear at the same position in all seven input strings, enklact, to perform the mapping in Python 2 requiring 𝟤𝟫 bytes. After their answer, a lot of answers appeared using the same magic string in different languages to reduce the number of bytes needed. Yet nobody reduced the byte count in Python.

Trying to solve the problem on my own, my first attempt was using only the input strings’ last decimal digit to perform the mapping, though this approach did not save on bytes (read my PCG answer for more on this 𝟥𝟢 byte solution).

After a few more hours of working on this problem, however, I achieved to bring down the byte count by one entire byte.

I did so by using a simple brute-force algorithm to check for Python expressions which can be used to perform the sought after mapping. To do so, I use Python’s apostrophes (`...`) to turn the found expression into a string — str(...) is three whole bytes longer — and index that string with the input strings’ lengths. It sure is not very readable, but only takes 𝟤𝟪 bytes — and that is all that matters.

lambda S:`6793**164`[len(S)]

After finding the 𝟤𝟪 byte function which uses a 𝟫 byte expression (6793**164), I attempted to find an even shorter expression. And even though I did not yet find one, I did write a more general brute-force Python program (source code shown below; can also be downloaded) than the one I linked to in my PCG answer.

Brute-forcing takes exponentially more time the more digits you have to check, so my brute-forcer still requires the user to decide for themselves which expressions should be tried.
There are three parameters that define the search; a regex pattern that should be contained in the expression’s string, an offset that pattern should ideally have and a target length. If an expression is found that takes as many bytes as or less bytes than the target length, an exclamation point is printed
Though this program did not prove useful in this case, there may come another challenge where an arithmetic expression golfer could come in handy.

My program may not have found shorter expressions, but definitely some impressive ones (the +... at the end refers to an additional offset from the string index which — unsurprisingly — take additional bytes):

I also considered using division to generate long strings of digits which may match; the only problem is that Python floating-point numbers only have a certain precision which does not produce long enough strings. Again, using exponentiation (**) and bitshifting (<<) I could not come up with a working expression that takes less bytes.

Source code: arithmetic-golfing_brute.py

brainfuck X

2017-08-26, post № 178

brainfuck, programming, Python, #ANSI, #argparse, #braindraw, #color, #esoteric, #interpreter, #PPM, #tape, #Turing machine

While browsing StackExchange PCG [1] questions and answers, I came across a challenge regarding drawing the swiss flag. In particular, I was interested in benzene’s answer, in which they showcased a brainfuck dialect capable of creating two-dimensional 𝟤𝟦-bit color images. In this post I present this dialect with slight changes of my own, as well as an interpreter I wrote in Python 2.7 (source code is listed below and can also be downloaded).

brainfuck-x_swiss.png

Urban Müller’s original brainfuck (my vanilla brainfuck post can be found here) works similar to a Turing machine, in that the memory consists of a theoretically infinitely large tape with individual cells which can be modified. What allows brainfuck X (or braindraw, as benzene called their dialect) to create color images is, that instead of a one-dimensional tape, a three-dimensional tape is used. This tape extends infinitely in two spacial dimensions and has three color planes. Each cell’s value is limited to a byte (an integer value from 𝟢 to 𝟤𝟧𝟧) which results in a 𝟤𝟦-bit color depth.

Adding to brainfuck’s eight commands (+-<>[].,), there are two characters to move up and down the tape (^v) and one character to move forwards in the color dimension (*). Starting on the red color plane, continuing with the green and ending in the blue. After the blue color plane, the color planes cycle and the red color plane is selected. benzene’s original language design which I altered slightly had three characters (rgb) to directly select a color plane. Whilst this version is supported by my interpreter )the flag --colorletters is necessary for that functionality(, I find my color star more brainfucky — directly calling color planes by their name seems nearly readable.
brainfuck’s vanilla eight characters still work in the same way, brainfuck X can thereby execute any vanilla brainfuck program [2]. Also, there still is a plaintext output — the tape’s image is a program’s secondary output.

Having executed the final brainfuck instruction, the interpreter prints out the tape to the terminal — using ANSI escape codes. Because of this, the color depth is truncated in the terminal view, as there are only 𝟤𝟣𝟨 colors supported. [3]
For the full 𝟤𝟦-bit color depth output, I use the highly inefficient Portable Pixmap Format (.ppm) as an output image file format. To open .ppm files, I recommend using the GNU Image Manipulation Program; specifying the output file name is done via the --output flag.

The Swiss flag image above was generated by benzene’s braindraw code (see their StackExchange answer linked to above); the resulting .ppm file was then scaled and converted using GIMP.
Interpreter command: python brainfuckx.py swiss.bfx -l -o swiss.ppm

Usage

  • Being written in pure Python, the interpreter is completely controlled via the command line. The basic usage is python brainfuck-x.py <source code file>; by using certain flags the functionality can be altered.
  • --input <input string>, -i <input string> specifies brainfuck’s input and is given as a byte stream (string).
  • --simplify, -s outputs the source code’s simplified version; the source code with all unnecessary characters removed.
  • --colorstar selects the color star color plane change model which is the default.
  • --colorletters, -l selects the color letter color plane change model.
  • --silent stops the interpreter from outputting warnings, infos and the final tape.
  • --maxcycles <cycles>, -m <cycles> defines the maximum number of cycles the brainfuck program can run; the default is one million.
  • --watch, -w allows the user to watch the program’s execution.
  • --watchdelay <delay> defines the time in seconds the interpreter sleeps between each watch frame.
  • --watchskip <N> tells the interpreter to only show every 𝑁th cycle of the execution.
  • --output <output file name>, -o <output file name> saves the final tape as a .ppm image file.
Source code: brainfuck-x.py

Cyclic Quine

2017-08-12, post № 177

Programming, Python, Quine, #self-referential

A classic quine is a program which outputs its own source code.
At first, such a program’s existence seems weird if not impossible, as it has to be so self-referential that it knows about itself everything, including how to know about itself. However, writing quines is possible, if not [1] trivial.

A cyclic quine therefore is a program which outputs source code which differs from its own source code, yet outputs the original source code when run (the cycle length could be greater than one). So when running source codes 𝛹 and 𝛷, they output source codes 𝛷 and 𝛹.

Therefore, when one saves the first program as q0.py and the second as q1.py, one can create both source codes from one another (the following bash commands [2] will not change the files’ contents).

$ python q0.py > q1.py
$ python q1.py > q0.py
$ python q0.py | python > q0.py
$ python q0.py | python | python | python | python > q1.py
Source code: cyclic-quine_q0.py
Source code: cyclic-quine_q1.py

Mandelbrot Set III

2017-07-29, post № 176

Java, mathematics, programming, #dragging, #fractal, #fractals, #Mandelweed, #mouse control, #multithreading, #scrolling, #thread, #viewer

I wrote my first ever Mandelbrot Set renderer back in 2015 and used Python to slowly create fractal images. Over a year later, I revisited the project with a Java version which — due to its code being actually compiled — ran much faster, yet had the same clunky interface; a rectangle the user had to draw and a key they had to press to the view change to the selected region.
In this post, over half a year later, I present my newest Mandelbrot Set fractal renderer (download the .jar), written in Java, which both runs fast and allows a much more intuitive and immersive walk through the complex plane by utilizing mouse dragging and scrolling.
The still time demanding task of rendering fractals — even in compiled languages — is split up into a low quality preview rendering, a normal quality display rendering and a high quality 4K (UHD-1 at 𝟥𝟪𝟦𝟢 ⨉ 𝟤𝟣𝟨𝟢 pixels to keep a 𝟣𝟨 : 𝟫 image ratio) rendering, all running in seperate threads.

mandelbrot-set-iii_3840x2160_Zre-0.5205350147993641_Zim0.5100929168733932_Zom0.0012332949501170922_Clr256.png
Rainbow spiral

The color schemes where also updated, apart from the usual black-and-white look there are multiple rainbow color schemes which rely on the HSB color space, zebra color schemes which use the iterations taken modulo some constant to define the color and a prime color scheme which tests if the number of iterations taken is prime.

mandelbrot-set-iii_3840x2160_Zre-0.7452258168703411_Zim0.09808819723426355_Zom0.001532723995106021_Clr256.png
Zebra spiral

Apart from the mouse and keyboard control, there is also a menu bar (implemented using Java’s JMenuBar) which allows for more conventional user input through a proper GUI.

Controls

  • Left mouse dragging: pan view,
  • Left mouse double click: set cursor’s complex number to image center,
  • Mouse scrolling: zoom view,
  • Mouse scrolling +CTRL: pan view,
  • ‘p’: render high definition fractal,
  • ‘r’: reset view to default,
  • ‘w’, ‘s’: zoom frame,
  • Arrow keys: pan view,
  • Arrow keys +CTRL: zoom view,
  • Menu bar

    • “Fractal”: extra info about current fractal rendering,
    • “Color Scheme”: change color scheme and maximum iteration depth,
    • “HD”: controls for high definition rendering,
    • “Extra”: help and about.
mandelbrot-set-iii_3840x2160_Zre-0.5800917751448121_Zim0.4886972982384146_Zom0.002865014852390483_Clr256.png
Blue spiral

A bit more on how the three threads are implemented.
Whenever the user changes the current view, the main program thread renders a low quality preview and immediately draws it to the screen. In the background, the normal quality thread (its pixel dimensions match the frame’s pixel dimensions) is told to start working. Once this medium quality rendering is finished, it is preferred to the low quality rendering and gets drawn on the screen.
If the user likes a particular frame, they can initiate a high quality rendering (4K UHD-1, 𝟥𝟪𝟦𝟢 ⨉ 𝟤𝟣𝟨𝟢 pixels( either by pressing ‘q’ or selecting “HD ❯ Render current frame”. This high quality rendering obviously takes some time and a lot of processing power, so this thread is throttled by default to allow the user to further explore the fractal. Throttling can be disabled through the menu option “HD ❯ Fast rendering”. There is also the option to tell the program to exit upon having finished the last queued high definition rendering (“HD ❯ Quit when done”).
The high definition renderings are saved as .png files and named with their four defining constants. Zim and Zre define the image’s complex center, Zom defines the complex length above the image’s center. Clr defines the number of maximum iterations.

mandelbrot-set-iii_3840x2160_Zre-0.1512120462044041_Zim0.6496162871560206_Zom1.8511095575145676E-4_Clr256.png
Another blue spiral

Just to illustrate how resource intensive fractal rendering really is.
A 4K fractal at 𝟥𝟪𝟦𝟢 ⨉ 𝟤𝟣𝟨𝟢 pixels with an iteration depth of 𝟤𝟧𝟨 would in the worst case scenario (no complex numbers actually escape) require 3840\cdot 2160\cdot 256\cdot 4=8493465600 double multiplications. If you had a super-optimized CPU which could do one double multiplication a clock tick (which current CPUs definitely cannot) and ran at 𝟦.𝟢𝟢 GHz, it would still take that massively overpowered machine \frac{8493465600}{4\cdot 10^9}=2.123 seconds. [1] Larger images and higher maximum iterations would only increase the generated overhead.
The program’s source code is listed below and can also be downloaded (.java), though the compiled .jar can also be downloaded.

mandelbrot-set-iii_3840x2160_Zre-1.7860326891575045_Zim-1.0938146063987995E-5_Zom4.328865274911016E-4_Clr256.png
Green self-similarity

Unrelated to algorithmically generating fractal renderings, I recently found a weed which seemed to be related to the Mandelbrot Set and makes nature’s intertwined relationship with fractals blatently obvious. I call it the Mandel Weed.

mandelbrot-set-iii_mandel-weed.jpg
Mandel Weed
Source code: MandelbrotSetExplorer.java
Extra assets: mandelbrot-set-iii_mk.sh, mandelbrot-set-iii_3840x2160_Zre-0.7451119128483528_Zim0.0980857735166215_Zom0.0015327239951060212_Clr256.png, mandelbrot-set-iii_3840x2160_Zre-1.3091727307741978_Zim0.0798169059094002_Zom4.0442023600132944E-8_Clr256.png, mandelbrot-set-iii_3840x2160_Zre-1.3203805004961733_Zim0.08763686694940123_Zom6.070130782934575E-9_Clr256.png, mandelbrot-set-iii_3840x2160_Zre-1.3716081094521764_Zim0.0758961984414791_Zom8.091648167718248E-4_Clr256.png, mandelbrot-set-iii_3840x2160_Zre-1.99181475237973_Zim-3.982534247685533E-7_Zom3.0401148534882E-6_Clr256.png, mandelbrot-set-iii_3840x2160_Zre-1.9963766117665036_Zim-8.944472584209038E-7_Zom7.062361041362908E-6_Clr256.png, mandelbrot-set-iii_3840x2160_Zre-1.996376645843924_Zim3.7516943208736117E-6_Zom2.859109719673806E-10_Clr256.png, mandelbrot-set-iii_3840x2160_Zre0.1907127228021151_Zim-0.5505294608801341_Zom5.308930362839995E-4_Clr256.png, mandelbrot-set-iii_3840x2160_Zre0.2428560033896521_Zim0.5073262577216794_Zom5.8089779367154003E-5_Clr256.png, mandelbrot-set-iii_3840x2160_Zre0.3411570576688668_Zim0.5750384378500207_Zom1.100339471773629E-4_Clr256.png, mandelbrot-set-iii_3840x2160_Zre0.35167270726867744_Zim-0.06184250438694921_Zom9.837555143550957E-5_Clr256.png, mandelbrot-set-iii_3840x2160_Zre0.35940345191970086_Zim-0.06373194819231794_Zom8.586182289196699E-7_Clr256.png, mandelbrot-set-iii_3840x2160_Zre0.4443477244180509_Zim-0.38358230098736174_Zom1.4540775100673537E-6_Clr256.png, mandelbrot-set-iii_3840x2160_Zre0.44503254319557867_Zim-0.38339448301161105_Zom2.4624930313254245E-6_Clr256.png

Asciify

2017-07-15, post № 175

art, ascii, programming, Python PIL, #ANSI, #ASCII, #characters, #image processing, #low resolution, #retro look

Most images nowadays are represented using pixels. They are square, often relatively small and numerous, come in (2^8)^3 different colors and thereby do a good job being the fundamental building block of images. But one can imagine more coarse-grained and differently shaped pixels.
An interesting fact is, that in most monotype fonts two characters placed right next to each other (for example ‘$$’) occupy roughly a square area. So simple ASCII characters can indeed be used to approximately describe any ordinary image.
Asciify does exactly this; it takes in an image and some optional parameters and maps the pixels’ intensity onto a character set. Both the large and small default character sets are taken from a post by Paul Bourke.

In conjunction with asciify.py, I wrote index.py, which asciifies a bunch of images and results in their html form; it also creates an index. All images asciified for this post can be viewed through this index [1].

Converting an image to its asciified form works best when there is a lot of contrast in the image. Because of this, some pre-processing of the image may be required for best results (all images shown where only cropped or rotated). The built-in color functionality also only knows of 𝟪 colors, so bright and different colors look the best, as they interestingly differentiate from one another. The asciified image’s size also plays a role, the larger it is, the better the characters blend into one and appear to be one image.

Asciify is operated on a command prompt; python asciify.py img.png. To parse arguments, the built-in Python module argparse is used. The images are opened and read using the Python Imaging Library module PIL, which needs to be installed for this program to work.
Optional arguments include --size N, where the maximum size can be specified, --invert and --smallcharset, which can sometimes increase the asciified image’s visual appeal and --html, which will output an html file to be viewed in a browser. To see the program’s full potential, simply run python asciify.py --help [2].
Source code for both asciify.py and index.py can be downloaded, the first is also listed below.

The two examples above use the color mode, though certain images also work in default black and white mode, such as this spider I photographed.

Then again, the colored text also has its charm, especially when the source image has bright colors and a lot of contrast.

Source code: asciify.py

Seventeen

2017-07-01, post № 174

mathematics, #17, #2017, #integer, #integer sequences, #July, #sequences

Today it is the first day of July in the year 2017. On this day there is a point in time which can be represented as 1.7.2017, 17:17:17.
To celebrate this symbolically speaking 17-heavy day, I created a list of 17 integer sequences which all contain the number 17.
All sequences were generated using a Python program; the source code can be viewed below or downloaded. Because the following list is formatted using LaTex, the program’s plaintext output can also be downloaded.

  1. Prime numbers 𝑛.

    \{2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, \dots\}
  2. Odd positive integers 𝑛 whose number of goldbach sums (all possible sums of two primes) of 𝑛 + 𝟣 and 𝑛 - 𝟣 are equal to one another.

    \{5, 7, 15, 17, 19, 23, 25, 35, 75, 117, 177, 207, 225, 237, 321, 393, 453, 495, 555, 567, \dots\}
  3. Positive integers n who are part of a Pythagorean triple excluding 𝟢: n^2=a^2+b^2 with integers a,b>0.

    \{5, 10, 13, 15, 17, 20, 25, 26, 29, 30, 34, 35, 37, 39, 40, 41, 45, 50, 51, 52, \dots\}
  4. Positive integers 𝑛 where \lfloor(n!)^{\frac{1}{n}}\rfloor is prime

    \{4, 5, 6, 7, 8, 12, 13, 17, 18, 19, 28, 29, 33, 34, 35, 44, 45, 46, 49, 50, \dots\}
  5. Positive integers 𝑛 with distance 𝟣 to a perfect square.

    \{1, 2, 3, 5, 8, 10, 15, 17, 24, 26, 35, 37, 48, 50, 63, 65, 80, 82, 99, 101, \dots\}
  6. Positive integers 𝑛 where the number of perfect squares including 𝟢 less than 𝑛 is prime.

    \{2, 3, 4, 5, 6, 7, 8, 9, 17, 18, 19, 20, 21, 22, 23, 24, 25, 37, 38, 39, \dots\}
  7. Prime numbers 𝑛 where either 𝑛 - 𝟤 or 𝑛 + 𝟤 (exclusive) are prime.

    \{3, 7, 11, 13, 17, 19, 29, 31, 41, 43, 59, 61, 71, 73, 101, 103, 107, 109, 137, 139, \dots\}
  8. Positive integers 𝑛 whose three-dimensional vector’s (n,n,n) floored length is prime, \lfloor\sqrt{3\cdot n^2}\rfloor is prime.

    \{2, 3, 8, 10, 11, 17, 18, 24, 25, 31, 39, 41, 46, 48, 60, 62, 63, 76, 91, 100, \dots\}
  9. Positive integers 𝑛 who are the sum of a perfect square and a perfect cube (excluding 𝟢).

    \{2, 5, 9, 10, 12, 17, 24, 26, 28, 31, 33, 36, 37, 43, 44, 50, 52, 57, 63, 65, \dots\}
  10. Positive integers 𝑛 whose decimal digit sum is the cube of a prime.

    \{8, 17, 26, 35, 44, 53, 62, 71, 80, 107, 116, 125, 134, 143, 152, 161, 170, 206, 215, 224, \dots\}
  11. Positive integers 𝑛 for which \text{decimal\_digitsum}(n)+n is a perfect square.

    \{2, 8, 17, 27, 38, 72, 86, 135, 161, 179, 216, 245, 275, 315, 347, 432, 467, 521, 558, 614, \dots\}
  12. Prime numbers 𝑛 for which \text{decimal\_digitsum}(n^4) is prime.

    \{2, 5, 7, 17, 23, 41, 47, 53, 67, 73, 97, 103, 113, 151, 157, 163, 173, 179, 197, 199, \dots\}
  13. Positive integers 𝑛 where \text{decimal\_digitsum}(2 \cdot n) is a substring of 𝑛.

    \{9, 17, 25, 52, 58, 66, 71, 85, 90, 104, 107, 115, 118, 123, 137, 142, 151, 156, 165, 170, \dots\}
  14. Positive integers 𝑛 whose decimal reverse is prime.

    \{2, 3, 5, 7, 11, 13, 14, 16, 17, 20, 30, 31, 32, 34, 35, 37, 38, 50, 70, 71, \dots\}
  15. Positive integers 𝑛 who are a decimal substring of n^n.

    \{1, 5, 6, 9, 10, 11, 16, 17, 19, 21, 24, 25, 28, 31, 32, 33, 35, 36, 37, 39, \dots\}
  16. Positive integers 𝑛 whose binary expansion has a prime number of 𝟣’s.

    \{3, 5, 6, 7, 9, 10, 11, 12, 13, 14, 17, 18, 19, 20, 21, 22, 24, 25, 26, 28, \dots\}
  17. Positive integers 𝑛 whose 7-segment representation uses a prime number of segments.

    \{1, 2, 3, 5, 7, 8, 12, 13, 15, 17, 20, 21, 26, 29, 30, 31, 36, 39, 47, 48, \dots\}
Source code: seventeen.py

Tau Day MMXVII

2017-06-28, post № 173

C, mathematics, programming, #constant

Today it is June the 28th which means that it is 𝜏 day!
The irrational and transcendental constant 𝜏 is what defines \pi=\frac{\tau}{2}, which obviously makes it an important constant. To celebrate this day, I created a C program which calculates 𝜏 by randomly creating 𝟫-dimensional points inside the 𝟫-dimensional hypercube and testing if they are inside the 𝟫-dimensional hypersphere with its center located at (0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5) [1].

Today’s 𝜏 time is 3:18:53 as \tau=6.2831853\dots. As one does not know if the time is specified as ante or post meridiem, there are actually two perfectly acceptable 𝜏 times.

            ;b$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$b
        h$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
     .$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$+ 
    $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$  
  ,$$$$$$$$$kn}~I"`($$$$$$$x````````````````   
 )$$$m             $$$$$$$$                    
+$$c              q$$$$$$$>                    
$$                $$$$$$$$                     
                 $$$$$$$$.                     
                '$$$$$$$$                      
                $$$$$$$$                       
               ?$$$$$$$a                       
               $$$$$$$$                        
              Q$$$$$$$f                        
              $$$$$$$$                         
             $$$$$$$$i                         
             $$$$$$$$                          
            $$$$$$$$;                          
            $$$$$$$$                           
           J$$$$$$$w                           
           $$$$$$$$;             $$~           
           $$$$$$$$             $$$            
          '$$$$$$$$C          "$$$             
           $$$$$$$$$w       '$$$$              
           $$$$$$$$$$$$$$$$$$$$B               
            $$$$$$$$$$$$$$$$$$.                
             $$$$$$$$$$$$$$$'                  
                8$$$$$$$q                      

The formula used for calculating 𝜏 is derived from a 𝟫-dimensional hypersphere’s hypervolume formula V=\frac{32\cdot\pi^4}{945}\cdot R^9 (see this Wikipedia article).

\begin{cases}
    V = \frac{2^5 \cdot \tau^4 \cdot R^9}{945 \cdot 2^4} = \frac{2 \cdot R^9}{945} \cdot \tau^4\\
    \tau^4 = \frac{V \cdot 945}{2 \cdot R^9}; R = 0.5\\
    \tau^4 = \frac{V \cdot 945 \cdot 2^9}{2}\\
    \tau = \sqrt[4]{V \cdot 241920}
\end{cases}

The constant gets calculated to \tau^*=6.293700. The real value is approximately 𝜏 = 𝟨.𝟤𝟪𝟥𝟣𝟪𝟧…, which makes the percent error

\left|\frac{\tau^*-\tau}{\tau}\right|=\left|\frac{6.293700-6.283185}{6.283185}\right|=0.001673514 = 0.167\%.

Thereby, this C program’s approximation is not too far off. [2] The source code is listed below and can also be downloaded here. Instructions on how to compile it using GCC can be seen below or in the source code.

gcc tau.c -o tau -lm; ./tau
tau = 6.293700

Resources worth checking out regarding 𝜏 are The Tau Manifesto and 2 Pi or Not 2 Pi? I wish everybody a happy 𝜏 day.

Source code: tau-day-mmxvii.c

Mandelbrot Set ASCII Viewer

2017-06-17, post № 172

ascii, HTML, JavaScript, mathematics, programming, Python, #animation, #fractal

The Mandelbrot Set is the set of all complex points which, when one iteratively and infinitely applies the function f_c(z)=z^2+c, converge to a value. This simple rule results in stunning complexity and beauty.
Many Mandelbrot Set animations use regularly colored pixels to represent the number of iterations needed at the fractal’s edges to escape converging. Yet this mathematical object can also be represented as ASCII characters — similar to what I did in my Curses Cam post. The characters are chosen according to their opaqueness. A full stop (‘.’) looks lighter than a dollar sign (‘$’), so they represent a smaller or larger number of iterations needed. The order of characters used is taken from this post by Paul Bourke.
As there are only 𝟩𝟢 characters used, each frame is being rendered twice to determine the minimum number of iterations needed by every point in that frame. Thereby the full visual character range is used.

The characters shown below represent a Mandelbrot Set still. To see the zoom in action, either run the program (listed below) or take a look at this Mandelbrot Set ASCII journey.

      ..................''''''''``"">>II``''''......                          
    ..................''''''''``^^,,ii::^^``''''......                        
  ..................''''''''``^^::ww$$++,,````''''......                      
................''''''''``^^^^""::$$$$$$::""^^``''''......                    
..............''''''````""{{;;XX$$$$$$$$uuUU,,,,""''......                    
............''''``````^^,,rr$$$$$$$$$$$$$$$$<<$$--``........                  
........''``````````^^""LL$$$$$$$$$$$$$$$$$$$$__""``''......                  
..''''''^^!!"""",,""""::__$$$$$$$$$$$$$$$$$$$$$$ll""''........                
''''````^^::__IIYYii::ll$$$$$$$$$$$$$$$$$$$$$$$$pp^^''........                
''``````"";;[[$$$$$$++__$$$$$$$$$$$$$$$$$$$$$$$$$$^^''''......                
``^^^^,,;;>>$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ww``''''......                
"",,,,II$$nn$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$""``''''......                
"",,,,II$$nn$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$""``''''......                
``^^^^,,;;>>$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ww``''''......                
''``````"";;[[$$$$$$++__$$$$$$$$$$$$$$$$$$$$$$$$$$^^''''......                
''''````^^::__IIYYii::ll$$$$$$$$$$$$$$$$$$$$$$$$pp^^''........                
..''''''^^!!"""",,""""::__$$$$$$$$$$$$$$$$$$$$$$ll""''........                
........''``````````^^""LL$$$$$$$$$$$$$$$$$$$$__""``''......                  
............''''``````^^,,rr$$$$$$$$$$$$$$$$<<$$--``........                  
..............''''''````""{{;;XX$$$$$$$$uuUU,,,,""''......                    
................''''''''``^^^^""::$$$$$$::""^^``''''......                    
  ..................''''''''``^^::ww$$++,,````''''......                      
    ..................''''''''``^^,,ii::^^``''''......                        

The fractal viewer is written in Python 2.7 and works by determining the terminal’s size and then printing a string of according size. This creates the illusion of a moving image, as the terminal will hopefully always perfectly scroll so that only one frame is visible at a time. In the code’s first non-comment line one can change the complex point at the image’s center, (really, its conjugate, which is partially irrelevant as the set is symmetric along the real axis) the initial zoom value (complex distance above the image’s center), the zoom factor (the factor by which the zoom value gets multiplied after a frame), the total number of frames ( - 𝟣 means there is no upper limit), the delay between frames (in seconds, can be floating-point) and the color characters used.

The program’s source code may not be particularly easy to read, yet it does its job and only requires seven non-comment lines! The code is shown below, though the .py file can also be downloaded.
To achieve the JavaScript animation linked to above, I wrote a simple Python converter which takes in the fractal renderer’s output and it spits out an HTML page. This converter’s code is not listed, though the .py file can be downloaded. Instructions on how to use the converter can be seen in its source code.

mandelbrot-set-ascii-viewer.py; Python2, 593 bytes, 8 lines
# Python 2.7 Code; Jonathan Frech, 15th and 16th of June 2017
P,Z,F,N,D,K=-.707+.353j,3,.9,-1,.1," .'`^\",:;Il!i><~+_-?][}{1)(|\\/tfjrxnuvczXYUJCLQ0OZmwqpdbkhao*#MW&8%B@$"
import os,time,sys;H,W,S,n=map(int,os.popen("stty size").read().split())+[sys.stdout,0];W/=2
def C(c):
	global m;z,i=0j,-1
	while abs(z)<=2 and i<len(K)-1+M:z,i=z*z+c,i+1
	m=min(m,i);return K[i-M]*2
while n<N or N==-1:h=Z*2.;w=h*W/H;R=lambda:"\n\n"*(n!=0)+"\n".join("".join(C(P-complex(w/2-w*x/W,h/2-h*y/H))for x in range(W))for y in range(H));M,m=0,len(K);R();M=max(M,m);S.write(R());S.flush();Z,n=Z*F,n+1;time.sleep(D)

A285494

2017-06-03, post № 171

mathematics, OEIS, programming, Python, Wolfram Language, #decimal, #decimal sum, #digit sum, #distinct prime, #int, #integer, #number, #number theory, #prime, #prime factors, #primes

The On-Line Encyclopedia of Integer Sequences gets regularly updated with new integer sequences. One of the recent updates was contributed by me, A285494.

A285494 is the list of all numbers 𝑘 so that its digit sum equals its number of distinct prime factors.
A number’s digit sum is the sum of all of its decimal digits. The number 𝟨𝟤𝟪𝟥𝟣𝟪𝟧𝟥, for example, has a digit sum of 𝟨 + 𝟤 + 𝟪 + 𝟥 + 𝟣 + 𝟪 + 𝟧 + 𝟥 = 𝟥𝟨.
A number’s number of distinct prime factors is the number of different prime numbers that multiply together to result in the original number. As an example, 62831853=3^2\cdot 7\cdot 127\cdot 7853, so it has five prime factors of which four are distinct.
Thereby one can conclude that 𝟨𝟤𝟪𝟥𝟣𝟪𝟧𝟥 is not an entry in this sequence, as 36\neq 4.

The sequence is certainly infinite, as the number k=2\cdot 10^n with n\in \mathbb{N}^* has a digit sum of 2+(0\cdot n)=2 and — because k=2^{n+1}\cdot 5^n — exactly two distinct prime factors.

In the encyclopedia entry, I provided a Mathematica one-liner to compute the first few entries of this sequence. Since then, I have also written a Python two-liner to achieve the same goal.

(* Mathematica *)
Select[Range[2,10000],Total[IntegerDigits[#]]==Length[FactorInteger[#]]&]
Out = {20, 30, 102, 120, 200, 300, 1002, 1200, 2000, 2001, 2002, 3000, 3010}
# Python 2.7
>>> def p(n):exec"i,f=2,set()\nwhile n>1:\n\tif n%i==0:f.add(i);n/=i;i=1\n\ti+=1";return len(f)
>>> print filter(lambda n:p(n)==sum(map(int,str(n))),range(2,10001))
[20, 30, 102, 120, 200, 300, 1002, 1200, 2000, 2001, 2002, 3000, 3010]