Revision history [back]

For your problem, since zeta(2n) has an explicit expression in terms of Bernoulli number I would rather go with

sage: [((-1)**(1+n//2)*2**(n-1) * bernoulli(n)/factorial(n)).factor() for n in range(2,20,2)]
[2^-1 * 3^-1,
 2^-1 * 3^-2 * 5^-1,
 3^-3 * 5^-1 * 7^-1,
 2^-1 * 3^-3 * 5^-2 * 7^-1,
 3^-5 * 5^-1 * 7^-1 * 11^-1,
 3^-6 * 5^-3 * 7^-2 * 11^-1 * 13^-1 * 691,
 2 * 3^-6 * 5^-2 * 7^-1 * 11^-1 * 13^-1,
 2^-1 * 3^-7 * 5^-4 * 7^-2 * 11^-1 * 13^-1 * 17^-1 * 3617,
 3^-9 * 5^-3 * 7^-3 * 11^-1 * 13^-1 * 17^-1 * 19^-1 * 43867]