OFFSET
1,3
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..65537
MATHEMATICA
f[p_, e_] := If[p == 2, 1, (p + 1)*p^(e - 1)]; psiOdd[1] = 1; psiOdd[n_] := Times @@ f @@@ FactorInteger[n]; a[n_] := IntegerExponent[NestWhile[psiOdd, n, # != 2^IntegerExponent[#, 2] &], 2]; Array[a, 100] (* Amiram Eldar, Aug 31 2021 *)
PROG
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Aug 31 2021
STATUS
approved