OFFSET
1,2
COMMENTS
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 1..12916 (a(n) < 10^7)
Vaclav Kotesovec, Plot of a(n)/n^(3/2) for n = 1..12916
EXAMPLE
lcm(tau(6), sigma(6), pod(6)) = lcm(4, 12, 36) = 36 = pod(6).
MATHEMATICA
Select[Range[3000], LCM @@ {(d = DivisorSigma[0, #]), DivisorSigma[1, #], (pod = #^(d/2))} == pod &] (* Amiram Eldar, Oct 24 2020 *)
PROG
(Magma) [m: m in [1..10^5] | LCM([#Divisors(m), &+Divisors(m), &*Divisors(m)]) eq &*Divisors(m)]
(PARI) isok(m) = my(d=divisors(m), prd=vecprod(d)); lcm([#d, vecsum(d), prd]) == prd; \\ Michel Marcus, Oct 24 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Jaroslav Krizek, Oct 23 2020
STATUS
approved