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, , so it has five prime factors of which four are distinct.
Thereby one can conclude that 𝟨𝟤𝟪𝟥𝟣𝟪𝟧𝟥 is not an entry in this sequence, as .

The sequence is certainly infinite, as the number with has a digit sum of and — because  — 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]

JSweeper

2017-05-20, post № 170

games, Java, programming, #2016, #classic, #clone, #mine, #Minesweeper, #pixel, #sweeping, #Windows

Adding to my collection of clones of popular, well-known games, I created back in November of 2016 a Java-implementation of the all-time Windows classic game, Minesweeper.

Minesweeper was pre-installed on every installation of Windows up to and including Windows 7 and has been ported to a variety of different systems. Because of this, nearly everyone has at least once in their life played Minesweeper or at least heard of it.
In Minesweeper you are presented with a square grid of covered tiles containing either numbers or mines. Your task is it to uncover all tiles which are not mines in the least amount of time. When you uncover a mine, it explodes and the game is lost. To aid in figuring out which tiles are mines and which are not, every tile that is not a mine tells you how many mines are in the neighbouring eight tiles. Tiles which have no neighbouring mines are drawn gray and uncover neighbouring non-mine tiles once uncovered.
More on Minesweeper can be found in this Wikipedia article — I am linking to the German version, as the current English version has major flaws and lacks crucial information. If you are so inclined, feel free to fix the English Minesweeper Wikipedia article.

In my clone, there are three pre-defined difficulty levels, directly ported from the original Minesweeper game, and an option to freely adjust the board’s width and height as well as the number of bombs which will be placed. Gameplay is nearly identical to the original, as my clone also uses a square grid and the tile’s numbers correspond to the number of bombs in the eight tiles surrounding that tile.
The game has a purposefully chosen pixel-look using a self-made font to go along with the pixel-style.

To play the game, you can either download the .jar file or compile the source code for yourself. The source code is listed below and can be downloaded as a .java file.

Pinhole Photographs MMXVII

2017-05-06, post № 169

art, #light, #nature, #photography, #picture, #tulip, #World-Wide Pinhole Day, #WWPD

Purple Tulip

Crimson Tulip

Light Interference

Multibrot Set

2017-04-22, post № 168

Java, mathematics, programming, #animated gif, #animation, #Cartesian, #complex, #complex arithmetic, #fractal, #generalization, #gif, #Mandelbrot set, #multi-threading, #polar, #reals, #threading

The Mandelbrot Set is typically defined as the set of all numbers for which — with , and  — the limit converges. Visualizations of this standard Mandelbrot Set can be seen in three of my posts (Mandelbrot Set, Mandelbrot Set Miscalculations and Mandelbrot Set II).

However, one can extend the fractal’s definition beyond only having the exponent 𝟤 in the function to be with ^[1]. The third post I mentioned actually has some generalization as it allows for , although the approach used cannot be extended to real or even rational numbers.

The method I used in the aforementioned post consists of manually expanding for each 𝑛. The polynomial , for example, would be expanded to .
This method is not only tedious, error-prone and has to be done for every exponent (of which there are many), it also only works for whole-number exponents. To visualize real Multibrots, I had to come up with an algorithm for complex number exponentiation.

Luckily enough, there are two main ways to represent a complex number, Cartesian form and polar form . Converting from Cartesian to polar form is simply done by finding the number’s vector’s magnitude and its angle to the 𝑥-axis . (The function is used in favor of to avoid having to divide by zero. View this Wikipedia article for more on the function and its definition.)
Once having converted the number to polar form, exponentiation becomes easy, as

With the exponentiated in polar form, it can be converted back in Cartesian form with

Using this method, converting the complex number to perform exponentiation, I wrote a Java program which visualizes the Multibrot for a given range of exponents and a number of frames.
Additionally, I added a new strategy for coloring the Multibrot Set, which consists of choosing a few anchor colors and then linearly interpolating the red, green and blue values. The resulting images have a reproducible (in contrast to randomly choosing colors) and more interesting (in contrast to only varying brightness) look.

The family of Multibrot Sets can also be visualized as an animation, showing the fractal with an increasing exponent. The animated gif shown below was created using ImageMagick’s convert -delay <ms> *.png multibrot.gif command to stitch together the various .png files the Java application creates. To speed up the rendering, a separate thread is created for each frame, often resulting in 𝟣𝟢𝟢％ CPU-usage. (Be aware of this should you render your own Multibrot Sets!)

To use the program on your own, either copy the source code listed below or download the .java file. The sections to change parameters or the color palette are clearly highlighted using block comments (simply search for /*).
To compile and execute the Java application, run (on Linux or MacOS) the command javac multibrot.java; java -Xmx4096m multibrot in the source code’s directory (-Xmx4096m tag optional, though for many frames at high quality it may be necessary as it allows Java^[2] to use more memory).
If you are a sole Windows user, I recommend installing the Windows 10 Bash Shell.

Easter MMXVII

2017-04-16, post № 167

art, ascii, haiku, poetry, #ascii egg, #celebration, #easter egg, #egg

                           vunnnnnnnxvvczYX                         
                       uuxrjjft///tfjrnuvvcXXUU                     
                    cuxxrf/\|(|||/tttrnxuuvcczYYJL                  
                 cuuvnxrjt/\\\/\//tffrxnuvvcccXXYYCLC               
               cccvuunxjrttttttjjfrrxnnnvvcczXzXUUUJL0              
             cczvvvuuuxnrxrrrjrxrxnnuuuvvcczXXXXYXUJJLQZ            
           czzzXzvccvnuunnnnnnnnuuuvuvvccczzzXYXYUYUJCLQO           
          zzzzXzccvvvvvuunnnnunuunnuuucvzczzXXYXUYJUJJCQ0Oq         
         XXXXzzzzcccvvvvuuvunnnnnuuuvvvcczzXzXYXXYYUUCCL0Qmp        
        XXXzXzzzzzcccccvuunnnnnnnnunvucvzzzXzYXXYUUUUJJCL0OZq       
       XYYYYXXzczzzcczcvuuunxxnxxxnuuuvzzXXXXYXUYUUUUUCJCQ0Owp      
      YYUYXYYXXXXzczczcvvuunxxrrrxnnuvuzzzzXXXYYUUUUUUJJCCL0Omd     
     YYUUYUYUXXXXczzzzcvvvunxrxrrrxxnuvcczzXzYYXYXUYUUJJLQQ0OZmk    
    YUUJJJYUYYXUXYXzczzczuuxrjrrrjrxxccccczzXXXYUYUYUJUCCCL00mwdw   
    UJCJJJUYXYXXXzzXXzzvvnxxrjjjrjxxuvvczXzXYYYXXUYUUUJCJLLL0Zmpb   
   CLLCJJJCUUUXUUXXzzzccvvunxxrrxrnuuvzzzXYXXYYYUUUYJUJJJJLL0OZwda  
   QLQLCJCUUUUUUYYXXXXXzvuuxxxrrnxxuuvvccXzXYUYUYYYUUUUCCCLLQ0Zmpk  
  00QQLLLCUJUUUUUYYYXzzzccvuuxxxxnnuvczzcXXXXUUUUYYJUYUUJCLL0QOwqk* 
  Z0OQQLLJCJUUUYUYYYXXXzzcvvnunnuuvvvvczzzXXXYYYYUYUYUJJCJLL0QOmqpa 
 0ZmZ00QLLCCJCJUUUUYYXXzzcccvvuuuucvczXXXYXYYUUUYYYUUUJJJJC0L0ZOqpkM
 ZOZO0QLLLCCCJJJUUYYXYXXYXzccvvccvczzzzXXXYYYYYUYYUYUUUUJJCCQQOZwdk*
 wwmZ0OQCQCCJJJJUUUYUYYYYXXzzzcccczXXXzXXXXYYXUYYYUUJUUJJCULQQ0Zmpbo
 qmZZZ00QQLCJUJUUJUUUUUYYXXXXzcYXXzXUXYXYYYXUUUYUUUUUUUJCJCCQQOZmqb#
 dqpZZ00QQLLLLCCJJCUUUYUYUYUXYXYYYYYYYYYYUYYUYUYYUUUYUUUJJLLQ0Omwqb*
qpqwZZOO0QLQCCCJCUJYUUUUUUYXXYYYUUUUYUYYUYJUYUYUYYYUUUJJJULLQ00mwqba
pddpwmZ000QQCQCCLCJJCJUUJUJUUUUUUYUYJUUUUUYUYJUJUUUJUJUUCCLLQ0OOmpba
qkppwwmOOO0LQLLCCCCJJJJJJUJJJUCJUUUUUUUUUJUJJUJUJUUJUCCCCCLQ00OZqqh*
 kdbpmmZZO00QQQLQLCCCJCCUCJJJUJCUJJUUUUUJJUUUUJJYUJUUJJJCLLL00OZwqk#
 ahbdqwmmOZ0Q00LLLJCCCJJJCJJJJJCUJJJYUUJUUUUUUUUUUUJUJCJCLLQOQZmqdh#
 oohdqwwZZOO0QQLLQLQCCCCLLCLCCJJJJCUJUJJJJUUUUYUJJJCCCCCCCQ0OOmmpk#*
 d*abdqqwmmOOZ00OQLQ0LCLQCCCJCLCCCJCCJJUUUJUJUUJUJUJCCCCLCQOOmwpdh# 
  M#hkbqwmZmZOOO0Q00L0QLCLLQLLQLJCLCCJJCJCUJJJJJJJJCJJLLCQ0ZmwqpkaM 
   oohkddqwmmmOZOOO0LQL0LLQCL0LCLLLLLCJCJLJJJUUJJJJCCCLQ00OZmwqba*  
   WM*akbqppwmwZOZOQOO000LQLQQLCCLLLCCCCJJCCCCJCJCLCCCLQ0OZmppboo   
    MWohhdpwqwwmwOZZZOOO00QLQLQQLLLLLCCCLLLCLCLCLLJLLQ0ZOmmwdbaa    
     &W#hkbpqwqwqZZZZZZO000QQQ0QLLQLCLLCQCLQLLLLCLLQQL0ZZmwpbkap    
      8&#ohkbddpqwwmmmZOZOOO0OO0Q0Q0QLLQQLQLLQQ0QQQ00ZOmmqbdha      
       88WMohobdbpqpwwmqZZmZZZOQZO000Q0OQQQ00QQQ000OZmmqppkk*       
        M%8M#*okkbppdqqppmwZmmmwOOZZOO00O0O0O0O0OmZqqwqdbkho        
          B@8&M**ohkbbdqppqwqwwwwwZZZmZmOOOOOZZmmwwppbbkka          
            @@8&MW*ookkkbdpdpdqpqqwqpwZwZwwmqwqpppdpdkk*            
              $$%%&WM#*oahkkhkddppdpbqwpppppbpdkkbbhho              
                @$@B8&WM**#aoaahbkkhkhkbbkkbhhkhkkab                
                   $@%%8&WWM#M#*o**oooaoooaaaaook                   
                       @%%W&WWW#*##o*oa#**#M*                       
                             #88&&%8%&                              

T-3PO — Tic-Tac-Toe Played Optimally

2017-04-08, post № 166

games, HTML, JavaScript, mathematics, programming, Python, #computer player, #game ai, #perfect, #perfect play, #three in a row, #three to win, #three

Tic-Tac-Toe, noughts and crosses, Xs and Os, three in a row or whatever you want to call it may be the simplest perfect information game that is enjoyable by humans. Two players set their pieces (X or O) on an 𝟥 ⨉ 𝟥 grid, alternating their turns. The first player to get three of their pieces in a line, wins. If no player succeeds to get a line, the game ends in a draw.

Tic-Tac-Toe’s simplicity may become clear, if you consider that skilled players — people who have played a few rounds — can reliably achieve a draw, thereby playing perfectly. Two perfect players playing Tic-Tac-Toe will — whoever starts — always tie, so one may call the game virtually pointless, due to there practically never being a winner.
Because of its simple rules and short maximal number of turns (nine) it is also a game that can be solved by a computer using brute-force and trees.

The first Tic-Tac-Toe-playing program I wrote is a Python shell script. It lets you, the human player, make the first move and then calculates the best possible move for itself, leading to it never loosing. On its way it has a little chat whilst pretending to think about its next move. The Python source code can be seen below or downloaded here.

The second Tic-Tac-Toe-playing program I wrote uses the exact same method of optimizing its play, though it lets you decide who should begin and is entirely written in JavaScript. You can play against it by following this link.

Both programs look at the entire space of possible games based on the current board’s status, assumes you want to win and randomly picks between the moves that either lead to a win for the computer or to a draw. I did not include random mistakes to give the human player any chance of winning against the computer. Other Tic-Tac-Toe-playing computers, such as Google’s (just google the game^[1]), have this functionality.

Second Anniversary

2017-03-28, post № 165

art, #2 years, #celebration, #collage, #j, #J-Blog, #jblog, #two years

J-Blog celebrates its second anniversary!
Exactly two years ago, on the twenty-eighth of March 2015, the very first post on this blog appeared, appropriately named “Hello World”. Since then, including this one, 𝟣𝟨𝟦 other posts have been posted. Here a few of the image highlights from both years, with their corresponding post.

                     -   
            ?][)0Zqpkdc  
         QmphM&88888%av  
        Z&88%%%8M*%B8X   
       JW&Mbw\1] ZB@b    
       ?-        #BB%c   
                 m@@@kc  
                 ?o@BBC  
                  ]&@@*v 
                   L%@%Y 
                    o@BZ 
                    LB@a 
                    ]%BM?
                     *B&v
                     k%&v
                     mBWv
                     wBWv
                     ZB* 
                     w@O 
                     W8z 
                    [BO  
                    ao]  
Ckc                k&z   
b@#Z?            [a8L    
X%$@*Q[?      ?[q%&}     
 {&%%8&kwm00mdM%%w?      
  ?ZW8%%%%%%B8&q]        
    -}\wpkpft[           

Rectangles

Prime Spiral

Gradient Drawer

Pattern

Stacking Stones

Pastel

White Fireworks

Jappy Jird

Curses Cam

Maze Solving

Bifurcation Diagram

2017-03-25, post № 164

mathematics, PIL, programming, Python, #alpha, #chaos, #chaos theory, #delta, #Feigenbaum, #fractal, #iterations, #Mandelbrot set, #modelling, #population

Generating the famous fractal, which can be used to model populations with various cycles, generate pseudo-random numbers and determine one of nature’s fundamental constants, the Feigenbaum constant 𝛿.
The fractal nature comes from iteratively applying a simple function, with , and looking at its poles.
The resulting image looks mundane at first, when looking at , though the last quarter section is where the interesting things are happening (hence the image below only shows the diagram for ).
From 𝜆 = 𝟥 on, the diagram bifurcates, always doubling its number of poles, until it enters the beautiful realm of chaos and fractals.

For more on bifurcation, fractals and 𝛿, I refer to this Wikipedia entry and WolframMathworld.