Let f(x) = x^101 + x^94 + x^57 + x^33 - 1 and let g(x) = (x-1) * f(x) = x^102 - x^101 + x^95 - x^94 + x^58 - x^57 + x^34 - x^33 - x + 1
Start with the assumption that x^2+x+1 is a factor of f(x). If that is true then x^3-1 will be a factor of g(x).
Rearrange g(x) into (x^102-x^57) - (x^101-x^95) - (x^94-x^58) + (x^34-x) - (x^33-1). Each of the five binomials in this arrangement has x^3-1 as a factor, therefore g(x) has x^3-1 as a factor, which then implies that x^2+x+1 is a factor of f(x).
Another approach
