main :: IO ()
main = mapM_ print [count 1 1 1 1, count 0.5 0.5 0.5 1]
count :: Double -> Double -> Double -> Double -> Int
count x y z r = length
  [(a, b, c)
    | a <- map fromIntegral [floor $ x - r .. floor $ x + r + 1]
    , b <- map fromIntegral [floor $ y - r .. floor $ y + r + 1]
    , c <- map fromIntegral [floor $ z - r .. floor $ z + r + 1]
    , inside a b c]
  where
    inside x0 y0 z0 = (x - x0) ^ 2 + (y - y0) ^ 2 + (z - z0) ^ 2 <= r ^ 2
 

bWFpbiA6OiBJTyAoKQptYWluID0gbWFwTV8gcHJpbnQgW2NvdW50IDEgMSAxIDEsIGNvdW50IDAuNSAwLjUgMC41IDFdCmNvdW50IDo6IERvdWJsZSAtPiBEb3VibGUgLT4gRG91YmxlIC0+IERvdWJsZSAtPiBJbnQKY291bnQgeCB5IHogciA9IGxlbmd0aAogIFsoYSwgYiwgYykKICAgIHwgYSA8LSBtYXAgZnJvbUludGVncmFsIFtmbG9vciAkIHggLSByIC4uIGZsb29yICQgeCArIHIgKyAxXQogICAgLCBiIDwtIG1hcCBmcm9tSW50ZWdyYWwgW2Zsb29yICQgeSAtIHIgLi4gZmxvb3IgJCB5ICsgciArIDFdCiAgICAsIGMgPC0gbWFwIGZyb21JbnRlZ3JhbCBbZmxvb3IgJCB6IC0gciAuLiBmbG9vciAkIHogKyByICsgMV0KICAgICwgaW5zaWRlIGEgYiBjXQogIHdoZXJlCiAgICBpbnNpZGUgeDAgeTAgejAgPSAoeCAtIHgwKSBeIDIgKyAoeSAtIHkwKSBeIDIgKyAoeiAtIHowKSBeIDIgPD0gciBeIDIK