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,5

Links

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

Formula

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

Prog

(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)));

Crossrefs

Cf. A163511, A331410, A365385.

Cf. also A334873.

Sequence in context: A287647 A331749 A274438 * A087565 A079163 A108536

Adjacent sequences: A365384 A365385 A365386 * A365388 A365389 A365390

Keyword

nonn

Author

Antti Karttunen, Sep 03 2023

Status

approved