import sympy
sympy.sqrt(3)

x, y = symbols('x y')
expr = x + 2*y
expr

from sympy import expand, factor
expanded_expr = expand(x*expr)
expanded_expr

factor(expanded_expr)