OFFSET
1,2
COMMENTS
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
Michael De Vlieger, Log log scatterplot of a(n), n = 1..10^5.
FORMULA
EXAMPLE
a(1) = 1 since s(1) = 1 = 1*1. (We count divisor 1 only once, since it is a square root.)
a(2) = 2 since s(2) = 4 = 2*2, and omega(2) = omega(4). Alternatively, rad(2) = rad(4).
a(3) = 6 since s(3) = 8 = 2*4, omega(2) = omega(4) = omega(8), i.e., rad(2) = rad(4) = rad(8), and 2+4 = 6.
a(4) = 3 since s(4) = 9 = 3*3, and omega(3) = omega(9); rad(3) = rad(9).
a(5) = 14 since s(5) = 16 = 2*8 = 4*4; 2+4+8 = 14.
a(6) = 5 since s(6) = 25 = 5*5.
a(7) = 12 since s(7) = 27 = 3*9; 3+9 = 12.
a(8) = 30 since s(8) = 32 = 2*16 = 4*8; 2+4+8+16 = 30.
a(9) = 6 since s(9) = 36 = 6*6, etc.
MATHEMATICA
nn = 1500;
s = Union@ Flatten@ Table[a^2*b^3, {b, Surd[nn, 3]}, {a, Sqrt[nn/b^3]}];
Map[Function[{n, r},
DivisorSum[n, # &, PrimeNu[#] == PrimeNu[n/#] == r &]] @@
{#, PrimeNu[#]} &, s]
CROSSREFS
KEYWORD
nonn
AUTHOR
Michael De Vlieger, Jan 15 2025
STATUS
approved