OFFSET
0,3
COMMENTS
A number is a perfect power iff it is 1 or its prime exponents (signature) are not relatively prime.
The n-th superprimorial number is A006939(n) = Product_{i = 1..n} prime(i)^(n - i + 1).
LINKS
Andrew Howroyd, Table of n, a(n) for n = 0..200
FORMULA
a(n) = 1 + Sum_{k=2..n} mu(k)*(1 - Product_{i=1..n} 1 + floor(i/k)). - Andrew Howroyd, Aug 30 2020
EXAMPLE
The a(0) = 1 through a(4) = 15 divisors:
1 2 12 360 75600
-------------------------
1 1 1 1 1
4 4 4
8 8
9 9
36 16
25
27
36
100
144
216
225
400
900
3600
MATHEMATICA
chern[n_]:=Product[Prime[i]^(n-i+1), {i, n}];
perpouQ[n_]:=Or[n==1, GCD@@FactorInteger[n][[All, 2]]>1];
Table[Length[Select[Divisors[chern[n]], perpouQ]], {n, 0, 5}]
PROG
(PARI) a(n) = {1 + sum(k=2, n, moebius(k)*(1 - prod(i=1, n, 1 + i\k)))} \\ Andrew Howroyd, Aug 30 2020
CROSSREFS
A000325 is the uniform version.
A336416 gives the same for factorials instead of superprimorials.
A000217 counts prime power divisors of superprimorials.
A000961 gives prime powers.
A006939 gives superprimorials or Chernoff numbers.
A022915 counts permutations of prime indices of superprimorials.
A091050 counts perfect power divisors.
A181818 gives products of superprimorials.
A294068 counts factorizations using perfect powers.
A317829 counts factorizations of superprimorials.
KEYWORD
nonn
AUTHOR
Gus Wiseman, Jul 24 2020
EXTENSIONS
Terms a(10) and beyond from Andrew Howroyd, Aug 30 2020
STATUS
approved