(PARI) isok(n) = if(n, factorback(digits(n)), 0) * sumdigits(n) == 2*n \\
Mohammed Yaseen, Jul 22 2022
(Python)
from math import prod
def s(n): return sum(map(int, str(n)))
def p(n): return prod(map(int, str(n)))
for n in range(0, 10**6):
if p(n)*s(n)==2*n: