OFFSET
1,2
COMMENTS
a(2) = 2 is the only prime term possible, since the product of tau(p^i) is always even, and 2 is the only even prime. - Michael De Vlieger, Mar 27 2016
The corresponding k are: 1, 1, 2, 2, 2, 2, 2, 3, 4, 4, 4, 5, 5, 5. - Michel Marcus, Apr 08 2016; updated by Max Alekseyev, Jun 11 2025
LINKS
Max Alekseyev, Table of n, a(n) for n = 1..14
EXAMPLE
d(4050) * d(4050^2) = 30 * 135 = 4050;
d(66528) * d(66528^2) = 96 * 693 = 66528.
MAPLE
with(numtheory): P:=proc(q) local a, k, n;
for n from 1 to q do a:=tau(n); k:=1;
while a<n do k:=k+1; a:=a*tau(n^k); od;
if n=a then print(n); fi; od; end: P(10^6);
MATHEMATICA
Select[Insert[Complement[Range@ #, Prime@ Range@ PrimePi@ #] &[2 10^5], 2, 2], Function[k, AnyTrue[Range@ 3, Product[DivisorSigma[0, k^i], {i, #}] == k &]]] (* Michael De Vlieger, Mar 25 2016 *)
PROG
(PARI) isok(m) = my(k = 1, prd = 1); while (prd < m, prd *= numdiv(m^k); k++); prd == m; \\ Michel Marcus, Apr 08 2016, Jun 12 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Paolo P. Lava, Mar 22 2016
EXTENSIONS
a(8)-a(10) from Hiroaki Yamanouchi, Apr 07 2016
a(11)-a(14) from Max Alekseyev, Jun 10 2025
STATUS
approved