Floyd-Steinberg Dithering

Given a 64x64 grayscale image in plain PGM format, convert it to a palette of black (0) and white (255) using the Floyd-Steinberg algorithm. Print the result, using "#" for black and spaces for white.

The input is 3 lines of header information, followed by 4096 pixel values, each on their own line.

    P2          # magic number
    64 64       # image dimensions
    255         # maximum pixel value
    19          # (0,0) value
    18          # (1,0) value
    ...

• Let a be the pixel value.
• Let b be the nearest palette value to a (either 0 or 255).
• Set the current pixel to b. Now a "debt" of a-b needs to be spread to 4 of the neighboring pixels. For each of these 4 neighbors and its corresponding fraction f (as seen below), add floor((a-b)*f) to its value.

     --   --   --
     --   []  7/16
    3/16 5/16 1/16

• Some pixel values may temporarily become less than 0 or greater than 255.
• When processing pixels on the edge, any adjustments made to pixels outside of the image should be ignored.

External links: Wikipedia