This is more complicated than just nffactor above, a few hours in this
case and you'll probably get no improvement. It might be better to start
from T = polredbest(A) [ will be monic ], and iterate polcompositum(T,T)
and polredbest() to obtain a "simple" (monic) polynomial.