OFFSET
1,7
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..10000
PROG
(Python)
from math import prod
from sympy import factorint, primepi
def A364558(n): return (1<<primepi(max(f:=factorint(n)))+sum(e-1 for e in f.values())-1)-prod(p**(e-1)*(p-1) for p, e in f.items()) if n>1 else 0 # Chai Wah Wu, Jul 29 2023
CROSSREFS
KEYWORD
sign
AUTHOR
Antti Karttunen, Jul 28 2023
STATUS
approved