main=mapM print[n|n<-s,7==sum[p-gcd(p^2)n+1|p<-s,1>mod n p]]
s=[2..762]