One more question about if something will be welcomed. Tonelli-Shanks algorithm?

In my vision it belongs to Boost.Multiprecision because it requires for the modular adaptor and Jacobi symbol computation to be used, which is impossible to do from Boost.Integer without introducing the dependency on the Boost.Multiprecision.