|
%I #19 Jun 20 2024 22:50:25
%S 1,1,1,1,1,1,1,1,1,1,1,6,1,1,1,1,1,6,1,10,1,1,1,36,1,1,1,14,1,1,1,1,1,
%T 1,1,1,1,1,1,100,1,1,1,22,15,1,1,216,1,10,1,26,1,36,1,196,1,1,1,30,1,
%U 1,21,1,1,1,1,34,1,1,1,6,1,1,15,38,1,1,1,1000
%N a(n) = A062760(n)*A066636(n).
%C If n has unequal prime exponents (a term in A059404), then a(n) > 1; otherwise a(n) = 1.
%H Michael De Vlieger, <a href="/A373703/b373703.txt">Table of n, a(n) for n = 1..10000</a>
%F For n > 1, a(n) = A007947(n)^k where k is the difference between the greatest and least exponents in the prime power factorization of n.
%e a(12) = A062760(12)*A066636(12) = 2*3 = 6.
%e a(45) = A062760(45)*A066636(45) = 3*5 = 15.
%t Table[Function[{r, s}, r^(Max[s] - Min[s])] @@ {Times @@ #[[All, 1]], #[[All, -1]]} &@ FactorInteger[n], {n, 120}] (* _Michael De Vlieger_, Jun 13 2024 *)
%o (PARI) a(n) = if (n==1, 1, my(f=factor(n)); (factorback(f[, 1]))^(vecmax(f[,2])-vecmin(f[,2]))); \\ _Michel Marcus_, Jun 14 2024
%Y Cf. A007947, A059404, A062760, A066636.
%K nonn
%O 1,12
%A _David James Sycamore_, Jun 13 2024
|
