// factorization of polynomials over finite fields
k<w> := GF(2, 100);
P<x> := PolynomialRing(k);
time Factorization(x^100 + x + 1);