read a b 'y?y=x%(x=y),_:x'
for((x=a,y=b,a/=_,b/=x,l=a*a+b*b;c<x*x*2;++c)){
for((i=0;i<l;++i)){
echo $[(a*(i+c/x/x)%l+c%x*l)*x] $[(b*i%l+c/x%x*l)*x]
}
}