OFFSET
0,3
COMMENTS
A number k has "equal prime multiplicities" (or is "uniform") iff its prime signature is constant, meaning that k is a power of a squarefree number.
LINKS
David A. Corneth, Table of n, a(n) for n = 0..9999
Jon Maiga, Computer-generated formulas for A336415, Sequence Machine.
FORMULA
a(n) = A327527(n!).
EXAMPLE
The a(n) uniform divisors of n for n = 1, 2, 6, 8, 30, 36 are the columns:
1 2 6 8 30 36
1 3 6 15 30
2 4 10 16
1 3 8 15
2 6 10
1 5 9
4 8
3 6
2 5
1 4
3
2
1
In 20!, the multiplicity of the third prime (5) is 4 but the multiplicity of the fourth prime (7) is 2. Hence there are 2^3 - 1 = 3 divisors with all exponents 3 (we subtract |{1}| = 1 from that count as 1 has no exponent 3). - David A. Corneth, Jul 27 2020
MATHEMATICA
Table[Length[Select[Divisors[n!], SameQ@@Last/@FactorInteger[#]&]], {n, 0, 15}]
PROG
(PARI) a(n) = sumdiv(n!, d, my(ex=factor(d)[, 2]); (#ex==0) || (vecmin(ex) == vecmax(ex))); \\ Michel Marcus, Jul 24 2020
(PARI) a(n) = {if(n<2, return(1)); my(f = primes(primepi(n)), res = 1, t = #f); f = vector(#f, i, val(n, f[i])); for(i = 1, f[1], while(f[t] < i, t--; ); res+=(1<<t - 1) ); res }
val(n, p) = my(r=0); while(n, r+=n\=p); r \\ David A. Corneth, Jul 27 2020
CROSSREFS
The version for distinct prime multiplicities is A336414.
The version for nonprime perfect powers is A336416.
Uniform partitions are counted by A047966.
Numbers with distinct prime multiplicities are A130091.
Divisors with distinct prime multiplicities are counted by A181796.
Maximum divisor with distinct prime multiplicities is A327498.
Uniform divisors are counted by A327527.
Maximum uniform divisor is A336618.
1st differences are given by A048675.
KEYWORD
nonn
AUTHOR
Gus Wiseman, Jul 22 2020
EXTENSIONS
Terms a(31) and onwards from David A. Corneth, Jul 27 2020
STATUS
approved