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

Author

Antti Karttunen, Dec 11 2019

Keywords

• nonn

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

Crossrefs

Cf. A002322, A051903, A327295.

Programs

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