Finding a billion factorials in 60 ms with SIMD