$ρ$-regularization subproblems: Strong duality and an eigensolver-based algorithm