zeta = 1+2i: monomials of h = [(3, 2, 0), (0, 3, 2), (4, 0, 1)], deg h = 5, deg f = 10, forms with [3, 3, 3] terms
expected from BDJ's recursion (Prop. 2.7): e_1..e_6 = [10, 66, 454, 3114, 21368, 146488]   (d_j = [5, 8, 22, 48, 117, 205])
random line x_i(t) = p_i + q_i t modulo 2^31 - 1, seed 1
  n = 1: raw degree 10, common factor of degree 0, deg f^1 = 10  = e_1   (0.0 s)
  n = 2: raw degree 100, common factor of degree 34, deg f^2 = 66  = e_2   (0.0 s)
  n = 3: raw degree 660, common factor of degree 206, deg f^3 = 454  = e_3   (0.0 s)
  n = 4: raw degree 4540, common factor of degree 1426, deg f^4 = 3114  = e_4   (0.2 s)
  n = 5: raw degree 31140, common factor of degree 9772, deg f^5 = 21368  = e_5   (3.2 s)
  n = 6: raw degree 213680, common factor of degree 67192, deg f^6 = 146488  = e_6   (114.8 s)
[PASS] deg f^n on this line for n = 1..6: [10, 66, 454, 3114, 21368, 146488]; BDJ recursion: [10, 66, 454, 3114, 21368, 146488]
