Best of show

最優秀賞

Line wrap

Highlight

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>

#define H 20

#define X for
#define F(v,b) X(k = 0; k < 3; k++) v[k] = b;

typedef double R, V[3], U[4], (*J)(R);
typedef void Z;
typedef int _;

V o, d, v;
R (*q)[H][6], C[4*H][H][6], sr = 15, br = 25, A, s, t;
_ P, Q, g, S, k, e = 256, (*u)(_) = putchar, D[64];
J f = cos, Y = sin;
char B[256];

Z
T(_ *p, _ k, _ n)
{
	_ g;

	*p = n;

	X(g = 0; k && g < 4; g++)
		T(p + (9 & (g * 4|g))*(1<<k - 1), k - 1,
		  n + (1<<6-2 * k)*(4 - g&3));
}

Z
O(_ c)
{
	B[++S] = c,

	(S > e-2 || c==129) && (*B = S, S = fwrite(B, S+1, 1, stdout) - 1);
}





    
					                       Z N
					                  (R*v){R w=sqrt
					                (*v**v +1[v]*v[1]
					              +v[2]*2[v]);F(v,v[k]
					             /w)}Z         K(_ j){V
					             n,p            ;R*b=(*
					            q)              [j],w=f
					                             (A),x=
					                             Y(A);
					                             F(p,
					                           br*o[
						                 k]+sr
		      * (n[                  k]               =f(t)*k[v]+
		Y(t)*(v[(k+1)%3]*           d[(            k+2)%3]-v[(k+2)%
	     3]*d[(k+1)%3])))*b++=p[1      ]+e/                2,*b++=-*p*x+
	   p[2]*w-e*e,*b++=(*p*w+p[2]*x)+e/2;F(                    b,fabs(*n
	  *!!k+n[1]*(        k<2)+n[2])/sqrt(2+                      !(k-1)))
	}Z M(R*a,_                s,Z(*m)(_)){_                        i;R d=
       pow(25/log                   (36),4)/s/                         377;X(
     *a=i=0;i<s                      ;i++,*a+=d                        )m(i);
      }Z G(_ i                         ){R w=-Y                        (P*s),
     x=f(P*s),                          y=Y(Q*s                        ),z=f
     (Q*s);F(                            v,k[o]                       =(k<2
    )*(x+3)*(                            k?y:z)     -!(k-2          )*w)N(
    v);F(d,P*                             (k?k-    1?x:w*y:w*z   )+Q*(k?
    k-1?0:*o:-                            o[1])     )N(d),w=*d**v+d[1]
    *v[1]+d[2]                             *v[2       ];F(v,k[v]-w
    *d[k])N(v),q                           =&C[
    i],M(&t,H,K);                          }Z W
     (char*s){*s&&(                         W(s
     +1),u(*s-98*(*s>
      97)));}Z E(_ z){R
       *r,*a,*b;U*q,l,x,d
	,I[256];_ i,j,m,c,y,
	 w,h,Y=sizeof l;fputc(
	  46,stderr),M(&s,4*H,G),
	    W("ibcbcbbbbb,"),g=130,
	     O(e/2);X(y=0;y<e;y++){X(i
		=e;i;)*(I[--i])=0;X(;i<4
		  *H*H;i++)X(*l=k=0;k-5;k++
		     ,a=b)if(b=C[(i/H+((k^k/
		       2)&1))%(4*H)][(i+(k/2&1
			 ))%H],k&&y<*a^y<*b){X(
			    h=0;h<4;h++)h[x]=a[h
			      +2]+(b[h+2]-a[h+2])
				*(y-*a)/(*b-*a),h?
				   *l&&(d[h]=(l[h]
     -h                             [x])/w):(w=1+
     fabs                            (*l-(*x=(_)*x)
     ) );                              X(q=I+(_)*x;
     *l&&                                w--;q+=2*(
     *x<*                                l)-1)X(**q
     >C[i                                /H][i%H][1
     ]&&(                                memcpy(*q,
     x,Y),*                              *q=C[i/H][
     i%H][1                              ]),h=1;h<
     4;h++)x                            [h]+=h[d];
     memcpy(l                           ,x,Y);}X(
     i=0;i<e;i++                        ,O(*r?c:
     127),g=g<e-1                    ?g+1:(O(e/
     2),130))X(r=&I[               i][3],c=0,j
     =2;j+1;j--,r--)k=3|!(    j-1)*4,m=*r*k+
     (63*( *r*k-(_)(*r*k))>D[(y*8&56)+(i&7)
     ]),        c<<=k/4+2,c|=(m|-(m>k))&
     k;}            O(129);u(0);}_
    
    
	
main(_ c, char **v)
{
	if (c < 3)
		return fprintf(stderr, "Usage: %s P Q [F]\n", *v), 1;

	P = atoi(v[1]),
	Q = atoi(v[2]),

	T(D, 3, 0);

	W("bb\346cbcba98FIG");

	X(; k < 384; k++)
		u(k / 3<<("gdb"[k%3] - 97) | 31 | (k % 3!=1)<<5);

	W("bbbce0.2EPACSTEN\x0b\xff!"),

	M(&A, c > 3 ? atoi(v[3]) : 40, E),

	u(59);
	
	return 0;
}

/* vi:set ts=8: */

引用元：https://www.ioccc.org/2005/persano/persano.c

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

動作

トーラス結び目のアニメGIFを生成する。

$ gcc -o persano persano.c -lm

$ ./persano 5 2 > knot5-2.gif
........................................

解説

./persano p qで(p,q)型トーラス結び目の画像を生成する。 (p,q)型トーラス結び目とは、ドーナツ形状（トーラス）の穴の周りをp周しつつ、ドーナツの断面方向にq周して元の位置に戻るときの軌跡のような形。

生成されるGIFはLZW圧縮をしていないとのこと。無圧縮GIFのようなフォーマットになっていると思われる。

コード形状は”S^3”。4次元空間における3次元球面のこと。 typedef double V[3],U[4]で、3次元と4次元のベクトルの型を定義している（詳細な計算方法は解読していない）。

バイナリを標準出力に吐き出すので、テキストモードがあるシステムでは動かない。作者は「今日にそんな環境を使っている人はいないので些細な問題 ;)」と言っているが、Windowsのことを指している。

読者の演習問題が3つある。

ポリゴンを増やしてスムーズにせよ（1点）
チューブの断面の半径を変える方法を示せ（2点）
有名な超越数の定数、πの位置を示せ（5点）

問1の答え：#define H 20の数字を変えれば良い。審査員が冒頭で./persano 3 16 > knot3-16.gifを例にあげているが、これはデフォルトのHでは小さすぎておかしくなる。

#define H 100とした結果がこちら。

問2の答え：sr = 15をsr = 5などとすれば細くなる。./persano 5 2の例。

問3の答え：pow(25/log(36),4)/377が2πの近似値になっている。