main=mapM print[x|x<-[2..6^5],all(`elem`"2357")$show x,gcd(3^x-3)x==x]