login
A366902
The number of exponentially evil divisors of n.
9
1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1
OFFSET
1,8
COMMENTS
First differs from A050361 at n = 128.
The number of divisors of n that are exponentially evil numbers (A262675), i.e., numbers having only evil (A001969) exponents in their canonical prime factorization.
The sum of these divisors is A366904(n) and the largest of them is A366906(n).
LINKS
FORMULA
Multiplicative with a(p^e) = A159481(e).
a(n) >= 1, with equality if and only if n is a cubefree number (A004709).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Product_{p prime} Sum_{k>=1} 1/p^A262675(k) = 1.241359937856... .
MATHEMATICA
f[p_, e_] := Floor[e/2] + If[OddQ[e] || OddQ[DigitCount[e + 1, 2, 1]], 1, 0]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]
PROG
(PARI) s(n) = n\2 + (n%2 || hammingweight(n+1)%2); \\ after Charles R Greathouse IV at A159481
a(n) = vecprod(apply(x -> s(x), factor(n)[, 2]));
CROSSREFS
KEYWORD
nonn,easy,mult
AUTHOR
Amiram Eldar, Oct 27 2023
STATUS
approved