|
|
A072696
|
|
a(n) = lcm(d(n^3), d(n)), where d(n) = A000005, the number of divisors of n.
|
|
2
|
|
|
1, 4, 4, 21, 4, 16, 4, 20, 21, 16, 4, 84, 4, 16, 16, 65, 4, 84, 4, 84, 16, 16, 4, 40, 21, 16, 20, 84, 4, 64, 4, 48, 16, 16, 16, 441, 4, 16, 16, 40, 4, 64, 4, 84, 84, 16, 4, 260, 21, 84, 16, 84, 4, 40, 16, 40, 16, 16, 4, 336, 4, 16, 84, 133, 16, 64, 4, 84, 16, 64, 4, 420, 4, 16
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
If n is a product of k distinct primes, then a(n) = 4^k.
|
|
LINKS
|
|
|
MATHEMATICA
|
Table[LCM[DivisorSigma[0, n^3], DivisorSigma[0, n]], {n, 80}] (* Wesley Ivan Hurt, Nov 25 2017 *)
|
|
PROG
|
(PARI) a(n) = {my(e = factor(n)[, 2]); lcm(vecprod(apply(x -> 3*x+1, e)), vecprod(apply(x -> x+1, e))); } \\ Amiram Eldar, Dec 02 2023
|
|
CROSSREFS
|
|
|
KEYWORD
|
easy,nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|