#define/**/Alan/**/(fflush(0),j=c=0;++c<b[6]+7;v=b[c]^(b[c+ char*H="close\r\nContent-Type:text/html\r\n\r\n<canvas id=\ 'c'width='128'height='128'style='width:256px;height:256px'\ ><script>x=y=-1;v=window.c;v.onmousemove=function(e){x=e.p\ ageX-v.offsetLeft;y=e. pageY-v.offsetTop};w=\ new WebSocket('ws'+ location.href.subs\ tr(4));w.binary" "Type='arraybuff\ er';w.onmess" "age=function(\ e){c=v.getC" "ontext('2d')\ ;b=c.crea" "teI" "mageData(12\ 8,128);u= new Uint8Array(d=e. data);b.dat\ a.set(u) ;c.putImageData(b,0,0) ;u[0]=x/2;\ u[1]=y/2 ;w.send(d.slice(0,x<0?0: 2));x=y=-\ 1}</scr" "ipt>",b[9999];float*e,d[ 65536],u,v ,lu,lv, z;void*f;typedef unsigned long long l;l*p,t[ 99]={0x67452301,0xEFCDAB89, 0x98BADCFE ,308438* 881,3285377520 },i,j,k ,n,m[204]; #include /*io*/ <stdio.h> /*turing*/ void s(){ ; for( k =0;k++<7 * 2 *9;)for(i= 0;++i<127; ) e=d + i* 4 +512*k,1[e] +=(i-u)*(i -u)+(k-v)* (k- v )<20?*e+=UV_DROP:0;}l*q= m, g=1LL<<32; l(h)(l v,l (a)){ return(v<<a)|(v%g>>(32-a ) );}int main (int c/*ha r*/, char**y){for(e=d;fgets(b ,m[97]=480,\ stdin)&&*b >31; )for( i=j=0;n=b[i],i<18 ? (l)"sec-webs\ ocket-key: "[i] ==(n|32):n-32? m[j/4] |=(j< 24?n:(l) "258EAFA5-" "E914-47DA-95CA-C5AB0DC85B" "11\x80"[j-24]) <<(3-j%4)*8 ,++j<61:1;)i++;for(j=3;--j; ){for(i=4;++i<74 ;)i>9?p=q+i- 2 ,p[8]=h(p[5 ]^*p^p[-6]^p[-8] #define/*hah aha <--mouth*/ laplacian(u)l##u ,1):(t[i]=t[9 -i] );for(p=t+7 ,i=-1;++i<82;k =*p++,p[2]=h(p[1],5) +p[-3]+n+ (i/+ 20%2?*p^k^p[- 2]:(*p&k)|((i< 20?~*p:*p|k)&p [-2 ]))+*q++,t[+ 81-i]+=(i>76)* (* p= h(* p,30*(i<80) )))for(;n-k;n=((( 1LL +(9<<i/20 )/8)<<60)/k+k)/2) k=n ;} for (i =0;i<57; i++)i<30?q[i%5* 16/ 3+ i/5] |=t[i%5]%g<<4 >>(30+i%5%3*2-i/5 *6)& 63,0 :(n =*q++,b[i-30]= n<62?n-37+"fl!" [n /26]:n *4- 205);for(prin\ tf("HTTP/1.1 "/* | */"%d " "OK\r\nConnecti\ on:%s%.27s=\r\n\r\n" ,101 +99*! *m,!*m?H:"upgrade\ \r\nUpgrade:webso" "c" "ket\r\nSec-WebSocket-Acce\ pt:",b);e-d<65536 ; e+= 4)e[1]=(*e=UV_BACKGROUND);for (;*m;c=0){for(f= c< 2 ?0:fopen(y[1],"r");f&&fscanf(f, "%f,%f,",&u,&v)>0;)s ();for(;j<SPEED*4;j+=2)for(i=0;i< 63504;e+=5-j%4*2,z= TIMESTEP,*e=u+z*(elta_u),e[1]=v+z* (elta_v),i+=4)for(c= 0;c+2;u=v,v=e[4],lu=lv,lv=e[-508]+*e+ e[8]+e[516]-v*4)e=d+i+512+i/504*8+j%4-c--;for(;n=c%4,65544> c;c++)putchar(c>7?n>2||(u=d[c-n-7]*4*(n<2?!n?RGB))>255?255: #define/*x*/Turing(x)4]=getchar()&127))u=v;if(b[6]>1)s();}} u<0?0:u:c<0?124-c*3:c==5);/* -- The word "genius" is */ for Alan Turing (1912--1954)

引用元：https://www.ioccc.org/2015/endoh1/prog.c

審査員・作者による説明：https://www.ioccc.org/2015/endoh1/hint.html

動作

反応拡散系のシミュレータ。

コンパイル方法が難しい。 たとえば次のようにビルドする。

$ gcc -o prog prog.c \
    -Delta_u="Du*laplacian(u)-u*v*v+F*(1-u)" \
    -Delta_v="Dv*laplacian(v)+u*v*v-v*(F+K)" \
    -DF=0.040 -DK=0.060 -DDu=0.200 -DDv=0.100 \
    -DUV_BACKGROUND=1,0 -DUV_DROP=-0.5,0.5 \
    -DTIMESTEP=1.0 -DSPEED=2 -DRGB=255:128:192

そして、ucspi-tcpのtcpserverコマンドを使って次のように起動する。

$ tcpserver -v 127.0.0.1 10333 ./prog ioccc.txt

ブラウザでhttp://localhost:10333を開くと、IOCCCと書かれた画面がじわっと変化していく。

マウスカーソルをブラウザの画面の上で動かすことでかき混ぜることができる。

他にもコンパイルオプションでいろいろな挙動をさせられる。

詳しくはhint.textを参照のこと。

解説

反応拡散系は、何種類かの化学物質がお互いに反応しつつ、全体に拡散していくという様子を表現する数理モデル。 身近なところでは、豹やシマウマなどの柄は反応拡散系でモデル化できるとされている。

このプログラムは、128 x 128のマスのあるフィールドで、2種類の化学物質の量をシミュレーションする。 各マスの中で物質が化学反応したり、隣接するマスから流入または流出することで、単位時間ごとに物質の量が増減していく。

ブラウザをUIとして利用する点が新しい。 WebSocketとcanvas要素を使っている。 このプログラムは、サーバとしてブラウザからのリクエストを待ち受け、JavaScriptやCanvas要素を含んだhtmlを返す。 ブラウザはそれを見て接続をWebSocketにアップグレードすることを要求するのでそれに応える。 WebSocketに反応拡散系のシミュレーション結果を流し、ブラウザ側はそれをCanvasに表示し続ける。 また、ブラウザ側のマウスイベントはWebSocket経由でプログラムに通知し、シミュレーションに影響を与える。

TCPの待受自体はtcpserverにまかせるが、HTTPリクエストの処理、WebSocketのアップグレード対応は自力で行っている。 WebSocketにはSHA1計算やBase64エンコーディングが必要で、これらもprog.cの中で実装されている。

コードの形状はアラン・チューリングの肖像画。 チューリングと言えばエニグマの暗号解読、チューリング機械、チューリングテストなど多数の業績で知られるが、反応拡散系の提唱もチューリングの業績の1つ。 コードには絵画らしく、最下行にAlan Turing (1912--1954)と書いてある。 一方で、コメントでいろいろ落書きをしている。 口のあたりにhahahaや<--mouthなど。

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