When 3712 is in you will get

```
>>> solve(x-1/sqrt(1+x**2))
[sqrt(-2 + 2*sqrt(5))/2, -sqrt(-2*sqrt(5) - 2)/2]
```

Perhaps you could review https://github.com/sympy/sympy/pull/1925 ?