R=range S=set(R(2000)) for i in R(3,2000): S-=set(R(i+i,2000,i)) if(i in S)*(i+2 in S):print"%d,%d"%(i,i+2)