|
%I #13 Jan 26 2022 17:58:30
%S 0,1,1,1,1,1,1,1,1,1,1,2,3,4,5,6,7,8,9,10,1,3,1,5,3,7,2,9,5,11,1,4,5,
%T 2,7,8,1,10,11,4,1,5,3,7,1,9,5,11,3,13,1,6,7,8,9,2,11,12,13,14,1,7,2,
%U 1,5,11,1,13,7,5,1,8,9,10,11,12,13,2,15,16
%N Numerator of sum of reciprocals of all nonzero digits of n in decimal representation.
%H Reinhard Zumkeller, <a href="/A214949/b214949.txt">Table of n, a(n) for n = 0..10000</a>
%F a(A037264(n)) = a(A037268(n)) = a(A214958(n)) = a(A214959(n)) = 1;
%F a(n) = a(A004719(n)).
%t nsr[n_] := Numerator[Total[1/Select[IntegerDigits[n], # > 0 &]]]; nsr /@ Range[0, 79] (* _Jayanta Basu_, Jul 13 2013 *)
%o (Haskell)
%o import Data.Ratio ((%), numerator)
%o a214949 = f 0 where
%o f y 0 = numerator y
%o f y x = f (y + if d == 0 then 0 else 1 % d) x'
%o where (x',d) = divMod x 10
%o (PARI) a(n) = my(d=digits(n)); numerator(sum(k=1, #d, if (d[k], 1/d[k]))); \\ _Michel Marcus_, Jan 26 2022
%Y Cf. A004719, A037264, A037268, A214958, A214959.
%Y Cf. A214950 (denominators).
%K nonn,base,frac
%O 0,12
%A _Reinhard Zumkeller_, Aug 02 2012
|
