login
Decimal expansion of Sum_{n >= 1, b(n) != 1} 1/b(n)^n, where b(n) = A217626(n).
2

%I #60 Jan 17 2018 17:51:02

%S 1,3,8,1,5,9,8,0,6,6,1,0,6,1,3,8,7,3,5,1,4,8,7,0,9,0,0,1,0,6,3,0,8,6,

%T 8,3,0,4,4,8,5,5,2,7,7,2,4,4,3,5,1,0,6,6,5,9,7,7,6,6,9,3,8,2,3,0,0,7,

%U 6,5,3,1,7,8,4,9,6,1,7,3,6,5,5,6,6,9,9,2,7,4,9,4,9,2,6,6,8,1,5,0,5,1,8,5,1,6,0,3

%N Decimal expansion of Sum_{n >= 1, b(n) != 1} 1/b(n)^n, where b(n) = A217626(n).

%C Decimal expansion of the limit to which tends the difference between the series for 1/A217626(x) and the Gamma function evaluated for M, where 2<M and 1<=x<M!. - _R. J. Cano_, May 25 2017

%H R. J. Cano, <a href="/A219995/b219995.txt">Table of n, a(n) for n = 0..4999</a>

%H R. J. Cano, <a href="/A219995/a219995_2.txt">Sequencer program.</a>

%e 0.138159806610613873514870900106308683044855277244...

%t First@ RealDigits@ N[#, 108] &@ Total@Map[1/#1^#2 & @@ # &, DeleteCases[ MapIndexed[{#1, First@ #2} &, #], _?(First@ # == 1 &)]] &@ Differences@ Flatten@ Table[(Drop[#, If[m == 1, 0, (m - 1)!]] - First@ #)/9 &@ Map[FromDigits, Permutations@ Range@ m], {m, 7}] (* _Michael De Vlieger_, May 25 2017, after _Jean-François Alcover_ at A215940 *)

%o (PARI) See LINKS.

%Y Cf. A215940, A217626.

%K nonn,base,easy,cons

%O 0,2

%A _R. J. Cano_, Dec 02 2012

%E Edited by _N. J. A. Sloane_, Dec 03 2012

%E Offset corrected by _Rick L. Shepherd_, Jan 14 2014