isqrt
ed = dx + dy == 0 ? 1: isqrt(dx * dx + dy * dy);
return (isqrt(dp / di));
ed = dx + dy == 0 ? 1 : isqrt(dx * dx + dy * dy);