OFFSET
1,1
EXAMPLE
Sequence includes all primes with sum of digits > 2, i.e., all primes not of the form 10^d + 1 for nonnegative integers d (the only known primes of this form are 2, 11, and 101; see A185121).
Sequence also includes all squared primes with sum of digits > 3.
MATHEMATICA
sud[x_] := Apply[Plus, IntegerDigits[x]] Do[s=sud[n]-DivisorSigma[0, n]; If[s>0, Print[n]], {n, 1, 256}]
Select[Range[100], Total[IntegerDigits[#]]>DivisorSigma[0, #]&] (* Harvey P. Dale, Mar 07 2020 *)
CROSSREFS
KEYWORD
base,easy,nonn
AUTHOR
Labos Elemer, Sep 26 2002
EXTENSIONS
Edited by Jon E. Schoenfield, Sep 23 2018
STATUS
approved