STATUS
reviewed
approved
reviewed
approved
proposed
reviewed
editing
proposed
(Python 3.8+)
proposed
editing
editing
proposed
return sum(s*e//p for p, e in sum(((lambda x: Counter({k:x[k]*m for k in x}))(factorint(m)) for m in range(2, n+1)), start=Counter({2:0})).items()) if n > 1 else 0 # Chai Wah Wu, Jun 12 2022
(Python)
from math import prod
from collections import Counter
from sympy import factorint
def A260619(n):
s = prod(i**i for i in range(2, n+1))
return sum(s*e//p for p, e in sum(((lambda x: Counter({k:x[k]*m for k in x}))(factorint(m)) for m in range(2, n+1)), start=Counter({2:0})).items()) if n > 1 else 0 # Chai Wah Wu, Jun 12 2022
approved
editing
reviewed
approved
proposed
reviewed
editing
proposed