for(x in 3:2e3)if(all((x*x+x+x)%%(3:x-1)))cat(x,",",x+2,"
",sep="")