main=mapM(\n->print$sum$map(gcd$2^n)[1..n])[1..500]