[Tex/LaTex] drawing patterns based on prime factorization

diagramspstrickstikz-pgf

can someone help me generate code to type the diagram for any given natural number based on its prime factorization such as these diagrams

consider the next number 36 and its prime factorization, $36=3^2\cdot 2^2$, writing the primes in descending order… then, start with the pattern for a 3, which is 3 circles 120 degrees apart, now, replace each circle with 3 more circles in the same pattern to account for the $3^2$, this part would look like the pattern for 9 which is also $3^2$. now.. to finish the pattern for 36, we take the pattern for 9 and in each circle replace it with the pattern for 4… that is each of the 9 circles gets replaced with 4 tinier circles arranged in the same pattern as the pattern for 4

Best Answer

(See update at the end for an alternative output)

Done with Lualatex and tikz, simply converting the code from this javascript fiddle

This is a sample result for number 36:

Which was generated with:

\begin{tikzpicture}[y=-1cm]
\draw[red] (0,0) rectangle (5,5);
\primediagram{36}{2.5}{2.5}{2.5};
\end{tikzpicture}

The macro \primediagram receives four parameters. The first one is the number to represent. The remaining three are the center and radius of the diagram to be generated.

A more interesting example:

Produced by a tikz loop:

\begin{tikzpicture}[y=-1cm]
\foreach \y in {0,...,6} {
  \foreach \x in {0,...,4} {
    \draw(3.1*\x,3.1*\y) rectangle +(3,3);
    \pgfmathtruncatemacro{\number}{\x+5*\y}
    \primediagram{\number}{3.1*\x+1.55}{3.1*\y+1.55}{1.3};
    \node[black!50, below right] at (3.1*\x, 3.1*\y)  {\number};
  }
}
\end{tikzpicture}

The complete code

To produce these figures, three files are required:

primediagram.sty

This one simply defines the tex macro which interfaces with the lua code:

% This is primediagram.sty
\directlua{dofile("primediagram.lua")}
\newcommand{\primediagram}[4]{
    \directlua{draw(#1,#2,#3,#4)}
}

primediagram.lua

This one contains the lua code which actually performs the computations and outputs the required \draw commands for tikz. It is a direct translation to lua of the mentioned javascript code.

-- This is primediagram.lua
local smallfirst = false
local off2 = 0

function circle(cx, cy, s) 
    tex.print(string.format("\\fill (%f, %f) circle(%f);", cx, cy, s))
end

function draw(N, cx, cy, s)
  if N==0 then
    return
  end
  if N==1 then
     circle(cx,cy,s)
  else
    local f, r, d, oy, x, y
    f = primefactor(N)
    if f == 2 then
        oy = 0;
        if N % 4 == 0 then
            f = 4;
            r = 2 * s / (f + 2);
            d = f * s / (f + 2);
        else 
            f = 2;
            r = 0.75 * 2 * s / (2 + 2);
            d = 2 * s / (2 + 2);
        end
    else
        r = 2 * s / (f + 2);
        d = f * s / (f + 2);
        oy = d / 2 * (1 - math.cos(math.pi / f));
    end
    for i = 0, f do
        x = math.sin(math.pi + 2 * math.pi * (i + 0.5) / f + off2);
        y = math.cos(math.pi + 2 * math.pi * (i + 0.5) / f + off2);
        draw(N / f, cx + x * d, cy - y * d + oy, r);
    end
  end
end

local primes = {2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 
    47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 101, 103, 107, 109, 
    113, 127, 131, 137, 139, 149, 151, 157, 163, 167, 173, 179, 181, 
    191, 193, 197, 199, 211, 223, 227, 229, 233, 239, 241, 251, 257, 
    263, 269, 271}

function primefactor(N) 
    local ans = N;
    for i,pi in pairs(primes) do
        if N % pi == 0 then
            ans = pi;
            if smallfirst then
                return ans
            end
        end
    end
    return ans
end

The main document

You only need to include the sty package, as for example:

\documentclass{article}
\usepackage{nopageno}
\usepackage[margin=1cm]{geometry}
\usepackage{tikz}
\usepackage{primediagram}

\begin{document}
\begin{tikzpicture}[y=-1cm]
\foreach \y in {0,...,6} {
  \foreach \x in {0,...,4} {
    \draw(3.1*\x,3.1*\y) rectangle +(3,3);
    \pgfmathtruncatemacro{\number}{\x+5*\y}
    \primediagram{\number}{3.1*\x+1.55}{3.1*\y+1.55}{1.3};
    \node[black!50, below right] at (3.1*\x, 3.1*\y)  {\number};
  }
}
\end{tikzpicture}
\end{document}

Which produces the figure already shown.

Update: Rotating sub-groups

If each sub-groups is rotated according to the angle at which it is drawn, the output is different, nicer, and more close to the one shown by the OP. I also removed the black filling of the dots:

New result

This is the new code (only primediagram.lua has to be replaced):

-- This is primediagram.lua
local smallfirst = false
local off2 = 0

function circle(cx, cy, s) 
    tex.print(string.format("\\draw (%f, %f) circle(%f);", cx, cy, s))
end

function draw(N, cx, cy, s, a)
  if N==0 then
    return
  end
  if N==1 then
     circle(cx,cy,s)
  else
    local f, r, d, x, y, off
    off = 0
    a = a or 0
    tex.print(string.format("\\begin{scope}[xshift=%fcm, yshift=-%fcm, rotate=%f]", cx, cy, a))
    f = primefactor(N)
    if f == 2 then
        oy = 0;
        if N % 4 == 0 then
            f = 4;
            r = 2 * s / (f + 2);
            d = f * s / (f + 2);
        else 
            f = 2;
            r = 0.75 * 2 * s / (2 + 2);
            d = 2 * s / (2 + 2);
            off = math.pi/2
        end
    else
        r = 2 * s / (f + 2);
        d = f * s / (f + 2);
    end
    for i = 0, f-1 do
        local angle = math.pi + 2 * math.pi * (i + 0.5) / f + off
        x = math.sin(angle);
        y = math.cos(angle);
        draw(N / f, x * d, y * d , r, angle*180/math.pi+180);
    end
    tex.print("\\end{scope}")
  end
end

local primes = {2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 
    47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 101, 103, 107, 109, 
    113, 127, 131, 137, 139, 149, 151, 157, 163, 167, 173, 179, 181, 
    191, 193, 197, 199, 211, 223, 227, 229, 233, 239, 241, 251, 257, 
    263, 269, 271}

function primefactor(N) 
    local ans = N;
    for i,pi in pairs(primes) do
        if N % pi == 0 then
            ans = pi;
            if smallfirst then
                return ans
            end
        end
    end
    return ans
end

Related Solutions

[Tex/LaTex] How to generate all possible Venn diagrams (with the case below) efficiently

Right, here's some code:

\documentclass{standalone}
%\url{https://tex.stackexchange.com/q/67395/86}
\usepackage{tikz}

\makeatletter

\def\venn@strip#1#2\venn@STOP{%
  \def\venn@next{#1}%
  \gdef\venn@rest{#2}%
}

\newcommand{\venn}[1]{%
\begin{tikzpicture}
\coordinate (A) at (0,0);
\coordinate (B) at (2,0);
\coordinate (C) at (1,{sqrt(3)});
\coordinate (S-SE) at (5,-3);
\coordinate (S-NW) at (-3,{sqrt(3)+3});
  \edef\venn@rest{#100000000}%
  \foreach \i in {0,...,7} {
  \begin{scope}[even odd rule]
    \expandafter\venn@strip\venn@rest\venn@STOP
    \ifnum\venn@next=1\relax
    \pgfmathparse{Mod(\i,2) == 1 ? "(S-SE) rectangle (S-NW)" : ""}
    \path[clip] \pgfmathresult (A) circle[radius=2];
    \pgfmathparse{Mod(floor(\i/2),2) == 1 ? "(S-SE) rectangle (S-NW)" : ""}
    \path[clip] \pgfmathresult (B) circle[radius=2];
    \pgfmathparse{Mod(floor(\i/4),2) == 1 ? "(S-SE) rectangle (S-NW)" : ""}
    \path[clip] \pgfmathresult (C) circle[radius=2];
    \fill[rounded corners,red] (S-SE) rectangle (S-NW);
    \fi
  \end{scope}
  }
    \draw[ultra thick] (A) circle[radius=2];
    \draw[ultra thick] (B) circle[radius=2];
    \draw[ultra thick] (C) circle[radius=2];
    \draw[ultra thick,rounded corners] (S-SE) rectangle (S-NW);
\end{tikzpicture}
}

\makeatother

\newcommand{\allvendiagrams}{
% To generate the lot:
\foreach \j in {0,...,255} {
  \def\venncode{}
  \foreach \k in {0,...,7} {
    \pgfmathparse{Mod(floor(\j/2^\k),2) == 1 ? "\venncode1" : "\venncode0"}
    \global\let\venncode=\pgfmathresult
  }
  \venn{\venncode}

}
}

\begin{document}
\venn{10000000}
\venn{01000000}
\venn{11000000}
\end{document}

And here's the result:

Sample venn diagrams

I almost certainly have used a different code for the different regions - I went for simpler code. The rubric is that a 1 in the kth place fills the kth region, and the correspondence between labels and regions is to write out k as a binary number, then if the bit is set, that circle is used inside and if not, outside. At the end, we draw the region and circles on top. I haven't gone for much customisability, but hopefully it's fairly obvious what to change to get it to look different.

Edit from the questioner: By using \documentclass[border=3pt,tikz]{standalone} and invoking \allvendiagrams instead of \venn{10000000}\venn{01000000}\venn{11000000}, the output will show all the Venn diagrams as follows. But not in GIF for sure.

enter image description here

[Tex/LaTex] Mindmap sibling angle in Tikz

The style set angles for level is used to append the correct sibling distance for level #1.

It can either used inside the level style or possibly in a .listed option.

Usage should vary depending on the use-case in a real document. It may preferable to append set angles for level to grow cyclic or create a custom special grow cyclic that used this solution.

Code

\documentclass[tikz,convert=false]{standalone}
\usetikzlibrary{mindmap,trees,graphs}
\tikzset{
  set angles for level/.style={level #1/.append style={sibling angle=360/\the\tikznumberofchildren}},
  level/.append style={set angles for level=#1}% solution 1
}
\begin{document}
\begin{tikzpicture}
  \path[
    mindmap,
    concept color=black,
    text=white,
    grow cyclic,
    nodes=concept,
%    set angles for level/.list={1,...,4}% solution 2
  ]
  node {Computer Science} 
    child[concept color=green!50!black] {%
      node {practical}
      child {node {algorithms} }
      child {node {data structures} }
      child {node {pro\-gramming languages} }
      child {node {software engineer\-ing} }
      child {node {and a fifth}}
    }  
    child[concept color=blue] { node[concept] {applied}
      child {node {databases} }
      child {node {WWW} }
    }
    child[concept color=red] {node {technical} }
    child[concept color=orange] {node {theoretical} };
\end{tikzpicture}
\end{document}

Output