sympy求多项式的最大公因数和最小公倍数

1.程序源码

from sympy import *

x = symbols('x')

a = x**3 + 1
b = x**2 - 1

#最大公约数
e = gcd(a, b)
#最小公倍数
f = lcm(a, b)

print(e)
print(f)

2.运行及结果

$ python3 example.py
x + 1
x**4 - x**3 + x - 1