|
|
A000587
|
|
Rao Uppuluri-Carpenter numbers (or complementary Bell numbers): e.g.f. = exp(1 - exp(x)).
(Formerly M1913 N0755)
|
|
122
|
|
|
1, -1, 0, 1, 1, -2, -9, -9, 50, 267, 413, -2180, -17731, -50533, 110176, 1966797, 9938669, 8638718, -278475061, -2540956509, -9816860358, 27172288399, 725503033401, 5592543175252, 15823587507881, -168392610536153, -2848115497132448, -20819319685262839
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,6
|
|
COMMENTS
|
Alternating row sums of Stirling2 triangle A048993.
Row sums of triangle A144185 = a signed, shifted version of A000587: (1, 0, 1, 1, 2, 9, -9, 50, -267, -413, -2180, ...). Triangle A144185 is generated from A118433, the self-inverse triangle. - Gary W. Adamson, Sep 13 2008
Closely linked to A000110 and especially the contribution there of Jonathan R. Love (japanada11(AT)yahoo.ca), Feb 22 2007, by offering what is a complementary finding.
Number of set partitions of 1..n with an even number of parts, minus the number of such partitions with an odd number of parts. - Franklin T. Adams-Watters, May 04 2010
After -2, the smallest prime is a(36) = -1454252568471818731501051, no others through a(100). What is the first prime >0 in the sequence? - Jonathan Vos Post, Feb 02 2011
a(723) ~ 1.9*10^1265 is almost certainly prime. - D. S. McNeil, Feb 02 2011
Stirling transform of a(n) = [1, -1, 0, 1, 1, ...] is A033999(n) = [1, -1, 1, -1, 1, ...]. - Michael Somos, Mar 28 2012
Negated coefficients in the asymptotic expansion: A005165(n)/n! ~ 1 - 1/n + 1/n^2 + 0/n^3 - 1/n^4 - 1/n^5 + 2/n^6 + 9/n^7 + 9/n^8 - 50/n^9 - 267/n^10 - 413/n^11 + O(1/n^12), starting from the O(1/n) term. - Vladimir Reshetnikov, Nov 09 2016
Named after Venkata Ramamohana Rao Uppuluri and John A. Carpenter of the Mathematics Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee. They are called "Rényi numbers" by Fekete (1999), after the Hungarian mathematician Alfréd Rényi (1921-1970). - Amiram Eldar, Mar 11 2022
|
|
REFERENCES
|
N. A. Kolokolnikova, Relations between sums of certain special numbers (Russian), in Asymptotic and enumeration problems of combinatorial analysis, pp. 117-124, Krasnojarsk. Gos. Univ., Krasnoyarsk, 1976.
Alfréd Rényi, Új modszerek es eredmenyek a kombinatorikus analfzisben. I. MTA III Oszt. Ivozl., Vol. 16 (1966), pp. 7-105.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
M. V. Subbarao and A. Verma, Some remarks on a product expansion. An unexplored partition function, in Symbolic Computation, Number Theory, Special Functions, Physics and Combinatorics (Gainesville, FL, 1999), pp. 267-283, Kluwer, Dordrecht, 2001.
|
|
LINKS
|
W. Asakly, A. Blecher, C. Brennan, A. Knopfmacher, T. Mansour, and S. Wagner, Set partition asymptotics and a conjecture of Gould and Quaintance, Journal of Mathematical Analysis and Applications, Volume 416, Issue 2 (15 August 2014), Pages 672-682.
Valerio De Angelis and Dominic Marcello, Wilf's Conjecture, The American Mathematical Monthly, Vol. 123, No. 6 (2016), pp. 557-573.
Peter J. Larcombe, Jack Sutton, and James Stanton, A note on the constant 1/e, Palest. J. Math. (2023) Vol. 12, No. 2, 609-619. See p. 617.
|
|
FORMULA
|
E.g.f.: exp(1 - e^x).
a(n) = Sum_{k=0..n} (-1)^k S2(n, k), where S2(i, j) are the Stirling numbers of second kind A008277.
G.f.: (x/(1-x))*A(x/(1-x)) = 1 - A(x); the binomial transform equals the negative of the sequence shifted one place left. - Paul D. Hanna, Dec 08 2003
With different signs: g.f.: Sum_{k>=0} x^k/Product_{L=1..k} (1 + L*x).
Recurrence: a(n) = -Sum_{i=0..n-1} a(i)*C(n-1, i). - Ralf Stephan, Feb 24 2005
Let P be the lower-triangular Pascal-matrix, PE = exp(P-I) a matrix-exponential in exact integer arithmetic (or PE = lim exp(P)/exp(1) as limit of the exponential); then a(n) = PE^-1 [n,1]. - Gottfried Helms, Apr 08 2007
Take the series 0^n/0! - 1^n/1! + 2^n/2! - 3^n/3! + 4^n/4! + ... If n=0 then the result will be 1/e, where e = 2.718281828... If n=1, the result will be -1/e. If n=2, the result will be 0 (i.e., 0/e). As we continue for higher natural number values of n sequence for the Roa Uppuluri-Carpenter numbers is generated in the numerator, i.e., 1/e, -1/e, 0/e, 1/e, 1/e, -2/e, -9/e, -9/e, 50/e, 267/e, ... . - Peter Collins (pcolins(AT)eircom.net), Jun 04 2007
The sequence (-1)^n*a(n), with general term Sum_{k=0..n} (-1)^(n-k)*S2(n, k), has e.g.f. exp(1-exp(-x)). It also has Hankel transform (-1)^C(n+1,2)*A000178(n) and binomial transform A109747. - Paul Barry, Mar 31 2008
G.f.: 1 / (1 + x / (1 - x / (1 + x / (1 - 2*x / (1 + x / (1 - 3*x / (1 + x / ...))))))). - Michael Somos, May 12 2012
G.f.: -1/U(0) where U(k) = x*k - 1 - x + x^2*(k+1)/U(k+1); (continued fraction). - Sergei N. Gladkovskii, Sep 28 2012
G.f.: 1/(U(0)+x) where U(k) = 1 + x - x*(k+1)/(1 + x/U(k+1)); (continued fraction). - Sergei N. Gladkovskii, Oct 12 2012
G.f.: 1+x/G(0) where G(k) = x*k - 1 + x^2*(k+1)/G(k+1); (continued fraction). - Sergei N. Gladkovskii, Nov 04 2012
G.f.: (1 - G(0))/(x+1) where G(k) = 1 - 1/(1-k*x)/(1-x/(x+1/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Jan 17 2013
G.f.: 1 + x/(G(0)-x) where G(k) = x*k + 2*x - 1 - x*(x*k+x-1)/G(k+1); (continued fraction). - Sergei N. Gladkovskii, Jan 28 2013
G.f.: G(0)/(1+x), where G(k) = 1 - x^2*(k+1)/(x^2*(k+1) + (x*k - 1 - x)*(x*k - 1)/G(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Feb 07 2014
a(n) = Sum_{k=0..n} (-1)^k*Bell(k)*A129062(n, k).
a(n) = Sum_{k=0..n} (-1)^k*k!*A130191(n, k). (End)
|
|
EXAMPLE
|
G.f. = 1 - x + x^3 + x^4 - 2*x^5 - 9*x^6 - 9*x^7 + 50*x^8 + 267*x^9 + 413*x^10 - ...
|
|
MAPLE
|
# Uses floating point arithmetic, increase working precision for n>60.
if n = 0 then 1 else a := [seq(2, i=1..n-1)]; b := [seq(1, i=1..n-1)];
-exp(-x)*hypergeom(a, b, x); round(evalf(subs(x=-1, %), 66)) fi end:
# second Maple program:
b:= proc(n, t) option remember; `if`(n=0, 1-2*t,
add(b(n-j, 1-t)*binomial(n-1, j-1), j=1..n))
end:
a:= n-> b(n, 0):
|
|
MATHEMATICA
|
Table[ -1 * Sum[ (-1)^( k + 1) StirlingS2[ n, k ], {k, 0, n} ], {n, 0, 40} ]
With[{nn=30}, CoefficientList[Series[Exp[1-Exp[x]], {x, 0, nn}], x] Range[0, nn]!] (* Harvey P. Dale, Nov 04 2011 *)
a[ n_] := If[ n < 0, 0, n! SeriesCoefficient[ Exp[ 1 - Exp[x]], {x, 0, n}]]; (* Michael Somos, May 27 2014 *)
a[ n_] := If[ n < 0, 0, With[{m = n + 1}, SeriesCoefficient[ Series[ Nest[ x Factor[ 1 - # /. x -> x / (1 - x)] &, 0, m], {x, 0, m}], {x, 0, m}]]]; (* Michael Somos, May 27 2014 *)
b[1] = 1; k = 1; Flatten[{1, Table[Do[j = k; k -= b[m]; b[m] = j; , {m, 1, n-1}]; b[n] = k; k*(-1)^n, {n, 1, 40}]}] (* Vaclav Kotesovec, Sep 09 2019 *)
|
|
PROG
|
(PARI) {a(n) = if( n<0, 0, n! * polcoeff( exp( 1 - exp( x + x * O(x^n))), n))}; /* Michael Somos, Mar 14 2011 */
(PARI) {a(n) = local(A); if( n<0, 0, n++; A = O(x); for( k=1, n, A = x - x * subst(A, x, x / (1 - x))); polcoeff( A, n))}; /* Michael Somos, Mar 14 2011 */
(PARI) Vec(serlaplace(exp(1 - exp(x+O(x^99))))) /* Joerg Arndt, Apr 01 2011 */
(PARI) a(n)=round(exp(1)*suminf(k=0, (-1)^k*k^n/k!))
vector(20, n, a(n-1)) \\ Derek Orr, Sep 19 2014 -- a direct approach
(PARI) x='x+O('x^66); Vec(serlaplace(exp(1 - exp(x)))) \\ Michel Marcus, Sep 19 2014
(Python) # The objective of this implementation is efficiency.
# n -> [a(0), a(1), ..., a(n)] for n > 0.
A = [0 for i in range(n)]
A[n-1] = 1
R = [1]
for j in range(0, n):
A[n-1-j] = -A[n-1]
for k in range(n-j, n):
A[k] += A[k-1]
R.append(A[n-1])
return R
(Python)
# Python 3.2 or higher required
from itertools import accumulate
A000587, blist, b = [1, -1], [1], -1
for _ in range(30):
blist = list(accumulate([b]+blist))
b = -blist[-1]
(Haskell)
a000587 n = a000587_list !! n
a000587_list = 1 : f a007318_tabl [1] where
f (bs:bss) xs = y : f bss (y : xs) where y = - sum (zipWith (*) xs bs)
|
|
CROSSREFS
|
|
|
KEYWORD
|
sign,easy,nice,changed
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|