login
A336921
a(n) = A331410(n) - A087436(n).
6
0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 1, 0, 1, 0, 2, 0, 2, 1, 0, 1, 1, 0, 2, 1, 0, 0, 3, 1, 0, 0, 1, 2, 1, 0, 3, 2, 1, 1, 2, 0, 2, 1, 1, 1, 1, 0, 0, 2, 2, 1, 3, 0, 2, 0, 2, 3, 3, 1, 1, 0, 0, 0, 2, 1, 3, 2, 1, 1, 2, 0, 4, 3, 2, 2, 1, 1, 2, 1, 0, 2, 2, 0, 3, 2, 3, 1, 4, 1, 1, 1, 0, 1, 3, 0, 2, 0, 1, 2, 4, 2, 2, 1, 1
OFFSET
1,17
COMMENTS
Totally additive because both A087436 and A331410 are.
LINKS
FORMULA
a(n) = A331410(A336467(n)).
a(n) = A331410(n) - A087436(n).
a(n) = A336922(n) - A046660(A000265(n)).
For all n >= 1, a(n) <= A336118(n).
PROG
(PARI)
A000265(n) = (n>>valuation(n, 2));
A331410(n) = if(!bitand(n, n-1), 0, 1+A331410(n+(n/vecmax(factor(n)[, 1]))));
A336467(n) = { my(f=factor(n)); prod(k=1, #f~, if(2==f[k, 1], 1, (A000265(f[k, 1]+1))^f[k, 2])); };
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Aug 09 2020
STATUS
approved