Best Computed Graphics

最高の計算されたグラフィックス

Line wrap

Highlight

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#define P printf
#define F(x,y) ((FILE*)Z(x,y))
#define U(x,y) ((unsigned char*)Z(x,y))
#define X(x,y,z) ((double*)x+y+7*z)
#define Y(x,y) ((double*)x+y)
#define Z(x,y) *(long*)Y(x,y)
#define R(x,y) (x=y,y^=(y<<11),y^=(y>>15),y^=(y<<3),(x+y)&0x7fffffff)

int main(int c,char**o)
{
   void* v=malloc(0170UL);
   memset(v,0,0170UL);
   if(c!=4)P("bad args\n"),exit(1);
   Z(v,2)=atoi(o[1]);
   Z(v,1)=atoi(o[2]);
   Z(v,0)=Z(v,2)*Z(v,2)>>3;
   Z(v,6)=(long)malloc(Z(v,0));
   memset((void*)Z(v,6),0,Z(v,0));
   Z(v,13)=(long)fopen(o[3],"r");
   fscanf(F(v,13),"%ld",(long*)Y(v,3));
   Z(v,8)=(long)malloc(050+Z(v,3)*070);
   Z(Z(v,8),0)=Z(v,3);
   fscanf(F(v,13),"%lf %lf %lf %lf",Y(Z(v,8),1),Y(Z(v,8),2),Y(Z(v,8),3),Y(Z(v,8),4));
   for(Z(v,3)=0;Z(v,3)<Z(Z(v,8),0);
       fscanf(F(v,13),"%lf %lf %lf %lf %lf %lf %lf",X(Z(v,8),5,Z(v,3)),
              X(Z(v,8),6,Z(v,3)),X(Z(v,8),7,Z(v,3)),X(Z(v,8),8,Z(v,3)),
              X(Z(v,8),9,Z(v,3)),X(Z(v,8),10,Z(v,3)),X(Z(v,8),11,Z(v,3))),
           ++Z(v,3));
   fclose(F(v,13));
   Z(v,14)=8675309;
   for(Z(v,3)=0;Z(v,3)<Z(v,1);++Z(v,3))
   {
      *Y(v,7)=(R(Z(v,4),Z(v,14))/(double)0x7fffffff),Z(v,5)=0;
      for(Z(v,4)=0;Z(v,4)<Z(Z(v,8),0);++Z(v,4),++Z(v,5))
         if(*Y(v,7)<*X(Z(v,8),11,Z(v,5)))
            break;
         else
            *Y(v,7)-=*X(Z(v,8),11,Z(v,5));
      *Y(v,11)=*X(Z(v,8),5,Z(v,5))**Y(v,9)+*X(Z(v,8),6,Z(v,5))**Y(v,10)+*X(Z(v,8),9,Z(v,5));
      *Y(v,12)=*X(Z(v,8),7,Z(v,5))**Y(v,9)+*X(Z(v,8),8,Z(v,5))**Y(v,10)+*X(Z(v,8),10,Z(v,5));
      *Y(v,9)=*Y(v,11);
      *Y(v,10)=*Y(v,12);
      if(Z(v,3)>10)
         if((*Y(v,9)>=*Y(Z(v,8),1))&&(*Y(v,9)<=*Y(Z(v,8),3))&&(*Y(v,10)>=*Y(Z(v,8),2))&&(*Y(v,10)<=*Y(Z(v,8),4)))
            *Y(v,11)=(*Y(v,9)+-*Y(Z(v,8),1))/(*Y(Z(v,8),3)-*Y(Z(v,8),1)),
      *Y(v,12)=(*Y(v,10)+-*Y(Z(v,8),2))/(*Y(Z(v,8),4)-*Y(Z(v,8),2)),
      *(U(v,6)+((((long)(Z(v,2)-*Y(v,12)*Z(v,2))*Z(v,2))+(long)(*Y(v,11)*Z(v,2)))>>3))|=
      1<<((((long)(Z(v,2)-*Y(v,12)*Z(v,2))*Z(v,2))+(long)(*Y(v,11)*Z(v,2)))&7);
   }
   P("#define x_width %ld\n#define x_height %ld\nstatic char x_bits[] = {\n",Z(v,2),Z(v,2));
   for(Z(v,3)=0;Z(v,3)<Z(v,0);++Z(v,3))
   {
      if(Z(v,3))
         P(",");
      P("0x%02x",*(U(v,6)+Z(v,3)));
   }
   P("};\n");
   free(Y(Z(v,6),0));
   free(Y(Z(v,8),0));
   free(v);
   return 0;
}

引用元：https://www.ioccc.org/2006/stewart/stewart.c

審査員・作者による説明：https://www.ioccc.org/2006/stewart/hint.text

動作

フラクタルの一種である反復関数系の描画を行う。

シェルピンスキーのギャスケット。

$ gcc -o stewart stewart.c

$ ./stewart 400 1000000 gasket > gasket.xbm

$ convert gasket.xbm gasket.png

ドラゴン曲線。

$ ./stewart 400 1000000 dragon > dragon.xbm

$ convert dragon.xbm dragon.png

リアルなシダ。

$ ./stewart 400 1000000 fern > fern.xbm

$ convert fern.xbm fern.png

リアルな木。

$ ./stewart 400 1000000 tree2 > tree2.xbm

$ convert tree2.xbm tree2.png

IOCCCのロゴを目指したもの。

$ ./stewart 400 1000000 ioccc > ioccc.xbm

$ convert ioccc.xbm ioccc.png

他にも多数の設定が添付されている。

解説

[[1992/imc]]や[[1994/tvr]]など、フラクタルを生成する作品はいくつも実績があるけれど、反復関数系（IFS）を扱うものは初とのこと。

設定ファイルのフォーマットは次の通り。

<num ifs> <x min> <y min> <x max> <y max>
<a> <b> <c> <d> <e> <f> <p>
<a> <b> <c> <d> <e> <f> <p>
<a> <b> <c> <d> <e> <f> <p>
.
.
.

aからfは、(prevx, prevy)から(x, y)に点を移すときの式の座標。

x = a * prevx + b * prevy + e;
y = c * prevx + d * prevy + f;

最後のpは、その式が実行される確率とのこと。

画像出力はXBMという古いフォーマット。

コードは、変数宣言がvoid*の1つだけで、long型の値もdouble型の値もポインタもすべてこれにキャストして扱うようになっている。実質的に、型がない状態。プリプロセスをするとキャストがたくさん出てくるので、みるのが大変。

賞名は[[2006/monge]]と対称的になっている。