A078547 a(n) = lcm(n, A052429(n)) - n.

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, 28, 493, 0, 62, 64, 0, 170, 70, 0, 740, 418, 78, 0, 123, 42, 473, 0, 135, 230, 1269, 0, 1715, 0, 204, 208, 742, 486, 0, 784, 1938, 1102, 2596, 0, 305, 124, 63, 128, 325, 0

Offset: 1

Author

Labos Elemer, Dec 05 2002

Keywords

• nonn
• base

Comments

a(n)=0 if n is divisible by each nonzero digit, i.e., if n in A002796.

Crossrefs

Cf. A052429, A078546, A078548, A002796.

Programs

Mathematica

lc[x_] := Apply[LCM, DeleteCases[IntegerDigits[x], 0]] Table[LCM[lc[w], w]-w, {w, 1, 128}]

PARI

lcnzd(n) = lcm(select(x->(x!=0), digits(n)));
a(n) = lcm(n, lcnzd(n)) - n; \\ Michel Marcus, Mar 18 2018