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

A002050 Number of simplices in barycentric subdivision of n-simplex.

Original entry on oeis.org

0, 1, 5, 25, 149, 1081, 9365, 94585, 1091669, 14174521, 204495125, 3245265145, 56183135189, 1053716696761, 21282685940885, 460566381955705, 10631309363962709, 260741534058271801, 6771069326513690645
Offset: 0

Views

Author

Keywords

Comments

Stirling transform of A052849(n)=[1,4,12,48,240,...] is a(n)=[1,5,25,149,1081,...]. - Michael Somos, Mar 04 2004
Stirling transform of A000142(n-1)=[0,1,2,6,24,...] is a(n-1)=[0,1,5,25,149,...]. - Michael Somos, Mar 04 2004
Stirling transform of 2*A005359(n-1)=[1,0,4,0,48,0,...] is a(n-1)=[1,1,5,25,149,...]. - Michael Somos, Mar 04 2004
"Stirling-Bernoulli transform" of A000225. - Paul Barry, Apr 20 2005
a(n) is the number of nonempty words that can be formed from an alphabet of nonempty subsets of [n] so that the letters in each word are pairwise disjoint. - Geoffrey Critzer, Apr 12 2009
Row sums of A053440. - Peter Bala, Jul 12 2014
Conjecture: Let k be a positive integer. The sequence obtained by reducing a(n) modulo k is eventually periodic with the period dividing phi(k) = A000010(k). For example, modulo 9 we obtain the sequence [0, 1, 5, 7, 5, 1, 5, 4, 5, 7, 5, 1, 5, 4, 5, 7, 5, 1, 5, 4, 5, 7, 5, ...], with an apparent period of 6 = phi(9) beginning at a(5). - Peter Bala, Aug 03 2023

References

  • R. Austin, R. K. Guy, and R. Nowakowski, unpublished notes, circa 1987.
  • Ulrike Sattler, Decidable classes of formal power series with nice closure properties, Diplomarbeit im Fach Informatik, Univ. Erlangen - Nuernberg, Jul 27 1994
  • 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).

Crossrefs

A000629, A000670, A002050, A052856, A076726 are all more-or-less the same sequence. - N. J. A. Sloane, Jul 04 2012
A diagonal of the triangle in A241168. Row sums of A053440.

Programs

  • Mathematica
    Table[Sum[Binomial[n, i]*Sum[StirlingS2[i, k]*k!, {k, 1, i}], {i, 1, n}], {n, 0, 20}] (* Geoffrey Critzer, Apr 12 2009 *)
    With[{nn=20},CoefficientList[Series[(Exp[2x]-Exp[x])/(2-Exp[x]),{x,0,nn}],x] Range[0,nn]!] (* Harvey P. Dale, May 28 2013 *)
    a[0] = 0; a[n_] := 2*Sum[k!*StirlingS2[n, k], {k, 2, n}] + 1; Table[a[n], {n, 0, 18}] (* Jean-François Alcover, Sep 27 2013, after Vladimir Kruchinin *)
  • PARI
    a(n)=if(n<0,0,n!*polcoeff(subst((y+y^2)/(1-y),y,exp(x+x*O(x^n))-1),n));

Formula

E.g.f.: (exp(2x)-exp(x))/(2-exp(x)).
a(n) = A000629(n) - 1.
a(n) = Sum_{k=0..n} (-1)^(n-k)k!*S2(n, k)(2^k-1). - Paul Barry, Apr 20 2005
a(n) = Sum_{k=1...n} binomial(n,k)*A000670(k). - Geoffrey Critzer, Apr 12 2009
a(n) ~ n!/log(2)^(n+1). - Vaclav Kotesovec, Jul 29 2013
a(n) = 1 + 2*Sum_{k=2..n} k!*Stirling2(n,k), n > 0, a(0)=1. - Vladimir Kruchinin, Sep 27 2013
G.f.: T(0)/(1-2*x) - 1/(1-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, Nov 29 2013
G.f.: Sum_{j>=1} j!*x^j / Product_{k=0..j} (1 - (k + 1)*x). - Ilya Gutkovskiy, Apr 04 2019

Extensions

More terms from James Sellers, Aug 22 2000

A002051 Steffensen's bracket function [n,2].

Original entry on oeis.org

0, 0, 1, 9, 67, 525, 4651, 47229, 545707, 7087005, 102247051, 1622631549, 28091565547, 526858344285, 10641342962251, 230283190961469, 5315654681948587, 130370767029070365, 3385534663256714251, 92801587319328148989, 2677687796244383678827, 81124824998504072833245, 2574844419803190382447051
Offset: 1

Views

Author

Keywords

Comments

a(n) is the number of ways to arrange the blocks of the partitions of {1,2,...,n} in an undirected cycle of length 3 or more, see A000629. - Geoffrey Critzer, Nov 23 2012
From Gus Wiseman, Jun 24 2020: (Start)
Also the number of (1,2)-matching length-n sequences covering an initial interval of positive integers. For example, the a(2) = 1 and a(3) = 9 sequences are:
(1,2) (1,1,2)
(1,2,1)
(1,2,2)
(1,2,3)
(1,3,2)
(2,1,2)
(2,1,3)
(2,3,1)
(3,1,2)
Missing from this list are:
(1,1) (1,1,1)
(2,1) (2,1,1)
(2,2,1)
(3,2,1)
(End)

Examples

			a(4) = 9. There are 6 partitions of {1,2,3,4} into exactly three blocks and one way to put them in an undirected cycle of length three. There is one partition of {1,2,3,4} into four blocks and 3 ways to make an undirected cycle of length four. 6 + 3 = 9. - _Geoffrey Critzer_, Nov 23 2012
		

References

  • 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).
  • Steffensen, J. F. Interpolation. 2d ed. Chelsea Publishing Co., New York, N. Y., 1950. ix+248 pp. MR0036799 (12,164d)

Crossrefs

A diagonal of the triangular array in A241168.
(1,2)-avoiding patterns are counted by A011782.
(1,1)-matching patterns are counted by A019472.
(1,2)-matching permutations are counted by A033312.
(1,2)-matching compositions are counted by A056823.
(1,2)-matching permutations of prime indices are counted by A335447.
(1,2)-matching compositions are ranked by A335485.
Patterns are counted by A000670 and ranked by A333217.
Patterns matched by compositions are counted by A335456.

Programs

  • Mathematica
    a[n_] := Sum[ k!*StirlingS2[n-1, k], {k, 0, n-1}] - 2^(n-2); Table[a[n], {n, 3, 17}] (* Jean-François Alcover, Nov 18 2011, after Manfred Goebel *)
    allnorm[n_]:=If[n<=0,{{}},Function[s,Array[Count[s,y_/;y<=#]+1&,n]]/@Subsets[Range[n-1]+1]];
    Table[Length[Select[Join@@Permutations/@allnorm[n],!GreaterEqual@@#&]],{n,0,5}] (* Gus Wiseman, Jun 24 2020 *)
  • PARI
    a(n) = sum(s=2, n-1, stirling(n,s+1,2)*s!/2); \\ Michel Marcus, Jun 24 2020

Formula

[n,2] = Sum_{s=2..n-1} Stirling2(n,s+1)*s!/2 (cf. A241168).
a(1)=0; for n >= 2, a(n) = A000670(n-1) - 2^(n-2). - Manfred Goebel (mkgoebel(AT)essex.ac.uk), Feb 20 2000; formula adjusted by N. J. A. Sloane, Apr 22 2014. For example, a(5) = 67 = A000670(4)-2^3 = 75-8 = 67.
E.g.f.: (1 - exp(2*x) - 2*log(2 - exp(x)))/4 = B(A(x)) where A(x) = exp(x)-1 and B(x) = (log(1/(1-x))- x - x^2/2)/2. - Geoffrey Critzer, Nov 23 2012

Extensions

Entry revised by N. J. A. Sloane, Apr 22 2014

A241170 Steffensen's bracket function [n,n-3].

Original entry on oeis.org

6, 25, 67, 145, 275, 476, 770, 1182, 1740, 2475, 3421, 4615, 6097, 7910, 10100, 12716, 15810, 19437, 23655, 28525, 34111, 40480, 47702, 55850, 65000, 75231, 86625, 99267, 113245, 128650, 145576, 164120, 184382, 206465, 230475, 256521, 284715, 315172
Offset: 3

Views

Author

N. J. A. Sloane, Apr 22 2014

Keywords

Crossrefs

A diagonal of the triangular array in A241168.

Programs

  • Magma
    [(n-2)*(n-1)*(24+7*n+3*n^2)/24: n in [3..50]]; // Vincenzo Librandi, Dec 12 2014
    
  • Magma
    I:=[6,25,67,145,275]; [n le 5 select I[n] else 5*Self(n-1)-10*Self(n-2)+10*Self(n-3)-5*Self(n-4)+Self(n-5): n in [1..50]]; // Vincenzo Librandi, Dec 12 2014
    
  • Maple
    with(combinat);
    T:=proc(n,k) add(stirling2(n,s+1)*s!/k!,s=k..n-1); end;
    [seq(T(n,n-3),n=3..16)];
  • Mathematica
    Table[(n-2)*(n-1)*(24 + 7*n + 3*n^2)/24,{n,3,20}] (* Vaclav Kotesovec, Apr 23 2014 *)
    CoefficientList[Series[(6 - 5 x + 2 x^2) / (1 - x)^5, {x, 0, 50}], x] (* Vincenzo Librandi, Dec 12 2014 *)
  • PARI
    for(n=3,30, print1((n-2)*(n-1)*(24+7*n+3*n^2)/24, ", ")) \\ G. C. Greubel, Feb 07 2018

Formula

See A241168.
a(n) = (n-2)*(n-1)*(24 + 7*n + 3*n^2)/24. - Vaclav Kotesovec, Apr 23 2014
From Vincenzo Librandi, Dec 12 2014: (Start)
G.f.: x^3*(6-5*x+2*x^2)/(1-x)^5.
a(n) = 5*a(n-1)-10*a(n-2)+10*a(n-3)-5*a(n-4)+a(n-5) for n>4.(End)

A241169 Steffensen's bracket function [n,3].

Original entry on oeis.org

0, 0, 0, 1, 14, 145, 1450, 15421, 180894, 2359225, 34072850, 540848341, 9363767974, 175619185105, 3547113529050, 76761061273261, 1771884886830254, 43456922321543785, 1128511554354422050, 30933862439582514181, 892562598747547111734, 27041608332832948251265, 858281473267724898703850
Offset: 1

Views

Author

N. J. A. Sloane, Apr 22 2014

Keywords

References

  • J. F. Steffensen, On a class of polynomials and their application to actuarial problems, Skandinavisk Aktuarietidskrift, 11 (1928), 75-97.

Crossrefs

A diagonal of the triangular array in A241168.

Programs

  • Maple
    with(combinat);
    T:=proc(n,k) add(stirling2(n,s+1)*s!/k!,s=k..n-1); end;
    [seq(T(n,3),n=1..16)];
  • Mathematica
    Flatten[{0,0,0,Table[Sum[StirlingS2[n, s+1]*s!/3!, {s,3,n-1}],{n,4,20}]}] (* Vaclav Kotesovec, Apr 22 2014 *)

Formula

See A241168.
a(n) ~ (n-1)! / (6*(log(2))^n). - Vaclav Kotesovec, Apr 22 2014
Showing 1-5 of 5 results.