for(i=1, 100, a = round(zeta(1-2*i)*(6*i)!)/(6*i)!;print("zeta(1-", 2*i,") = ", a," ", factor(denominator(a))))