sympy求多项式的最大公因式

1.源码实现

from sympy import *
import pprint

x = symbols("x")

f = x**4 + x**3 + 2*x**2 + x + 1
g = x**3 + x**2 + x + 1

p = solve(g, x)
q = [ i for i in p if f.subs({x:i}) == 0 ]

gcd = expand('*'.join(['(x-'+str(i)+')' for i in q ]))

pprint.pprint(gcd)

2.运行及其结果

python3 example.py
x**2 + 1