I want to compute gcdX((X−y)2−a,Xq−1−1) with respect to X(taking y as a field constant).
I can't see any direct implementation of this in sage.Can any one suggest how to implement it.
Here Arithmetic is over GF(p) and y is root of cyclotomic polynomial of degree r over GF(p) and q=pr