# include<stdio.h>
# include<math.h>
/*****************************************************/
void iternewton(ax,ay)
 double *ax,*ay;
{ double 
 x  = *ax,
 y  = *ay,
 ex = exp(x),
 ey = exp(y),
 xx = x*x,
 yy = y*y,
 f  = x*xx + 3*x*y + y*yy,
 g  = ex - 2*x + ey - y - 3,
 fx = 3*(xx + y),
 fy = 3*(x + yy),
 gx = ex - 2,
 gy = ey - 1,
 detr = 1/(fx * gy - fy * gx),
 detx = f * gy - fy * g,
 dety = fx * g - f * gx,
 dx = detx * detr,
 dy = dety * detr;
 *ax  = x - dx;
 *ay  = y - dy;
}
/*****************************************************/
void j(ax,ay)
 double *ax,*ay;
{
 int i,im=8;
 for (i=0;i<im;i++){
  iternewton(ax,ay);
//  printf("%G %G\n",*ax,*ay);
 };
} /* j */
/*****************************************************/
int main()
{
 double x,y,k1x=0,k1y=0,k2x=0,k2y=0,k3x=0,k3y=0,k4x=0,k4y=0,
  xa=-20,xb=20.001,ya=-20,yb=20.1002,v1,v2,v3,v4,vmin;
 int r,in,inmax=100000,ix,ixmax=1000,iy,iymax=1000;
 x= 1.6;y=-0.8;j(&x,&y); k1x=x; k1y=y;
 x=-0.3;y= 1.0;j(&x,&y); k2x=x; k2y=y;
 x=-0.8;y=-0.2;j(&x,&y); k3x=x; k3y=y;
 x=-0.5;y=-1.1;j(&x,&y); k4x=x; k4y=y;

 for (iy=0;iy<iymax;iy++)
 for (ix=0;ix<ixmax;ix++){
  x = xa + (xb-xa)*ix/(double)ixmax;
  y = ya + (yb-ya)*iy/(double)iymax;
  in=0;
  do {
   iternewton(&x,&y);
   v1 = (x-k1x)*(x-k1x)+(y-k1y)*(y-k1y);
   v2 = (x-k2x)*(x-k2x)+(y-k2y)*(y-k2y);
   v3 = (x-k3x)*(x-k3x)+(y-k3y)*(y-k3y);
   v4 = (x-k4x)*(x-k4x)+(y-k4y)*(y-k4y);
   r = 0;
                vmin = v1, r=1;
   if (v2<vmin) vmin = v2, r=2;
   if (v3<vmin) vmin = v3, r=3;
   if (v4<vmin) vmin = v4, r=4;
  in++;
  } while (in<inmax && vmin>1E-6);
  printf("%d\n",r);
 };
 return(0);
} /* main */
/*****************************************************/