reviewed
approved
reviewed
approved
proposed
reviewed
editing
proposed
return (31*prod((p**(5*(e+1))-1)//(p**5-1) for p, e in f)-70*(n+1)*prod((p**(3*(e+1))-1)//(p**3-1) for p, e in f) + (20*n*((n<<1)+3)+9)*prod((p**(e+1)-1)//(p-1) for p, e in f))//1920 # Chai Wah Wu, Jul 24 2024
(Python)
from math import prod
from sympy import factorint
def A374930(n):
f = factorint(n).items()
return (31*prod((p**(5*(e+1))-1)//(p**5-1) for p, e in f)-70*(n+1)*prod((p**(3*(e+1))-1)//(p**3-1) for p, e in f) + (20*n*((n<<1)+3)+9)*prod((p**(e+1)-1)//(p-1) for p, e in f))//1920 # Chai Wah Wu, Jul 24 2024
approved
editing
reviewed
approved
proposed
reviewed
editing
proposed
A374930[n_] := (31*DivisorSigma[5, n] - 70*(n + 1)*DivisorSigma[3, n] + (40*n^2 + 60*n + 9)*DivisorSigma[1, n])/1920; Array[A374930, 50, 3] (* _Paolo Xausa_, Jul 24 2024 *)
Array[A374930, 50, 3] (* Paolo Xausa, Jul 24 2024 *)
A374930[n_] := (31*DivisorSigma[5, n] - 70*(n + 1)*DivisorSigma[3, n] + (40*n^2 + 60*n + 9)*DivisorSigma[1, n])/1920; Array[A374930, 50, 3] (* Paolo Xausa, Jul 24 2024 *)
approved
editing