OFFSET
1,1
LINKS
Giovanni Resta, Table of n, a(n) for n = 1..95 (terms < 3.5*10^14)
EXAMPLE
8 is a term as p = 8 and 8*8 = 64 is divisible by 8+8 = 16.
3276 is a term as p = 3*2*7*6 = 252 and 3276*252 = 825552 is divisible by 3276+252 = 3528.
3787749 is a term as p = 3*7*8*7*7*4*9 = 296352 and 3787749*296352 = 1122506991648 is divisible by 3787749+296352 = 4084101.
MATHEMATICA
Select[Range[10^6], (p = Times @@ IntegerDigits@ #; p > 0 && Mod[# p, # + p] == 0) &] (* Giovanni Resta, May 08 2020 *)
PROG
(PARI) isok(m) = my(p=vecprod(digits(m))); p && !((m*p) % (m+p)); \\ Michel Marcus, May 08 2020
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Scott R. Shannon, May 08 2020
STATUS
approved