OFFSET
1,1
LINKS
PROG
(PARI)
A003961(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
A351555(n) = { my(s=sigma(n), f=factor(s), u=A003961(n)); sum(k=1, #f~, if((f[k, 1]%2) && 0!=(u%f[k, 1]), (valuation(n, f[k, 1])!=f[k, 2]), 0)); };
isA387160(n) = (isprime(core(n)) && (0==A351555(n)));
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Aug 24 2025
EXTENSIONS
Terms a(29)-a(34) from Giovanni Resta, Aug 25 2025
STATUS
approved
