OFFSET
1,4
COMMENTS
Completely determined by the exponents >=2 in the prime factorization of n (cf. A212172).
Not the same as the number of odd divisors of n (A001227(n)); see example.
Multiplicative because A000005 is multiplicative and A000265 is completely multiplicative. - Andrew Howroyd, Aug 01 2018
a(n) = 1 iff the number of divisors of n is a power of 2 (A036537). - Bernard Schott, Nov 04 2022
LINKS
FORMULA
From Antti Karttunen, Jan 14 2020: (Start)
(End)
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Product_{p odd prime} ((1 - 1/p)*(1 + Sum_{k>=1} a(k+1)/p^k)) = 2.076325817863586... . - Amiram Eldar, Oct 15 2022
EXAMPLE
48 has a total of 10 divisors (1, 2, 3, 4, 6, 8, 12, 16, 24 and 48). Since the largest odd divisor of 10 is 5, a(48) = 5.
MATHEMATICA
Table[Block[{nd=DivisorSigma[0, n]}, nd/2^IntegerExponent[nd, 2]], {n, 100}] (* Indranil Ghosh, Jul 19 2017, after PARI code *)
PROG
(PARI) a(n) = my(nd = numdiv(n)); nd/2^valuation(nd, 2); \\ Michel Marcus, Jul 19 2017
(Python)
from sympy import divisor_count, divisors
def a(n): return [i for i in divisors(divisor_count(n)) if i%2][-1]
print([a(n) for n in range(1, 101)]) # Indranil Ghosh, Jul 19 2017
CROSSREFS
KEYWORD
nonn,mult
AUTHOR
Matthew Vandermast, Jun 04 2012
STATUS
approved