login
a(n) = lcm(n, A052429(n)) - n.
2

%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