%I #23 Mar 07 2020 06:19:34
%S 2,12,12,120,20,504,168,720,180,1320,88,65520,10920,5040,720,24480,68,
%T 28728,3192,39600,27720,182160,1840,1965600,163800,39312,3024,97440,
%U 2320,3437280,229152,3769920,235620,42840,280,138181680,219336,35568,1872,3247200
%N Denominators of expansion of 1/x+1/log(1-x).
%C The numerator sequence is |A006232(n+1)|, n>=0.
%C |A006232(n+1)|= numerator(r(n)), n>=1, with r(n) := sum(|stirling1(n,k)|*B(k+1)/(k+1),k=1..n), n>=1 and B(n): =A027641(n)/A027642(n) (Bernoulli numbers) and stirling1(n,m)=A008275(n,m), n>=m>=1; r(0) := 1/2.
%F Denominators from e.g.f. 1/x + 1/log(1-x) (and of signed sequence from e.g.f. 1/x - 1/log(1+x)).
%F a(n) = denominator(r(n)), n>=0, with rational r(n) defined in one of the comments.
%e r(n) sequence, n>=0: 1/2, 1/12, 1/12, 19/120, 9/20, 863/504, 1375/168, 33953/720, 57281/180,...
%t With[{nn=40},Denominator[CoefficientList[Series[1/x+1/Log[1-x], {x,0,nn}] ,x] Range[0,nn]!]] (* _Harvey P. Dale_, Feb 18 2012 *)
%o (Sage)
%o def A075178_list(len):
%o f, R, C = 1, [0], [1]+[0]*len
%o for n in (1..len):
%o for k in range(n, 0, -1):
%o C[k] = -C[k-1] * k / (k + 1)
%o C[0] = -sum(C[k] for k in (1..n))
%o R.append((C[0]*f).denominator())
%o f *= n
%o return R[1:]
%o print(A075178_list(40)) # _Peter Luschny_, Feb 21 2016
%Y Cf. A006232, A075179, A006232 with A006233, A008275.
%K nonn,frac,easy
%O 0,1
%A _Wolfdieter Lang_, Sep 06, 2002