OFFSET
1,2
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
EXAMPLE
If n=prime, lcm(prime,1)=prime.
MATHEMATICA
Table[LCM[n, If[Abs@ n < 2, 0, n Total[#2/#1 & @@@ FactorInteger[Abs@ n]]]], {n, 59}] (* Michael De Vlieger, Feb 05 2017, after Michael Somos at A003415 *)
PROG
(PARI) a(n) = lcm(n, sum(i=1, #f=factor(n)~, n/f[1, i]*f[2, i]) ); \\ Michel Marcus, Dec 11 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, Jul 23 2003
STATUS
approved