static
void IrredCombine(ZZ_pEX& x, const ZZ_pEX& f, const ZZ_pEX& g)
{
if (deg(f) < deg(g)) {
IrredCombine(x, g, f);
return;
}
// deg(f) >= deg(g)...not necessary, but maybe a little more
// time & space efficient
long df = deg(f);
long dg = deg(g);
long m = df*dg;
vec_ZZ_pEX h(INIT_SIZE, dg);
long i;
for (i = 0; i < dg; i++) h[i].SetMaxLength(df);
h.SetLength(1);
set(h[0]);
vec_ZZ_pE a;
a.SetLength(2*m);
for (i = 0; i < 2*m; i++) {
a[i] = ConstTerm(h[0]);
if (i < 2*m-1)
MulByXPlusY(h, f, g);
}
MinPolySeq(x, a, m);
}
void DoMinPolyMod(zz_pX& h, const zz_pX& g, const zz_pXModulus& F, long m,
const vec_zz_p& R)
{
vec_zz_p x;
ProjectPowers(x, R, 2*m, g, F);
MinPolySeq(h, x, m);
}