main=mapM putStrLn[show a++',':show b|a:b:_<-1#[2..2000],a+2==b]
a#(p:x)=[a,p]:p#[q|q<-x,mod q p>0]