Suppose we have independently distributed $X_i \sim \text{Uniform}(0,a_i)$ where the $a_i>0$ are fixed numbers. I want to obtain the probability that $X_j=X_{n-i+1,n}$ where $X_{n-i+1,n}$ is the $n-i+1$-th order statistic of the sample $X_1,...,X_n$.
Obviously, $\{X_i\}_{i=1}^n\overset{d}{=}\{a_i U_i\}_{i=1}^{n}$ where $U_i$ are i.i.d standard uniform random variables. Is there a smarter way to calculate this probability other than integrating (and summing) over sets of the form $\{a_iU_i> a_jU_j>... \}$?