%I #4 Jul 14 2012 13:52:01
%S 4621,34361,2881861,63882851
%N 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.
%e 63882851=prime((6*3*8*8*2*8*5*1)*(6+3+8+8+2+8+5+1)).
%t 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_ *)
%Y Cf. A097223.
%K base,nonn
%O 1,1
%A _Farideh Firoozbakht_, Jun 09 2008