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)

Note that non-ascii characters in the above source code will be escaped (such as \x9f).

download

return to the top page