pub fn div255(x: u32) -> u32
round(x / 255) for x in 0..=65025: ((x + 128) * 257) >> 16 (pmulhuw), equal to (x + 128 + ((x + 128) >> 8)) >> 8