tl; dr: 當 $\color{black}{p=2}$ 時直接特判。當 $\color{black}{p\geq3}$ 時，我們一樣有 $\color{black}{x = \frac{-b+\sqrt{b^2-4ac}}{2a}}$，其中開平方根的部分可以用 Cipolla 或 Tonelli–Shanks 演算法，時間複雜度分別為 $\color{black}{O(\log p)}$ 與 $\color{black}{O((\log p)^2)}$。