cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A095403 Sum of digits of n minus the sum of digits of all distinct prime factors of n.

This page as a plain text file.
%I A095403 #4 Oct 15 2013 22:32:25
%S A095403 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,
%T A095403 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,
%U A095403 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
%N A095403 Sum of digits of n minus the sum of digits of all distinct prime factors of n.
%F A095403 a[n]=A007953[n]-A095402[n]
%e A095403 n=1000: A007953[1000]=1, prime set={2,5}, A095402[1000]=7, a[1000]=1-7=-6
%t A095403 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
%Y A095403 Cf. A007953, A051351.
%K A095403 base,sign
%O A095403 1,4
%A A095403 _Labos Elemer_, Jun 21 2004