As we now have Singular 4, I just made a little random test:
set_random_seed(0)
for i in range(1000):
P = 1
while random() < 0.9:
M = (randint(1, 5) * ([a, b, c, d, e][randint(0, 4)])^randint(1, 5))
A = (randint(1, 5) * ([a, b, c, d, e][randint(0, 4)])^randint(1, 5))
P = P*M + A
Q = 1
while random() < 0.9:
M = (randint(1, 5) * ([a, b, c, d, e][randint(0, 4)])^randint(1, 5))
A = (randint(1, 5) * ([a, b, c, d, e][randint(0, 4)])^randint(1, 5))
Q = Q*M + A
if i == 112:
print(P)
print()
print(Q)
F = factor(P*Q)
and this got stuck. Might be that it just needs more memory than I have; I
don't know how much it normally uses.
--
Jori Mäntysalo