[go: up one dir, main page]

login
A333793
a(n) = A333794(n) - A073934(n).
3
0, 0, 0, 0, 0, 1, 1, 0, 1, 3, 3, 3, 3, 8, 5, 0, 0, 6, 6, 9, 12, 18, 18, 7, 9, 18, 6, 22, 22, 18, 18, 0, 30, 15, 24, 16, 16, 33, 28, 21, 21, 37, 37, 48, 24, 69, 69, 15, 37, 36, 29, 48, 48, 25, 54, 50, 49, 77, 77, 44, 44, 73, 49, 0, 56, 83, 83, 45, 113, 75, 75, 36, 36, 71, 54, 87, 87, 81, 81, 45, 25, 84, 84, 87, 57, 128, 119, 108, 108, 71
OFFSET
1,10
LINKS
FORMULA
a(n) = A333794(n) - A073934(n).
a(p) = a(p-1), for all primes p.
a(A000079(n)) = a(A019434(n)) = 0, for all applicable n.
PROG
(PARI)
A073934(n) = if(1==n, n, n + A073934(n-(n/vecmin(factor(n)[, 1]))));
A333794(n) = if(1==n, n, n + A333794(n-(n/vecmax(factor(n)[, 1]))));
A333793(n) = (A333794(n)-A073934(n));
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Apr 05 2020
STATUS
approved