%I #10 Jul 16 2021 11:12:43
%S 3,115,6576,80855
%N Numbers k such that k = (d_1! + d_2! + ... + d_k!) - (d_1 + d_2 + ... + d_k) where d_1 d_2 ... d_k is the decimal expansion of k.
%e 80855 is in the sequence because 80855 = (8! + 0! + 8! + 5! + 5!) - (8 + 0 + 8 + 5 + 5).
%t okQ[n_]:=Module[{idn=IntegerDigits[n]},n==Total[idn!]-Total[idn]]; Select[Range[81000],okQ] (* _Harvey P. Dale_, Aug 26 2011 *)
%Y Cf. A097643.
%K base,fini,full,nice,nonn
%O 1,1
%A _Farideh Firoozbakht_, Aug 22 2004
|