x := 13*17*31;
y := 13*19*29;
for j from 1 to 200 do
   print(gcd(x^(300 + irem(j,181)), y^(200 + irem(j,183))))
end_for;

purge(x);
purge(y);