RosettaCodeData/Task/Smith-numbers/Mathematica/smith-numbers.math

4 lines
174 B
Plaintext

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