This approach is cool and equivalent to what I showed you, but it is more fragile in what it depends on, e.g.

```
>>> (x + z*(x + x*y)).expand()
x*y*z + x*z + x
>>> _.coeff(x)
y*z + z + 1
```

I would feel better if the more straightforward approach that I demonstrated is used.