tinypt_nee

#include"math.h" // sqrt
#include"stdio.h" // fopen
float a,b,c,x,y,z,u,v,w,i,j,d,e,s,t,X,Y,Z,A[1966080]={0};unsigned l,N=0,r=9,p,q
,n=1966080;float U(){r+=r*r|5;return r*4.7e-10f-1;}void I(){float h[]={a,b,c,x,
y,z};q=4,t=n;d=a*x+b*y+c*z;e=d*x-a,i=d*y-b,j=d*z-c;j=e*e+i*i+j*j;d=-sqrt(1-j)-d
;if(j<1&d<t&d>0)t=d,i=3;if(t>24)for(;q;){j=(1-h[--q])/h[q+3];if(j<t&j>0)i=q,t=j
;}}void P(){N++;for(p=0;p<n;){a=b=0,c=-6,x=p/3%1024/90.-7,y=p/3072/80.-4,z=9;j=
sqrt(x*x+y*y+81);x/=j,y/=j,z/=j;l=0,X=Y=Z=1.7;for(;I(),l<9;l++)if(t>24)l=9;else
{a+=t*x,b+=t*y,c+=t*z,u=a,v=b,w=c;if(i<3){u=v=w=0;if(i<2){Y/=4,Z/=4;if(!i)u=-1;
else v=-1,X/=4;}else w=-1;}j=(u*x+v*y+w*z)*2,x-=u*j+U()/6,y-=v*j+U()/6,z-=w*j+U
()/6;if(i<3|U()<0.4){if(i>2)X/=3,Y/=3;x=U()*5-2-a,y=-8-b,z=U()*3-3-c;s=sqrt(x*x
+y*y+z*z);x/=s,y/=s,z/=s;I();d=u*x+v*y+w*z;if(t>24&d>0)d/=s*s/80,A[p]+=d*X,A[p+
1]+=d*Y,A[p+2]+=d*Z;x=U(),y=U(),z=U();s=sqrt(x*x+y*y+z*z);if((d=x*u+y*v+z*w)<0)
s=-s;x/=s,y/=s,z/=s;d/=s/2;X*=d,Y*=d,Z*=d;}}p+=3;}}int main(){FILE*g=fopen("o",
"wb");int h[]={2<<16,2<<28,0,41944064,24},i=0;fwrite(h,2,9,g);for(;i<999;i++)P(
);for(p=n;p;){i=A[--p]/N*255;i=i>255?-1:i;fwrite(&i,1,1,g);}} // JB, Jan 8 2022