A336922 a(n) = A331410(n) - A005087(n).

0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 1, 0, 1, 0, 2, 1, 2, 1, 0, 1, 1, 0, 3, 1, 2, 0, 3, 1, 0, 0, 1, 2, 1, 1, 3, 2, 1, 1, 2, 0, 2, 1, 2, 1, 1, 0, 1, 3, 2, 1, 3, 2, 2, 0, 2, 3, 3, 1, 1, 0, 1, 0, 2, 1, 3, 2, 1, 1, 2, 1, 4, 3, 3, 2, 1, 1, 2, 1, 3, 2, 2, 0, 3, 2, 3, 1, 4, 2, 1, 1, 0, 1, 3, 0, 2, 1, 2, 3, 4, 2, 2, 1, 1

Offset: 1

Author

Antti Karttunen, Aug 09 2020

Keywords

• nonn

Links

• Antti Karttunen, Table of n, a(n) for n = 1..65537

Crossrefs

Cf. A000265, A005087, A054784 (positions of zeros), A331410, A336467, A336921, A336923.

Cf. also A336469.

Programs

PARI

A005087(n) = omega(n>>valuation(n,2));
A331410(n) = if(!bitand(n,n-1),0,1+A331410(n+(n/vecmax(factor(n)[, 1]))));
A336922(n) = (A331410(n) - A005087(n));

Formula

a(n) = A331410(n) - A005087(n).

a(n) = A336921(n) + A046660(A000265(n)).