Probability of one point being closer to fixed $x \in \mathbf{R}^2$ than another

Question

Let $x, x^{(0)}, x^{(1)} \in \mathbb{R}^2$, $r_1, r_2 \sim U[-1, 1]$. Point $x$ is fixed, $x = (0.32, 0)$, and $x^{(0)} = (0, r_1)$, $x^{(1)} = (1,r_2)$. What is the probability that the fixed point $x = (0.32, 0)$ will be closer to $x^{(1)}$ than to $x^{(0)}$, by euclidean distance?

Attempt

$x_0 \sim U[-1, 1]$, $x_1 \sim [0, 2]$, let $x_0' = (x_0 - x) \sim [-1.32, 0.68], x'_1 = (x_1 - x) \sim [-0.32, 1.68]$. Thus, we need to find $P(|x'_0| < |x'_1|)$. Any suggestions? I believe it could be dumbed down to a simple double integral.

Answer 1

The goal is to find $P(|x_0^\prime| - |x_1^\prime| \leq 0)$, where the distance is understood to be euclidian.

\begin{align} P(|x_0^\prime| - |x_1^\prime| \leq 0) & = P(((-.32)^2 +r_1^2)^{\frac{1}{2}} - ((.68)^2 +r_1^2)^{\frac{1}{2}} \leq 0)\\ &= P((.32^2 +r_1^2)^{\frac{1}{2}} \leq (.68^2 +r_1^2)^{\frac{1}{2}})\\ &= P(r_1^2 - r_2^2 \leq .68^2 - .32^2) \end{align}

The easiest way that I can think of answering this is by considering this the integration of the region $1\{r_1^2 - r_2^2 \leq .68^2 - .32^2\}$ over the joint density of $r_1^2$ and $r_2^2$. It is sufficient to find the density of the marginals since $r_1$ and $r_2$ are independent. For convenience I will write $X = r_1^2$ and $Y = r_2^2$

By the transformation method, we can find the CDF of $X$ in terms of the CDF of $r_1^2$, $P(r_1 \leq z) = \frac{z+1}{2}1\{z \in [-1,1]\}$. Analgously we can do the same for $Y$ and $r_2^2$, so we skip this step

Since $r_1$ has support $[-1,1]$, $X$ has support on $[0,1]$. Thus for $x \in [0,1]$

\begin{align} P(X \leq x) &= P(r_1^2 \leq x)\\ &= P(\{r_1 \leq x^{\frac{1}{2}}\} \cup \{-r_1 \leq x^{\frac{1}{2}}\})\\ &= P(r_1 \leq x^{\frac{1}{2}}) + 1 - P(r_1 \leq -x^{\frac{1}{2}})\\ &= \frac{x^{\frac{1}{2}} +1}{2} + 1 - \frac{1 -x^{\frac{1}{2}}}{2} \end{align}

Since the pdf is the derivative of the CDF with respect to x, the pdf of $Y$ is analagous, and $X \perp \!\!\! \perp Y$ we get

\begin{align} p_X(x) &= \frac{1}{2}x^{\frac{-1}{2}}1\{x \in [0,1]\}\\ p_{X,Y}(x,y) &= \frac{1}{4}x^{\frac{-1}{2}}y^{\frac{-1}{2}}1\{x \in [0,1]\}1\{y \in [0,1]\} \end{align}

Now using this joint density, we can express the original probability as an integral

\begin{align} P(r_1^2 - r_2^2 \leq .68^2 - .32^2) &= P(X -Y \leq .68^2 - .32^2)\\ &= \int\int 1\{X -Y \leq .68^2 - .32^2\}p_{X,Y}(x,y)dxdy\\ &= \int_0^{(.68^2 - .32^2)}\int_0^1 \frac{1}{4}x^{\frac{-1}{2}}y^{\frac{-1}{2}}dxdy + \int_{(.68^2 - .32^2)}^1\int_{y - (.68^2 - .32^2)}^1 \frac{1}{4}x^{\frac{-1}{2}}y^{\frac{-1}{2}}dxdy\\ &= .6 + .19775\\ &= .79775 \end{align}

To check and make sure I didn't make any obvious errors, here in a numerical simulation in R.

set.seed(682322)
a <- -1
b <- 1
n <- 1000000

r1 <- runif(n,a,b)
r2 <- runif(n,a,b)
c1 <- sqrt(.32^2 + r1^2)
c2 <- sqrt(.68^2 + r2^2)
mean(ifelse(c1 -c2 < 0,1,0)) 

0.797967