|
%I #12 Mar 18 2018 04:04:15
%S 0,0,0,0,0,0,0,0,0,9,10,10,36,24,10,42,112,64,162,18,40,20,132,20,40,
%T 72,364,48,504,27,90,90,30,192,90,30,756,432,108,36,160,80,504,40,160,
%U 264,1288,40,1728,45,250,250,780,520,50,810,1960,1120,2610,54,360,180
%N a(n) = lcm(n, A052429(n)) - A052429(n).
%t lc[x_] := Apply[LCM, DeleteCases[IntegerDigits[x], 0]] Table[LCM[lc[w], w]-lc[w], {w, 1, 128}]
%o (PARI) lcnzd(n) = lcm(select(x->(x!=0), digits(n)));
%o a(n) = my(lc=lcnzd(n)); lcm(n, lc) - lc; \\ _Michel Marcus_, Mar 18 2018
%Y Cf. A052429, A078546, A078547, A002796.
%K nonn,base
%O 1,10
%A _Labos Elemer_, Dec 05 2002
|
