from sympy import *
left, right, pln, a = symbols('left right pln a')
left = a
right = -2 * a + 1
pln = left - right
print(pln) # 3*a - 1
print(pln / polys.polytools.LC(pln)) # a - 1/3
print(solvers.solve(pln, a)) # [1/3]
он должен быть стабильным и иметь самые последние программы