from polymod import PolyMod, PMod
from xgcd import xgcd

Poly7 = PMod(7)
a = Poly7('5x^2+3x+4')
p = Poly7('x^3+x+1')
print 'xgcd param, mod 7:', a, 'and', p
g, s, t = xgcd(a, p)
print 'gcd:', g, ' with coefficients:', s, 'and', t
print 'so gcd is not 1, but another nonzero constant.'