https://pt13762104.github.io/cp-library

If it exists, it can be done using sqrt-decomposition.