until(print,for(i=!for(i=#r=[[0,input<q=""]],#F=input,for(j=#e=F[i][2],if(I-c=factor(I*p=F[i][1])[1,1],e,c=p;e%2&r=0;e\2),r=vecsort(vector(#r*2,k,x=r[k\/2];vecsort([abs(real(d=(x[1]+x[2]*I)*c=conj(c))),imag(d)])),1,8))),#r,print1(q,r[i][1]"^2 + "r[i][2]"^2");q=" = "))