f n=2^n-sum[f d|d<-[1..n-1],mod n d<1]
main=mapM(print.(div=<<f))[1..50]