Primer paso, deshacerse de varias raíces por determinar el mcd de a $p$ y su derivado $p'$ y si no es constante, dividiendo $p/gcd(p,p')$. Hay métodos para determinar un aproximado de mcd para no exacta de los coeficientes.
Ahora suponiendo que $p$ es cuadrado-libre, calcular el discriminante de $p$, es decir, una normalizado múltiples de la resultante de $p$$p'$. En la versión normalizada, haciendo caso omiso de los signos, uno tiene
$$
|disco(p)|=\prod_{i<j}|z_i-z_j|^2
$$
si $n=\deg p$ $z_1,...,z_n$ son las raíces de $p$. Para el siguiente, supongamos que $|z_1-z_2|$ es la distancia mínima. Todas las demás distancias será menor de lo $2R$ si $R$ es un límite en las magnitudes de las raíces (Cauchy, Lagrange, tal vez refinado a través de Graeffe iteraciones). Entonces
$$
|z_1-z_2|\ge r=\frac{\sqrt{disc(p)}}{(2R)^{(n+1)(n-2)/2}}
$$
En la no-exacto aritmética, si la precisión utilizado para el aproximado de mcd es alta, cerca de las raíces no pueden ser tratados de la misma raíz, por lo que la resultante discriminante y el resultado de la raíz de la separación de radio $r$ esto se reflejará en la pequeña distancia. Va para otro lado, permitiendo un amplio margen en el cómputo aproximado de la dpc resultará en una raíz de la separación de radio que se llamaría más bien de una raíz clúster de la separación de radio.