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);