!7

import java.awt.geom.*;

p {
  BufferedImage img = loadBufferedImage(#1006973);
  BWImage bw = getAlphaChannel(img);
  Image2B alpha = new Image2B(bw);
  
  int w = alpha.w(), h = alpha.h();
  
  Path2D.Float path = new Path2D.Float();
  
  int x = w/2;
  int y = -1;
  while (y < h && !isInteresting(getQuadrant(alpha, x, y))) ++y;
  
  int x1 = x, y1 = y;
  path.moveTo(x, y);
  
  new HashSet<Pt> seen;
  while (seen.add(new Pt(x, y))) {
    int q = getQuadrant(alpha, x, y);
    
    Pt next = null;
    for (Pt p : neighborhood(x, y)) 
      if (!seen.contains(p) && isInteresting(getQuadrant(alpha, p.x, p.y))) {
      next = p;
      break;
    }
    
    if (next == null) break;
      
    x = next.x;
    y = next.y;
    path.lineTo(x, y);
  }
  
  print("Have " + n(seen, "points"));

  showImage(alpha);
  showImage(renderShape(path, w, h));
}

sbool isInteresting(int q) {
  ret q != 0 && q != 0x1111;
}

static int getQuadrant(Image2B img, int x, int y) {
  ret (getQuadrant_getPixel(img, x, y) ? 0x1000 : 0) +
    (getQuadrant_getPixel(img, x+1, y) ? 0x0100 : 0) +
    (getQuadrant_getPixel(img, x, y+1) ? 0x0010 : 0) +
    (getQuadrant_getPixel(img, x+1, y+1) ? 0x0001 : 0);
}

static bool getQuadrant_getPixel(Image2B img, int x, int y) {
  ret x >= 0 && x < img.w && y >= 0 && y < img.h && img.getPixel(x, y);
}

static L<Pt> neighborhood(int x, int y) {
  ret ll(new Pt(x-1, y-1), new Pt(x, y-1), new Pt(x+1, y-1),
    new Pt(x-1, y), new Pt(x+1, y),
    new Pt(x-1, y+1), new Pt(x, y+1), new Pt(x+1, y+1));
}