给 u32 求逆的,可能是牛迭
constexpr u32 inv(u32 x){ u32 iv=2-x; for(int i=0;i<4;++i){ iv*=2-x*iv; } return iv; }