cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-5 of 5 results.

A000629 Number of necklaces of partitions of n+1 labeled beads.

Original entry on oeis.org

1, 2, 6, 26, 150, 1082, 9366, 94586, 1091670, 14174522, 204495126, 3245265146, 56183135190, 1053716696762, 21282685940886, 460566381955706, 10631309363962710, 260741534058271802, 6771069326513690646, 185603174638656822266, 5355375592488768406230
Offset: 0

Views

Author

N. J. A. Sloane, Don Knuth, Nick Singer (nsinger(AT)eos.hitc.com)

Keywords

Comments

Also the number of logically distinct strings of first order quantifiers in which n variables occur (C. S. Peirce, c. 1903). - Stephen Pollard (spollard(AT)truman.edu), Jun 07 2002
Stirling transform of A052849(n) = [2, 4, 12, 48, 240, ...] is a(n) = [2, 6, 26, 150, 1082, ...]. - Michael Somos, Mar 04 2004
Stirling transform of A000142(n-1) = [1, 1, 2, 6, 24, ...] is a(n-1) = [1, 2, 6, 26, ...]. - Michael Somos, Mar 04 2004
Stirling transform of (-1)^n * A024167(n-1) = [0, 1, -1, 5, -14, 94, ...] is a(n-2) = [0, 1, 2, 6, 26, ...]. - Michael Somos, Mar 04 2004
The asymptotic expansion of 2*log(n) - (2^1*log(1) + 2^2*log(2) + ... + 2^n*log(n))/2^n is (a(1)/1)/n + (a(2)/2)/n^2 + (a(3)/3)/n^3 + ... - Michael Somos, Aug 22 2004
This is the sequence of cumulants of the probability distribution of the number of tails before the first head in a sequence of fair coin tosses. - Michael Hardy (hardy(AT)math.umn.edu), May 01 2005
Appears to be row sums of A154921. - Mats Granvik, Jan 18 2009
This is the number of cyclically ordered partitions of n+1 labeled points. The ordered version is A000670. - Michael Somos, Jan 08 2011
A000670(n+1) = p(n+1) where p(x) is the unique degree-n polynomial such that p(k) = a(k) for k = 0, 1, ..., n. - Michael Somos, Apr 27 2012
Row sums of A154921 as conjectured above by Granvik. a(n) gives the number of outcomes of a race between n horses H1,...,Hn, where if a horse falls it is not ranked. For example, when n = 2 the 6 outcomes are a dead heat, H1 wins H2 second, H2 wins H1 second, H1 wins H2 falls, H2 wins H1 falls or both fall. - Peter Bala, May 15 2012
Also the number of disjoint areas of a Venn diagram for n multisets. - Aurelian Radoaca, Jun 27 2016
Also the number of ways of ordering n nonnegative integers, allowing for the possibility of ties, and also comparing the smallest integers with 0. Each comparison with 0 gives two possibilities, x > 0 or x=0. As such, without comparison with 0, we get A000670, the number of ways of ordering n nonnegative integers, allowing for the possibility of ties, or the number of ways n competitors can rank in a competition, allowing for the possibility of ties. For instance, for 2 nonnegative integers x,y, there are the following 6 ways of ordering them: x = y = 0, x = y > 0, x > y = 0, x > y > 0, y > x = 0, y > x > 0. - Aurelian Radoaca, Jul 09 2016
Also the number of ordered set partitions of subsets of {1,...,n}. Also the number of chains of distinct nonempty subsets of {1,...,n}. - Gus Wiseman, Feb 01 2019
Number of combinations of a Simplex lock having n buttons.
Row sums of the unsigned cumulant expansion polynomials A127671 and logarithmic polynomials A263634. - Tom Copeland, Jun 04 2021
Also the number of vertices in the axis-aligned polytope consisting of all vectors x in R^n where, for all k in {1,...,n}, the k-th smallest coordinate of x lies in the interval [0, k]. - Adam P. Goucher, Jan 18 2023
Number of idempotent Boolean relation matrices whose complement is also idempotent. See Rosenblatt link. - Geoffrey Critzer, Feb 26 2023

Examples

			a(2)=6: the necklace representatives on 1,2,3 are ({123}), ({12},{3}), ({13},{2}), ({23},{1}), ({1},{2},{3}), ({1},{3},{2})
G.f. = 1 + 2*x + 6*x^2 + 26*x^3 + 150*x^4 + 1082*x^5 + 9366*x^6 + 94586*x^7 + ...
From _Gus Wiseman_, Feb 01 2019: (Start)
The a(3) = 26 ordered set partitions of subsets of {1,2,3} are:
  {}  {{1}}  {{2}}  {{3}}  {{12}}    {{13}}    {{23}}    {{123}}
                           {{1}{2}}  {{1}{3}}  {{2}{3}}  {{1}{23}}
                           {{2}{1}}  {{3}{1}}  {{3}{2}}  {{12}{3}}
                                                         {{13}{2}}
                                                         {{2}{13}}
                                                         {{23}{1}}
                                                         {{3}{12}}
                                                         {{1}{2}{3}}
                                                         {{1}{3}{2}}
                                                         {{2}{1}{3}}
                                                         {{2}{3}{1}}
                                                         {{3}{1}{2}}
                                                         {{3}{2}{1}}
(End)
		

References

  • R. Austin, R. K. Guy, and R. Nowakowski, unpublished notes, circa 1987.
  • N. G. de Bruijn, Asymptotic Methods in Analysis, Dover, 1981, p. 36.
  • Eric Hammer, The Calculations of Peirce's 4.453, Transactions of the Charles S. Peirce Society, Vol. 31 (1995), pp. 829-839.
  • D. E. Knuth, personal communication.
  • J. D. E. Konhauser et al., Which Way Did the Bicycle Go?, MAA 1996, p. 174.
  • Charles Sanders Peirce, Collected Papers, eds. C. Hartshorne and P. Weiss, Harvard University Press, Cambridge, Vol. 4, 1933, pp. 364-365. (CP 4.453 in the electronic edition of The Collected Papers of Charles Sanders Peirce.)
  • Dawidson Razafimahatolotra, Number of Preorders to Compute Probability of Conflict of an Unstable Effectivity Function, Preprint, Paris School of Economics, University of Paris I, Nov 23 2007.

Crossrefs

Same as A076726 except for a(0). Cf. A008965, A052861, A008277.
Binomial transform of A000670, also double of A000670. - Joe Keane (jgk(AT)jgk.org)
A002050(n) = a(n) - 1.
A000629, A000670, A002050, A052856, A076726 are all more-or-less the same sequence. - N. J. A. Sloane, Jul 04 2012
Row sums of A028246.
A diagonal of the triangular array in A241168.
Row sums of unsigned A127671 and A263634.

Programs

  • Maple
    spec := [ B, {B=Cycle(Set(Z,card>=1))}, labeled ]; [seq(combstruct[count](spec, size=n), n=0..20)];
    a:=n->add(Stirling2(n+1,k)*(k-1)!,k=1..n+1); # Mike Zabrocki, Feb 05 2005
  • Mathematica
    a[ 0 ] = 1; a[ n_ ] := (a[ n ] = 1 + Sum[ Binomial[ n, k ] a[ n-k ], {k, 1, n} ])
    Table[ PolyLog[n, 1/2], {n, 0, -18, -1}] (* Robert G. Wilson v, Aug 05 2010 *)
    a[ n_] := If[ n<0, 0, PolyLog[ -n, 1/2]]; (* Michael Somos, Mar 07 2011 *)
    Table[Sum[(-1)^(n-k) StirlingS2[n,k]k! 2^k,{k,0,n}],{n,0,20}] (* Harvey P. Dale, Oct 21 2011 *)
    Join[{1}, Rest[t=30; Range[0, t]! CoefficientList[Series[2/(2 - Exp[x]), {x, 0, t}], x]]] (* Vincenzo Librandi, Jan 02 2016 *)
  • PARI
    {a(n) = if( n<0, 0, n! * polcoeff(subst( (1 + y) / (1 - y), y, exp(x + x * O(x^n)) - 1), n))} /* Michael Somos, Mar 04 2004 */
    
  • PARI
    {a(n)=polcoeff(sum(m=0, n, 2^m*m!*x^m/prod(k=1, m, 1+k*x+x*O(x^n))), n)} \\ Paul D. Hanna, Jul 20 2011
    
  • Python
    from math import comb
    from functools import lru_cache
    @lru_cache(maxsize=None)
    def A000629(n): return 1+sum(comb(n,j)*A000629(j) for j in range(n)) if n else 1 # Chai Wah Wu, Sep 25 2023

Formula

a(n) = 2*A000670(n) - 0^n. - Michael Somos, Jan 08 2011
O.g.f.: Sum_{n>=0} 2^n*n!*x^n / Product_{k=0..n} (1+k*x). - Paul D. Hanna, Jul 20 2011
E.g.f.: exp(x) / (2 - exp(x)) = d/dx log(1 / (2 - exp(x))).
a(n) = Sum_{k>=1} k^n/2^k.
a(n) = 1 + Sum_{j=0..n-1} C(n, j)*a(j).
a(n) = round(n!/log(2)^(n+1)) (just for n <= 15). - Henry Bottomley, Jul 04 2000
a(n) is asymptotic to n!/log(2)^(n+1). - Benoit Cloitre, Oct 20 2002
a(n) = Sum_{k=0..n} (-1)^(n-k)*Stirling2(n, k)*k!*2^k. - Vladeta Jovovic, Sep 29 2003
a(n) = Sum_{k=1..n} A008292(n, k)*2^k; A008292: triangle of Eulerian numbers. - Philippe Deléham, Jun 05 2004
a(1) = 1, a(n) = 2*Sum_{k=1..n-1} k!*A008277(n-1, k) for n>1 or a(n) = Sum_{k=1..n} (k-1)!*A008277(n, k). - Mike Zabrocki, Feb 05 2005
a(n) = Sum_{k=0..n} Stirling2(n+1, k+1)*k!. - Paul Barry, Apr 20 2005
A000629 = binomial transform of this sequence. a(n) = sum of terms in n-th row of A028246. - Gary W. Adamson, May 30 2005
a(n) = 2*(-1)^n * n!*Laguerre(n,P((.),2)), umbrally, where P(j,t) are the polynomials in A131758. - Tom Copeland, Sep 28 2007
a(n) = 2^n*A(n,1/2); A(n,x) the Eulerian polynomials. - Peter Luschny, Aug 03 2010
a(n) = (-1)^n*b(n), where b(n) = -2*Sum_{k=0..n-1} binomial(n,k)*b(k), b(0)=1. - Vladimir Kruchinin, Jan 29 2011
Row sums of A028246. Let f(x) = x+x^2. Then a(n+1) = (f(x)*d/dx)^n f(x) evaluated at x = 1. - Peter Bala, Oct 06 2011
O.g.f.: 1+2*x/(U(0)-2*x) where U(k)=1+3*x+3*x*k-2*x*(k+2)*(1+x+x*k)/U(k+1); (continued fraction). - Sergei N. Gladkovskii, Nov 14 2011
E.g.f.: exp(x)/(2 - exp(x)) = 2/(2-Q(0))-1; Q(k)=1+x/(2*k+1-x*(2*k+1)/(x+(2*k+2)/Q(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Nov 14 2011
G.f.: 1 / (1 - 2*x / (1 - 1*x / (1 - 4*x / (1 - 2*x / (1 - 6*x / ...))))). - Michael Somos, Apr 27 2012
PSUM transform of A162509. BINOMIAL transform is A007047. - Michael Somos, Apr 27 2012
G.f.: 1/G(0) where G(k) = 1 - x*(2*k+2)/( 1 - x*(k+1)/G(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Mar 23 2013
E.g.f.: 1/E(0) where E(k) = 1 - x/(k+1)/(1 - 1/(1 + 1/E(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Mar 27 2013
G.f.: T(0)/(1-2*x), where T(k) = 1 - 2*x^2*(k+1)^2/(2*x^2*(k+1)^2 - (1 - 2*x - 3*x*k)*(1 - 5*x - 3*x*k)/T(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Oct 29 2013
a(n) = log(2)*integral_{x>=0} (ceiling(x))^n * 2^(-x) dx. - Peter Bala, Feb 06 2015

Extensions

a(19) from Michael Somos, Mar 07 2011

A162508 A binomial sum of powers related to the Bernoulli numbers, triangular array, read by rows.

Original entry on oeis.org

-1, -2, 2, -4, 10, -6, -8, 38, -54, 24, -16, 130, -330, 336, -120, -32, 422, -1710, 3000, -2400, 720, -64, 1330, -8106, 21840, -29400, 19440, -5040, -128, 4118, -36414, 141624, -285600, 312480, -176400, 40320
Offset: 1

Views

Author

Peter Luschny, Jul 05 2009

Keywords

Comments

T(n,k) = sum_{v=0..k} (-1)^v*v*binomial(k,v)*(v+1)^(n-1)
for n >= 1, k >= 1; by convention T(0,0) = 1.
Gives a representation of the Bernoulli numbers B_{n} = B_{n}(1) (with B_1 = 1/2)
B_{n} = sum_{j=0..n} sum_{k=0..j} T(j,k)/(k+1)
T(n,1) = -2^(n-1) (n>=1)
T(n,n) = (-1)^n*n! (n>=1)
sum_{k=0..n} T(n,k) = -A000007(n-1) = -1,0,0,0,0,... (n>=1)
sum_{k=0..n} abs(T(n,k)) = A162509(n) = A073146(n,n-1) (n>=1)
sum_{k=0..n} T(n,k)/(k+1) = Bernoulli(n,1)-Bernoulli(n-1,1) (n>=1)
numer(sum(T(n,k)/(k+1),k=0..n)) = A051716(n) (n>=0)
denom(sum(T(n,k)/(k+1),k=0..n)) = A051717(n) (n>=0)
Contribution from Peter Luschny, Jul 08 2009: (Start)
More generally, define the polynomials (assume p[0,0](x)=1 and 0^0=1)
p[n,k](x) = sum_{v=0..k} (-1)^v*v*binomial(k,v)*(v+1+x)^(n-1)
[1], [0, -1], [0, -2-x, 2], [0, -4-4x-x^2, 10+4x, -6], ...
then T(n,k)=p[n,k](0) and (-1)^k*k!*Stirling2(n,k)=p[n,k](-1) (cf. A019538).
Assume now k >= 1 and read by rows. Then
p[n,k](1) = -1,-3,2,-9,14,-6,-27,74,-72,24,-81,350,-582,432,-120,...
(-1)^n*(-2)^(n-k)*p[n,k](-1/2))=1,3,2,9,16,6,27,98,90,24,81,544,924,576,120,..
(-1)^n*(-2)^(n-k)*p[n,k](-3/2))=1,1,2,1,8,6,1,26,54,24,1,80,348,384,120,... (End)
Variant of A199400.

Examples

			For n >= 1, k >= 1:
..................... -1
................... -2, 2
................. -4, 10, -6
.............. -8, 38, -54, 24
......... -16, 130, -330, 336, -120
..... -32, 422, -1710, 3000, -2400, 720
-64, 1330, -8106, 21840, -29400, 19440, -5040
		

Crossrefs

Programs

  • Maple
    T := proc(n,k) local v; if n=0 and k=0 then 1 else
    add((-1)^v*v*binomial(k,v)*(v+1)^(n-1),v=0..k) fi end:
    # Peter Bala's e.g.f. assuming offset 0:
    egf := (x, z) -> -((1-x)/exp(z) + x)^(-2):
    ser := series(egf(x, z), z, 10): coz := n -> n!*coeff(ser, z, n):
    row := n -> seq(coeff(coz(n), x, k), k = 0..n):
    seq(print(row(n)), n = 0..9); # Peter Luschny, Jan 28 2021
  • Mathematica
    t[n_, k_] := Sum[(-1)^v*v*Binomial[k, v]*(v + 1)^(n - 1), {v, 0, k}]; Table[t[n, k], {n, 1, 8}, {k, 1, n}] // Flatten (* Jean-François Alcover, Jul 26 2013, after Maple *)
  • Sage
    def A162508(n, k):
        if n==0 and k==0: return 1
        return add((-1)^v*v*binomial(k, v)*(v+1)^(n-1) for v in (0..k))
    for n in (1..8): [A162508(n, k) for k in (1..n)] # Peter Luschny, Jul 21 2014

Formula

From Peter Bala, Jul 21 2014: (Start)
T(n,k) = (-1)^k*k!*( Stirling2(n+1,k+1) - Stirling2(n,k+1) ), 1 <= k <= n.
T(n,k) = (-1)^k*(k + 1)*A038719(n+1,k+1).
E.g.f.: - B(-x,z)^2, where B(x,z) = 1/((1 + x)*exp(-z) - x) = 1 + (1 + x)*z + (1 + 3*x + 2*x^2)*z^2/2! + ... is an e.g.f. for A028246 (with an offset of 0).
Recurrence: T(n,k) = (k + 1)*T(n-1,k) - k*T(n-1,k-1).
The unsigned version of the triangle equals the matrix product A007318*A019538.
Assuming this triangle is a signed version of A199400 then the n-th row polynomial R(n,x) = 1/(1 - x)*( sum {k = 1..inf} k*(k + 1)^(n-1)*(x/(x - 1))^k ), valid for x in the open interval (-inf, 1/2). (End)

Extensions

More terms from Philippe Deléham, Nov 05 2011

A199335 Triangle T(n,k), read by rows, given by (0,1,0,2,0,3,0,4,0,5,0,6,0,7,0,8,0,9,...) DELTA (2,0,3,0,4,0,5,0,6,0,7,0,8,0,9,...), where DELTA is the operator defined in A084938.

Original entry on oeis.org

1, 0, 2, 0, 2, 4, 0, 2, 14, 8, 0, 2, 36, 66, 16, 0, 2, 82, 342, 262, 32, 0, 2, 176, 1436, 2416, 946, 64, 0, 2, 366, 5364, 16844, 14394, 3222, 128, 0, 2, 748, 18654, 99560, 156190, 76908, 10562, 256
Offset: 0

Views

Author

Philippe Deléham, Nov 05 2011

Keywords

Comments

Following an observation by Dale Gerdemann, it appears that T(n,k) = A120434(n+1,n-k) for n>=1, k>=1. - M. F. Hasler, Apr 18 2015
See also A144696. - Antti Karttunen, Apr 21 2015

Examples

			Triangle begins :
1
0, 2
0, 2, 4
0, 2, 14, 8
0, 2, 36, 66, 16
0, 2, 82, 342, 262, 32
0, 2, 176, 1436, 2416, 946, 64
		

Crossrefs

Formula

Sum_k{k, 0<=k<=n} T(n,k)*x^k = A000007(n), A000142(n+1), A162509(n+1) for x=0,1,2 respectively.
Sum_{k, 0<=k<=n} T(n,k)^2^(n-k) = A005649(n).

Extensions

Typo in 8th row corrected by Olivier Gérard, Oct 29 2012

A199400 Triangle T(n,k), read by rows, given by (2,0,3,0,4,0,5,0,6,0,7,0,8,0,9,...) DELTA (2,1,3,2,4,3,5,4,6,5,7,6,8,7,9,...) where DELTA is the operator defined in A084938.

Original entry on oeis.org

1, 2, 2, 4, 10, 6, 8, 38, 54, 24, 16, 130, 330, 336, 120, 32, 422, 1710, 3000, 2400, 720, 64, 1330, 8106, 21840, 29400, 19440, 5040, 128, 4118, 36414, 141624, 285600, 312480, 176400, 40320, 256, 12610, 158010, 853776, 2421720, 3900960, 3598560, 1774080, 362880
Offset: 0

Views

Author

Philippe Deléham, Nov 05 2011

Keywords

Comments

Variant of A162508.

Examples

			Triangle begins :
1
2, 2
4, 10, 6
8, 38, 54, 24
16, 130, 330, 336, 120
32, 422, 1710, 3000, 2400, 720
		

Crossrefs

Formula

T(n,k)=(k+1)!*A143494(n+2,k+2)
T(n,k)=(k+1)*T(n-1,k-1)+(k+2)*T(n-1,k).
Sum_{k, 0<=k<=n} T(n,k)=A162509(n+1).
T(n,n)=(n+1)!=A000142(n+1)
T(n,0)=2^n=A000079(n).

A299404 a(n) = 1 + Sum_{m >= 1} (m + 1)^n/2^(m - 1).

Original entry on oeis.org

3, 7, 23, 103, 599, 4327, 37463, 378343, 4366679, 56698087, 817980503, 12981060583, 224732540759, 4214866787047, 85130743763543, 1842265527822823, 42525237455850839, 1042966136233087207, 27084277306054762583, 742412698554627289063, 21421502369955073624919
Offset: 0

Views

Author

Joseph Wheat, Feb 20 2018

Keywords

Crossrefs

Programs

  • Mathematica
    Table[1 + LerchPhi[1/2, -n, 2], {n, 0, 20}] (* Vaclav Kotesovec, Apr 17 2018 *)
  • PARI
    a(n) = 1+ round(suminf(m=1, (m + 1)^n/2^(m - 1)));

Formula

a(n + 1) = 4*A162509(n + 1) + a(n).
a(n) = 2*A007047(n) + 1.
{a(4n - 3), a(4n - 2), a(4n - 1), a(4n)} mod 10 = {7, 3, 3, 9} for n > 0.
floor(log_2(a(n))) = A083652(n).
Lim_{n->infinity} (a(n)^(1/n))/n = 1/(e*log(2)). - Jon E. Schoenfield, Feb 24 2018
a(n)/n! ~ 4 / (log(2))^(n+1). - Vaclav Kotesovec, Apr 17 2018
Showing 1-5 of 5 results.