Factor[a^2 b + a b^2 + a + b - a b - 1]