%I #13 Dec 30 2017 09:27:11
%S 0,2,3,4,5,30,7,24,18,70,11,48,13,126,120,32,17,126,19,120,210,286,23,
%T 264,50,390,27,224,29,930,31,160,462,646,420,180,37,798,624,680,41,
%U 1722,43,528,585,1150,47,336,98,450,1020,728,53,162,880,1288,1254,1798,59
%N a(n) = lcm(n, A003415(n)).
%H Michael De Vlieger, <a href="/A086130/b086130.txt">Table of n, a(n) for n = 1..10000</a>
%e If n=prime, lcm(prime,1)=prime.
%t 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 *)
%o (PARI) a(n) = lcm(n, sum(i=1, #f=factor(n)~, n/f[1, i]*f[2, i]) ); \\ _Michel Marcus_, Dec 11 2013
%Y Cf. A003415.
%K nonn
%O 1,2
%A _Labos Elemer_, Jul 23 2003
|