RosettaCodeData/Task/Smith-numbers/Mathematica/smith-numbers.math
2023-07-01 13:44:08 -04:00

3 lines
174 B
Text

smithQ[n_] := Not[PrimeQ[n]] &&
Total[IntegerDigits[n]] == Total[IntegerDigits /@ Flatten[ConstantArray @@@ FactorInteger[n]],2];
Select[Range[2, 10000], smithQ]