ClearAll[SmarandachePrimeQ] SmarandachePrimeQ[n_Integer] := MatchQ[IntegerDigits[n], {(2 | 3 | 5 | 7) ..}] \[And] PrimeQ[n] s = Select[Range[10^5], SmarandachePrimeQ]; Take[s, UpTo[25]] s[[100]]