for(q,3,2e3,x:=q*(q+i:=2) while(x%i>0,i=i+1)i