replace the iteration-based `mulmod()` with `math/bits.Mul64` + `math/bits.Div64`
replace the iteration-based
mulmod()withmath/bits.Mul64+math/bits.Div64