シェルピンスキーカーペットの反復を作成する方法

Question 1

適切な場所に白いノードを配置するマクロを次に示します。

\sierpinski[options]{levels}

私が考えているオプションはですscale=が、これrotate=も同様に機能します。または、xscale=正方形ではない長方形が必要な場合。必ずを使用してくださいtransform shape。

\documentclass{article}

\usepackage{tikz}

\newcommand{\sierpinski}[2][]{\tikz[#1]{
  \draw[fill=black] rectangle(1,1);
    \foreach \n[evaluate=\n as \m using \n-1, evaluate=\n as \s using 1/3^\n, evaluate=\m as \p using 3^\m] in {1,...,#2}{
      \foreach \k[evaluate=\k as \x using (2*\k-1)/2/3^\m] in {1,...,\p}{
        \foreach \j[evaluate=\j as \y using (2*\j-1)/2/3^\m] in {1,...,\p}{
          \node[fill=white, minimum size=\s cm, inner sep=0] at (\x,\y){};
}}}}}

\begin{document}

\tikz{\draw[fill=black] rectangle(1,1);}\quad\sierpinski{1}\quad\sierpinski{2}\quad\sierpinski{3}\quad\sierpinski{4}

\end{document}

これが\sierpinski[scale=3, transform shape]{5}私のマシンが処理できるほぼすべてです。O(9ⁿ) なので、待つ準備をしてください。

Answer

適切な場所に白いノードを配置するマクロを次に示します。

\sierpinski[options]{levels}

私が考えているオプションはですscale=が、これrotate=も同様に機能します。または、xscale=正方形ではない長方形が必要な場合。必ずを使用してくださいtransform shape。

\documentclass{article}

\usepackage{tikz}

\newcommand{\sierpinski}[2][]{\tikz[#1]{
  \draw[fill=black] rectangle(1,1);
    \foreach \n[evaluate=\n as \m using \n-1, evaluate=\n as \s using 1/3^\n, evaluate=\m as \p using 3^\m] in {1,...,#2}{
      \foreach \k[evaluate=\k as \x using (2*\k-1)/2/3^\m] in {1,...,\p}{
        \foreach \j[evaluate=\j as \y using (2*\j-1)/2/3^\m] in {1,...,\p}{
          \node[fill=white, minimum size=\s cm, inner sep=0] at (\x,\y){};
}}}}}

\begin{document}

\tikz{\draw[fill=black] rectangle(1,1);}\quad\sierpinski{1}\quad\sierpinski{2}\quad\sierpinski{3}\quad\sierpinski{4}

\end{document}

これが\sierpinski[scale=3, transform shape]{5}私のマシンが処理できるほぼすべてです。O(9ⁿ) なので、待つ準備をしてください。

Question 2

最新のhttps://ctan.org/pkg/pst-fractal実行中 lualatex：

\documentclass[pstricks]{standalone}
\usepackage{pst-fractal} 
\begin{document}    

\begin{pspicture}(18,3)
\multido{\iA=1+1,\iB=0+12}{6}{%
  \psSierCarpet[scale=0.25,n=\iA](\iB,0.2)}
\end{pspicture}

\end{document}

オプション付きbasecolor=red,linecolor=blue:

Answer

最新のhttps://ctan.org/pkg/pst-fractal実行中 lualatex：

\documentclass[pstricks]{standalone}
\usepackage{pst-fractal} 
\begin{document}    

\begin{pspicture}(18,3)
\multido{\iA=1+1,\iB=0+12}{6}{%
  \psSierCarpet[scale=0.25,n=\iA](\iB,0.2)}
\end{pspicture}

\end{document}

オプション付きbasecolor=red,linecolor=blue:

Question 3

探していた tikz の回答が得られましたので、MetaPost/MetaFun で作成した回答を自由に追加します。おそらくコードは最適化でき、おそらく tikz でも同様のことができるでしょう (しかし、私にはできません)。最初のバージョンでは、正方形を描画してループし、白であるべき部分を塗りつぶしません。2 番目では再帰を使用します。

タイミングについて：

First version: 2.9s
Second version: 1.8s

私はそれを ConTeXt lmtx の MetaPost ページにラップしました。ファイルはでコンパイルできますcontext。

\starttext
\startMPpage[offset=1dk]
vardef sierpinski(expr w, n) =
image(
fill unitsquare scaled w ;
for i = 1 upto n :
    for j = 1 upto (3^(i-1)) :
        for k = 1 upto (3^(i-1)) :
            unfill unitsquare scaled (w/(3^i)) shifted ( (3*j-2)*w/(3^i), (3*k-2)*w/(3^i) ) ;
        endfor
    endfor
endfor)
enddef ;

for i = 1 upto 3 :
    draw sierpinski(3cm,i)   shifted (4*(i-1)*cm,  0  ) ;
    draw sierpinski(3cm,i+3) shifted (4*(i-1)*cm, -4cm) ;
endfor ;
\stopMPpage

\startMPpage[offset=1dk]
vardef Sierpinski(expr w,n) =
    save tmppic ;
    picture tmppic ;
    if n = 1 :
        image(
        fill unitsquare scaled w ;
        unfill unitsquare scaled (w/3) shifted (w/3,w/3) ;
        )
    else :
        tmppic := Sierpinski(w, n - 1) scaled 1/3 ;
        image(
        for i = 1 upto 3 :
            for j = 1 upto 3 :
                if ((i*j) <> 4) :
                    draw tmppic shifted (((i-1)/3)*w,((j-1)/3)*w) ;
                fi
            endfor
        endfor
        ) 
    fi
enddef ;

for i = 1 upto 3 :
    draw Sierpinski(3cm,i)   shifted (4*(i-1)*cm,  0  ) ;
    draw Sierpinski(3cm,i+3) shifted (4*(i-1)*cm, -4cm) ;
endfor ;
\stopMPpage
\stoptext

出力は両方のバリエーションで同じように見えるので (私が見る限り)、そのうちの 1 つだけを示します。

Answer