A101699 Numbers m such that phi(m) = d_1*d_1!+d_2*d_2!+...+d_k*d_k! where d_1 d_2 ... d_k is the decimal expansion of m.
1, 23, 66666, 13179592
Offset: 1
Examples
13179592 is in the sequence because phi(13179592) = 1*1!+3*3!+1*1!+7*7!+9*9!+5*5!+9*9!+2*2!.
Programs
-
Mathematica
Do[If[EulerPhi[n] == Apply[Plus, IntegerDigits[n]*IntegerDigits[n]! ], Print[n]], {n, 30000000}] Select[Range[132*10^5],EulerPhi[#]==Total[# #!&/@IntegerDigits[#]]&] (* Harvey P. Dale, Mar 06 2023 *)