#include <stdlib.h> #include <stdio.h> #define a typedef a long N; #define i =m(p(r,o,v),e,d N G,l,I,C,B,A,W,L,S , R,O,c,k,s = 80, U=13 ,T=169; a double P; #define F for( a P*E; a char*w; P sqrt(P); N H(){ F; O=scanf(" %[#]%*[^\n]",&O); ); scanf("%ld%*c",&O); #define D return D O; } #define V =malloc(sizeof( #define M M( #define X += P M N R){ s= s>=k?printf( "%ld%c" +3*(C<1),s-k, R&7?32:10):s; D R>0?B /2?H():getchar ():R?.9+.7*M R+1):0; } P t(N x,P K){ D x?(x>U?0:t(x+2,K)*(x-1) /x/K)+1/sqrt(K): t(2,K*K/U+1)*K; } N _(N J,N x){ F*(x?&W:& A)=W+J+x; (O=x=W/(J=k*k)) ||A<J||(( O=2*W/J) &&2*A/J<3); A X A-J+1){ F J*=x+O; S--*B&&x==O ; O=x=!J){ s X x+s; M 0); } W X W-J; S X 2; L X L-J*k+M S%8==B); } D 1; } E p(E J,E x,E O){ F R=T; R--; )O[R]=J[R]+.8*x[R]; D O; } P m(E W,E x,P s){ F R = U; --R; ){ E A=W+R*U; W[R]X s/ 12; *W X s; *A-=*x++; F; A>W; ){ *A/=*W; F O=R; O; O--)A[O]-= W[O]**A; A-=U; } W X 14; } D*W; } int main(int z,w*y){ N K=s,g; Y; M z); s= B=C=15-H(); s=G= M 8); s=l=M C); G*=g=C%3?1:3; if(I=s=M 8)){ N b=G+9,x=b*340+U,J=b*(l+5), * n V N ) * J ),*j=n+5*b; E h V P)*x),q=h+b,e=q+ b; P d=K; F; J--; )n[J]=J<b ? 0: I / 2 ; F; x; )h[--x]=0; s=M x); k=256; C =4-C; B=C/4; s=_(_(T,0),0); F ; -- z ; ){ R=atoi(*++y); c=R>0?K=R,c :-2*R; } z=G/K+1; z X g>z; I++; c++; F ; ++ J<l; ){ P u=0,C,y; E H =e+U,o= H+b*T,r =o-2*T,v=r-T; F; x<G; x++){ p( H , o -T,o); o X T; H X T; q[x]= .7*(q[x-1]+h [x] ); } p(H, H,r ); F; x; ) { E S= e ; w l = " !{ ,;lf6D@" ; j X 1-g; F; *++l; ){ *S++ =*j; j X*l%3*g- g+*l%5*b-3*b; } y =M x-G); y =M-J)* M-x)+ M-J-1)*y+.01; y= sqrt((u+q[--x]+.1)/y ); o-=T; C i)+ #define Z(x)(t(0,(x-C)/y)-l+1e-6*(x-Q)) .5; K=M-B); { N f=I,Q=0; F; f-Q>c; ){ P l=0 ; N H = C+c/2+1; H%=c; H X(f+Q+c-2*H)/2/c*c; l= Z(Q); O=(A-W)*Z(H)/Z(f); _(O,R=B?K>= H:L/k>O+W); *( R?&Q:&f)=H; } Q X f; *S=*j=Q/=2; s=B?s:Q+k; f=Q+n[x/z]; n[x/z]=x%z+J%(z*2/g)?f:!putc(x?l[4*g*f/z/z/I-8]:10,stderr); H-=U; F O=156; O--; ){ H--; *H=e[O%U]*e[O/U]/y+.8**H; } C-=Q+.5 ; y i/.9)-Q-C; p(p(r+T,r,r),H,r); h[x]=.7*h[x]+C*C; u X h[x]; u*=.7 ; d X C>0?-y:y; } } j X 9; } _(_(x,x),x); } D 0; }

引用元：https://www.ioccc.org/2000/bmeyer.c

審査員・作者による説明：https://github.com/ioccc-src/winner/blob/main/2000/bmeyer.hint

動作

グレイスケール画像の圧縮・展開ツール。

圧縮の例。端末に出てくるアスキーアートは飾り。

$ gcc -DY -o bmeyer bmeyer.c -lm

$ ./bmeyer < lenna.pgm > lenna.glic

 fffff666l;llllllllfffffffffffffffffffffflflllllffffffffff6DDflllllllllll
 ffff666fl;lllllllllfffffffffffffflfffflllllllllff666fffffffDD6llllllfl;,
 f66666ffl;llllllllllllfffffflllllllfffllllllllllffffffffffffDD6lllll,,,,
 666ff66fl;;llllllllllfflfffffff6666fffllllllllllfffffffffffff6DDll,,,,,,
 6fllf66fl;;lllllllllllllflfffff666666D6fllllllllffffffffffffff6Dl,,,,,,,
 fl;lf66fl;;llllllllllllllllfffff66666DDDD6flllllfffflfffffffffl,,,,,,,,,
 l;;lf66fl;;llllllllllllllllfffff6666DDDDDDD6l;;lffff;;lffffff;,,,,,,,;lf
 ;;;lf66fl;;lllllllfllllllfffffff66DDDDDDDDDDD6;lffff;,,lffff;,,,,,,,;fff
 ;;;lf66fl;;;llllll6llllllffflf666666DDDDDDDDDD6lffff;,l6DDDl,,,,,,,lff6f
 llllf666l;;llllll6fllllllllfff666666666DD66DDDDD6ffff6DDDDDf,,,,,,ffffff
 llllf666l;;llllllDfllllllfffff66666666666666666666DDDDDD6DD;,,,,;fffffff
 l;;lf666l;;llllllD6flllfffffffff666666666666666DDDDDDDD6f6f,,,,;ffffffff
 llllf666f;;llllllD6fflllffffffflllfllllllf666DD6DDDDD6flff,,,,;fffffffff
 llllf666f;;;llll;6Dffllfllfffll;;;;;;;;;f666666D66fflfff;,,,,;ff666fffff
 llllf666f;;;llllll66fllllll;;;,;,;,;;;f666666666;llf6f;,,,,,;ff66666ffff
 llllf666f;;;llllllfDlllll;;,,;,;;,,;lf666666DDDDf;;fl,,,,,,,lff666666fff
 llllf666f;;llllllllflll;;,,,,,,;;;f66f666666DDDD6;,ll,,,,,,lf6ffffffffff
 llllf666f;;llllllll;fl;,;,,,,,;,lf666fffff66DDD66l,,f,,,,,;fffff6fffffff
 ll;;f666f;;lll;ffff;;;;;,,,,,,;f6D6l;;;l;lf6D6l;;;,,f;;,,,ffffffffffffff
 ll;;f666f;;ll;lff;;;,;;l,,,,,;f66fllllf6flf6Dfll;,,,f;,,,lffffffffffffff
 ;;;;f666f;;llfl;,;l;,;ll;,,,;f6llfff6666ffffD6ffl;,,fl,,;fffffffffff6666
 ;;;,f666f;lll;l;;l;,,,;,l;,;f6;;lfff6666fflfD6ffl;,,ll,,lffffffffff6DDDD
 ;;,,f666f;lll;;;;ll,,,;l;,lff,,;llffff6ffflf66ffl,,,;f,;ffffffffff6DDDDD
 ;;;;f666f;;l;;;;;ll;;;,;llf;,,;;llffffffffff66fl,,;,,f;ffffffffff6DDDDDD
 ;;;;f666f;ll;;,,;;;ll;f6ll,,,,,;llfffffffff66fl,,,;;,llffffffffffDDDDDDD
 ,,,lf666flll;;;;l;;llf6fl,,,,;,;;llffffllffffl,,,,;;,lfffffffffffDDDDDDD
 ,,;lf666f;l;;,;,llllllf;,,,,,,,,;;llffffffffl,,,,,;;;lfffffffffl6DDDDDDD
 ,,,lf666fll,;,,,l;;llll,,;,,,,,,,;;;llfff66fl;,,,;;;;lflffffffflDDDDD6fl
 l;;,l666fl;,,;;;;lllllll;;,,,,,,,,lllfffff66DD6f;,;;;lffflllllllDDDDf;;;
 lff,l666ff,,,,;;,;;;;lfll;,,,,,,;,lllfffff6666DDD6;,,fffffffllffDD6l;;ll
 ,f6;l666fl,,,,;,,,,;l;fl;l,;,,,;l;lffffffff6666DDDD;,ffffffflf6DDDl;ll;l
 ,l6;l666f;,,,,;,,,,,;lfll;;;,,,ll;lfffffffff666DDDD6;ffffffllDDDDflll;ll
 ,l6ll666fl,,,,,,,,,,;;l;l,,,,,lll;ffffffffff6666DDDDfllfffffff6D6ll;;l;;
 ,;6f;666fl,,,,,,,,,,;l;;,,,,,;ll;lfffffffffff666DDDDD;;;;;;l6fD6;;;;;;l;
 ;l6flf66fl,,,,,;,,,,;;l;,,,,;;;;lfffffffffffff666DDDDflll;;;D6f;;;;;;ll;
 ,;66lf66f,,,,,,;,l,;,;l;,,,;llllfffffffffffffff666DDDD;;lllllll;;;;;ll;,

得られたファイルはpngより小さい。

$ wc -c lenna.pgm lenna.glic lenna.png
262159 lenna.pgm
127781 lenna.glic
150677 lenna.png

展開の例。おなじアスキーアートがでてくるけど省略。もとに戻ったことがわかる。

$ ./bmeyer < lenna.glic > lenna2.pgm
...

$ diff -s lenna.pgm lenna2.pgm
Files lenna.pgm and lenna2.pgm are identical

解説

標準エラー出力に出てくるアスキーアートは圧縮・展開の進捗を表示するためのもの。バー表示やパーセント表示より面白いでしょう、とのこと。

浮動小数点数演算を活用するので、FPUによってはうまく動かないことがあるみたい。いわゆる80 bit浮動小数点数演算問題。

コードは何らかのアスキーアートになっているが、意味はわからなかった。

圧縮アルゴリズムは独自のもの。 周囲の既知のピクセルから予測値を作り、それとの差分を載せるような感じだと思う。 作者が知る限り、どの画像圧縮よりも小さいデータにできたとのこと。