for (I=1; I<1000; I++) {
   print(fctr(x^50-y^50));
}
end$