from simpy import *

a,b=symbols('a b')
factor((-a**2+a)-(-b**2-b))