A365387 a(n) = A331410(n) * A331410(A163511(n)).

0, 0, 1, 0, 4, 1, 2, 0, 6, 4, 8, 1, 6, 2, 3, 0, 12, 6, 18, 4, 10, 8, 4, 1, 16, 6, 9, 2, 8, 3, 2, 0, 15, 12, 24, 6, 28, 18, 9, 4, 18, 10, 12, 8, 12, 4, 8, 1, 10, 16, 20, 6, 16, 9, 12, 2, 12, 8, 16, 3, 6, 2, 6, 0, 24, 15, 40, 12, 27, 24, 12, 6, 40, 28, 25, 18, 12, 9, 18, 4, 28, 18, 18, 10, 25, 12, 25, 8, 20, 12, 18, 4

Offset: 1

Author

Antti Karttunen, Sep 03 2023

Keywords

• nonn

Links

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

Crossrefs

Cf. A163511, A331410, A365385.

Cf. also A334873.

Programs

PARI

A163511(n) = if(!n, 1, my(p=2, t=1); while(n>1, if(!(n%2), (t*=p), p=nextprime(1+p)); n >>= 1); (t*p));
A331410(n) = if(!bitand(n,n-1),0,1+A331410(n+(n/vecmax(factor(n)[, 1]))));
A365387(n) = (A331410(n)*A331410(A163511(n)));

Formula

a(n) = A331410(n) * A365385(n).