A329885 a(n) = A051903(n) mod A002322(n), where A051903 gives the maximal prime exponent of n, and A002322 is Carmichael's lambda (also known as psi).

0, 0, 1, 0, 1, 1, 1, 1, 2, 1, 1, 0, 1, 1, 1, 0, 1, 2, 1, 2, 1, 1, 1, 1, 2, 1, 3, 2, 1, 1, 1, 5, 1, 1, 1, 2, 1, 1, 1, 3, 1, 1, 1, 2, 2, 1, 1, 0, 2, 2, 1, 2, 1, 3, 1, 3, 1, 1, 1, 2, 1, 1, 2, 6, 1, 1, 1, 2, 1, 1, 1, 3, 1, 1, 2, 2, 1, 1, 1, 0, 4, 1, 1, 2, 1, 1, 1, 3, 1, 2, 1, 2, 1, 1, 1, 5, 1, 2, 2, 2, 1, 1, 1, 3, 1

Offset

1,9

Comments

This differs from A051903 at n = 2, 4, 8, 12, 16, 24, 48, 80, 240. Are there any other such n? (None other found <= 201326592.)

Links

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

Prog

(PARI)
A002322(n) = lcm(znstar(n)[2]); \\ From A002322
A051903(n) = if((1==n), 0, vecmax(factor(n)[, 2]));
A329885(n) = (A051903(n)%A002322(n));

Crossrefs

Cf. A002322, A051903, A327295.

Sequence in context: A286133 A328248 A360158 * A085737 A364249 A191904

Adjacent sequences: A329882 A329883 A329884 * A329886 A329887 A329888

Keyword

nonn

Author

Antti Karttunen, Dec 11 2019

Status

approved