// factor polynomial over integer ring
P<x> := PolynomialRing(IntegerRing());
time Factorization((2*x+1)^80-1);