I want to solve a system of three equations to get the values of three variables. It doesn't matter whether I use Mathematica 13.1 or the Wolfram Engine 13.1, the results are the same. Here is the code:

Solve[{(a + 1)*(b + 1)*(c + 1) == 3 &&
           (a + 2)*(b + 2)*(c + 2) == -2 &&
           (a + 3)*(b + 3)*(c + 3) == -1},
           {a, b, c}, Reals] // N

I get six solutions:

{{a -> -3.1451, b -> -1.52398, c -> 1.66908}, {a -> -3.1451, 
  b -> 1.66908, c -> -1.52398}, {a -> -1.52398, b -> -3.1451, 
  c -> 1.66908}, {a -> -1.52398, b -> 1.66908, 
  c -> -3.1451}, {a -> 1.66908, b -> -3.1451, 
  c -> -1.52398}, {a -> 1.66908, b -> -1.52398, c -> -3.1451}}

Unfortunately, none of them is correct. And why do I get six solutions?

The correct answer is: a = 1.669079088, b = -1.523976396, c = -3.145102692

Do you have any suggestion how to improve the results?