\\ positive integer k only: this finite product is not valid for half-integers
default(realprecision, 140)
k = 3
prodeulerrat((1 - 1/p)^((k - 1)^2) *
sum(j = 0, k - 1, binomial(k - 1, j)^2/p^j), 1, 2)The generator computes the small-prime local factors directly with arb's hypergeometric function, computes the middle prime range from the exact power-series coefficients of the logarithm of the local factor in (2), and computes the remaining prime tail from prime-zeta values.