|
%I #18 Mar 26 2020 06:46:14
%S 0,0,0,0,0,0,0,0,0,0,0,0,26,14,0,32,102,54,152,0,21,0,115,0,25,52,351,
%T 28,493,0,62,64,0,170,70,0,740,418,78,0,123,42,473,0,135,230,1269,0,
%U 1715,0,204,208,742,486,0,784,1938,1102,2596,0,305,124,63,128,325,0
%N a(n) = lcm(n, A052429(n)) - n.
%C a(n)=0 if n is divisible by each nonzero digit, i.e., if n in A002796.
%t lc[x_] := Apply[LCM, DeleteCases[IntegerDigits[x], 0]] Table[LCM[lc[w], w]-w, {w, 1, 128}]
%o (PARI) lcnzd(n) = lcm(select(x->(x!=0), digits(n)));
%o a(n) = lcm(n, lcnzd(n)) - n; \\ _Michel Marcus_, Mar 18 2018
%Y Cf. A052429, A078546, A078548, A002796.
%K nonn,base
%O 1,13
%A _Labos Elemer_, Dec 05 2002
|
