cp's OEIS Frontend

A225480 a(n) = B2(n) * C(n) where B2(n) are generalized Bernoulli numbers and C(n) the Clausen numbers.

%I A225480 #17 Jan 15 2025 19:38:44
%S A225480 1,0,-2,0,14,0,-62,0,254,0,-5110,0,2828954,0,-114674,0,237036478,0,
%T A225480 -11499383114,0,183092554714,0,-3584085584926,0,3965530936622474,0,
%U A225480 -573989008898786,0,6375197353574922166,0,-9251189109760413581110,0,33111281730973040956798,0
%N A225480 a(n) = B2(n) * C(n) where B2(n) are generalized Bernoulli numbers and C(n) the Clausen numbers.
%C A225480 The Clausen numbers C(n) are T(n, 1) in A160014.
%H A225480 Peter Luschny, <a href="http://www.luschny.de/math/euler/StirlingFrobeniusNumbers.html">Stirling-Frobenius numbers</a>
%H A225480 Peter Luschny, <a href="http://www.luschny.de/math/euler/GeneralizedBernoulliNumbers.html">Generalized Bernoulli numbers</a>.
%F A225480 Let B(n,m) = sum_{k = 0..n} sum_{j = 0..k} sum_{v = 0..j} ((-1)^(n-v)/(j+1)) *binomial(n,k)*binomial(j,v)*(m*v)^k then a(n) = B(n,2)*A141056(n).
%F A225480 Let B2(n) = sum_{k=0..n} ((-1)^k*k!/(k+1)) S_{2}(n, k) where S_{2}(n, k) the Stirling-Frobenius subset numbers A039755(n, k) then a(n) = B2(n)*A141056(n).
%e A225480 The numerators of 1/1, 0/2, -2/6, 0/2, 14/30, 0/2, -62/42, 0/2, 254/30, 0/2, -5110/66, 0/2, 2828954/2730, ... (the denominators are the Clausen numbers).
%p A225480 B := (n, m) -> add(add(add(((-1)^(n-v)/(j+1))*binomial(n,k)*binomial(j, v)*(m*v)^k, v = 0..j), j = 0..k), k = 0..n);
%p A225480 C := proc(n) numtheory[divisors](n);map(i->i+1,%);select(isprime,%);mul(i,i=%) end:
%p A225480 A225480 := n -> B(n, 2)*C(n); seq(A225480(n), n = 0..33);
%t A225480 B[n_, m_] := Sum[((-1)^(n - v)/(j + 1))*Binomial[n, k]*Binomial[j, v]*If[k == 0, 1, (m*v)^k], {k, 0, n}, {j, 0, k}, {v, 0, j}];
%t A225480 c[n_] := Denominator[Sum[Boole[PrimeQ[d + 1]]/(d + 1), {d, Divisors[n]}]];
%t A225480 a[n_] := B[n, 2]*c[n];
%t A225480 Table[a[n], {n, 0, 33}] (* _Jean-François Alcover_, Aug 02 2019, from Maple *)
%o A225480 (Sage)
%o A225480 @CachedFunction
%o A225480 def EulerianNumber(n, k, m) :   # The Eulerian numbers
%o A225480     if n == 0: return 1 if k == 0 else 0
%o A225480     return (m*(n-k)+m-1)*EulerianNumber(n-1,k-1,m)+(m*k+1)*EulerianNumber(n-1,k,m)
%o A225480 @CachedFunction
%o A225480 def B(n, m):   # The generalized Bernoulli numbers
%o A225480     return add(add(EulerianNumber(n, j, m)*binomial(j, n - k)
%o A225480            for j in (0..n))*(-1)^k/(k+1) for k in (0..n))
%o A225480 def A225480(n):
%o A225480     if n == 0: return 1
%o A225480     C = mul(filter(lambda s: is_prime(s) , map(lambda i: i+1, divisors(n))))
%o A225480     return C*B(n, 2)
%o A225480 print([A225480(n) for n in (0..33)])
%Y A225480 Cf. A141056, A160014.
%K A225480 sign,frac
%O A225480 0,3
%A A225480 _Peter Luschny_, May 30 2013