STATUS
reviewed
approved
reviewed
approved
proposed
reviewed
editing
proposed
(Python)
from sympy import factorint
def a(n): return (n*(n+1)//2)%sum(p*e for p, e in factorint(n).items())
print([a(n) for n in range(2, 98)]) # Michael S. Branicky, Apr 24 2022
approved
editing
proposed
approved
editing
proposed
Robert Israel, <a href="/A352996/b352996.txt">Table of n, a(n) for n = 2..10000</a>
approved
editing
proposed
approved
editing
proposed
a(n) = n*(n+1)/2 mod the (sum (with multiplicity) of prime factors of n).