login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A095403 Sum of digits of n minus the sum of digits of all distinct prime factors of n. 3
1, 0, 0, 2, 0, 1, 0, 6, 6, -6, 0, -2, 0, -4, -2, 5, 0, 4, 0, -5, -7, 0, 0, 1, 2, 2, 6, 1, 0, -7, 0, 3, 1, -3, -4, 4, 0, -1, 5, -3, 0, -6, 0, 4, 1, 3, 0, 7, 6, -2, -5, 1, 0, 4, 3, 2, -1, 0, 0, -4, 0, 2, -1, 8, 2, 5, 0, 4, 7, -7, 0, 4, 0, -1, 4, 1, 5, 6, 0, 1, 6, 3, 0, 0, 0, 5, 1, 12, 0, -1, -1, 4, 5, 0, -1, 10, 0, 8, 13, -6, 0, -10, 0, -1, -9, -3, 0, 4, 0, -7, -10 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,4
LINKS
FORMULA
a[n]=A007953[n]-A095402[n]
EXAMPLE
n=1000: A007953[1000]=1, prime set={2,5}, A095402[1000]=7, a[1000]=1-7=-6
MATHEMATICA
ffi[x_] :=Flatten[FactorInteger[x]] lf[x_] :=Length[FactorInteger[x]] ba[x_] :=Table[Part[ffi[x], 2*j-1], {j, 1, lf[x]}] sd[x_] :=Apply[Plus, IntegerDigits[x]] tdp[x_] :=Flatten[Table[IntegerDigits[Part[ba[x], j]], {j, 1, lf[x]}], 1] sdp[x_] :=Apply[Plus, tdp[x]] a=Table[sd[w], {w, 1, 150}]; b=Table[sdp[w], {w, 1, 150}]; b-a
CROSSREFS
Sequence in context: A335022 A249732 A076694 * A011328 A048277 A059419
KEYWORD
base,sign
AUTHOR
Labos Elemer, Jun 21 2004
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 02:16 EDT 2024. Contains 371264 sequences. (Running on oeis4.)