如何优化岩浆多项式可简化性代码？

Question

我有以下岩浆代码：

Z := IntegerRing();
K<x> := PolynomialRing(Z);

for n := 9 to nmax do 
  for d := 3 to Floor(n/2) do 
    if n mod d eq 0 then
      break;
    end if;
    bmax := Ceiling(amax^(d/n));
    for b := 2 to bmax do
      lower := Ceiling(b^((n-d)/d)-b^((1-d)/d));
      upper := Floor(b^((n-d)/d)+b^((1-d)/d));
      if lower le upper then
        for c := lower to upper do
          F := Factorization(x^n+x+b*c);
          if #F ge 2 then 
            print n, b*c, F;
          end if;
        end for;
      end if;
    end for;
  end for;
end for;

其中nmax和amax是设置变量。该脚本检查多项式因式分解。如何优化它？