I want to solve the following optimization problem in variables $\theta_1, \theta_2, \dots, \theta_K$

\begin{equation}
    \begin{aligned}
    & \underset{\theta}{\text{minimize}}
    & & \sum_{k=1}^{K} \|g_k - \theta_k\|_2^2 \\
    & \text{subject to}
    & & \| \theta_m - \theta_l \|_2^2 \leq \| x_m - x_l \|_2^2 , \quad 1 \leq m < l \leq K,
    \end{aligned}
\end{equation}

where $x_i$ and $g_i$ are constants and $\| \cdot \|_2$ denotes the Euclidean norm.

Is there any closed-form for arbitrary $K$? Or an appropriate first-order iterative method which yields an approximate solution?