|
%I #17 Nov 30 2019 09:06:05
%S 1,1,2,4,6,6,6,8,90,10,6,12,210,14,30,16,30,18,42,20,770,22,6,24,
%T 13650,26,54,28,30,30,462,32,5610,34,210,36,51870,38,26,40,330,42,42,
%U 44,2070,46,6,48,324870,50,1122,52,30,54,43890,56,5510,58,6,60,930930
%N Denominator of ( A164555(n)/A027642(n) + 1/(n+1) ).
%C The sequence A174341(n)/a(n) = 2, 1, 1/2, 1/4, 1/6, 1/6, 1/6, ... becomes 2, -1, 1/2, -1/4, 1/6,.. under inverse binomial transform: an autosequence, where each second term flips the sign.
%H OEIS Wiki, <a href="https://oeis.org/wiki/Autosequence">Autosequence</a>
%o (PARI)
%o B(n)=if(n!=1, bernfrac(n), -bernfrac(n));
%o a(n)=denominator(B(n) + 1/(n + 1));
%o for(n=0, 60, print1(a(n),", ")) \\ _Indranil Ghosh_, Jun 19 2017
%o (Python)
%o from sympy import bernoulli, Rational
%o def B(n):
%o return bernoulli(n) if n != 1 else -bernoulli(n)
%o def a(n):
%o return (B(n) + Rational(1, n + 1)).as_numer_denom()[1]
%o [a(n) for n in range(61)] # _Indranil Ghosh_, Jun 19 2017
%Y Cf. A174341 (numerators).
%K nonn,frac
%O 0,3
%A _Paul Curtz_, Mar 16 2010
|
