NP-hardness of p-adic linear regression

$p$-adic linear regression is the problem of finding coefficients $\beta$ that minimise $\sum_i |y_i - x_i^\top\beta|_p$. We prove that computing an optimal solution is NP-hard via a polynomial-time reduction from Max Cut using a regularisation gadget.