a30 := primorial(42681)*primorial(151)*primorial(29)*primorial(11)
*primorial(7)*primorial(5)*primorial(4)*(primorial(3)^2)
*(primorial(2)^4)*(primorial(1)^9);
? 10000*Qr(a30)
%370 = 9999.10909330009471428528222805072598431749321026349669841
Qr(W) = sigma(W)/(harmonic(W)+log(harmonic(W))*exp(harmonic(W)))
harmonic(Z)=Euler+psi(Z+1)
primorial(W)=prod(Z=1,W,prime(Z))