login
A391074
a(n) is the number of nonuniform divisors of n.
1
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 2, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 2, 0, 0, 0, 1, 1, 0, 0, 3, 0, 1, 0, 1, 0, 2, 0, 2, 0, 0, 0, 3, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 4, 0, 0, 1, 1, 0, 0, 0, 3, 0, 0, 0, 3, 0, 0, 0, 2
OFFSET
1,24
COMMENTS
This sequence first differs from A062977, A335460 and A376514 at n = 36: A062977(36) = 0, A335460(36) = 4, A376514(36) = 3 while a(36) = 2 = card({12, 18}).
FORMULA
a(n) = A000005(n) - A327527(n).
a(n) = A000005(n) - A034444(n) - A368251(n).
a(n) <= A376514(n).
MAPLE
A391074 := proc(n) local f, e, k, t, m, M; f := ifactors(n)[2]; e := map(l -> l[2], f); k := nops(f); M := max(e); t := 0; for m from 2 to M do t := t + 2^nops(select(x -> m <= x, e)) - 1; end do; return NumberTheory:-tau(n) - 2^k - t; end proc: seq(A391074(n), n = 1 .. 88);
MATHEMATICA
a[n_] := Module[{e = FactorInteger[n][[;; , 2]]}, Times @@ (e+1) - 1 - Total[2^Accumulate[Count[e, #] & /@ Range[Max[e], 1, -1]] - 1]]; a[1] = 0; Array[a, 100] (* Amiram Eldar, Dec 03 2025 *)
KEYWORD
nonn,easy
AUTHOR
Felix Huber, Dec 03 2025
EXTENSIONS
Incorrect formula removed by Falk Hüffner, May 19 2026
STATUS
approved