A140175 Numbers n such that n=prime((d_1*d_2*...*d_k)*(d_1+d_2+...+d_k)) where d_1d_2...d_k is the decimal expansion of n.
4621, 34361, 2881861, 63882851
Offset: 1
Examples
63882851=prime((6*3*8*8*2*8*5*1)*(6+3+8+8+2+8+5+1)).
Crossrefs
Cf. A097223.
Programs
-
Mathematica
fQ[n_] := Block[{id = IntegerDigits@ n}, n == Prime[Plus @@ id Times @@ id]]; k = 1; lst = {}; While[ k < 10^8, If[ fQ@n, AppendTo[lst, k]; Print@ n]; k++ ]; lst (* Robert G. Wilson v *)