# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/

Search: id:a083687
Showing 1-1 of 1

%I A083687 #24 Dec 03 2019 03:54:21
%S A083687 1,5,7,761,671,4572347,1171733,518413759,32956355893,1949885751497,
%T A083687 21495895979,63715389517501781,22630025105469577,36899945775958445129,
%U A083687 517210776697519633301437,4518133367201930332907311663
%N A083687 Numerator of B(2n)*H(2n)/n*(-1)^(n+1) where B(k) is the k-th Bernoulli number and H(k) the k-th harmonic number.
%H A083687 Seiichi Manyama, <a href="/A083687/b083687.txt">Table of n, a(n) for n = 1..260</a>
%H A083687 Ira Gessel, <a href="https://doi.org/10.1016/j.jnt.2003.08.010">On Miki's identity for Bernoulli numbers</a> J. Number Theory 110 (2005), no. 1, 75-82.
%F A083687 Miki's identity : B(n)*H(n)*(2/n) = sum(i=2, n-2, B(i)/i*B(n-i)/(n-i)*(1-C(n, i)))
%t A083687 Table[ BernoulliB[2n] * HarmonicNumber[2n] / n // Numerator // Abs, {n, 1, 16}] (* _Jean-François Alcover_, Mar 24 2015 *)
%o A083687 (PARI) a(n)=numerator((-1)^(n+1)*bernfrac(2*n)*sum(k=1,2*n,1/k)/n)
%o A083687 (Python)
%o A083687 from sympy import bernoulli, harmonic, numer
%o A083687 def a(n):
%o A083687     return numer(bernoulli(2 * n) * harmonic(2 * n) * (-1)**(n + 1) / n)
%o A083687 [a(n) for n in range(1, 31)]  # _Indranil Ghosh_, Aug 04 2017
%Y A083687 Cf. A083688.
%K A083687 frac,nonn
%O A083687 1,2
%A A083687 _Benoit Cloitre_, Jun 15 2003

# Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE