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.

A057500 Number of connected labeled graphs with n edges and n nodes.

Original entry on oeis.org

0, 0, 1, 15, 222, 3660, 68295, 1436568, 33779340, 880107840, 25201854045, 787368574080, 26667815195274, 973672928417280, 38132879409281475, 1594927540549217280, 70964911709203684440, 3347306760024413356032, 166855112441313024389625, 8765006377126199463936000
Offset: 1

Views

Author

Qing-Hu Hou and David C. Torney (dct(AT)lanl.gov), Sep 01 2000

Keywords

Comments

Equivalently, number of connected unicyclic (i.e., containing one cycle) graphs on n labeled nodes. - Vladeta Jovovic, Oct 26 2004
a(n) is the number of trees on vertex set [n] = {1,2,...,n} rooted at 1 with one marked inversion (an inversion is a pair (i,j) with i > j and j a descendant of i in the tree). Here is a bijection from the title graphs (on [n]) to these marked trees. A title graph has exactly one cycle. There is a unique path from vertex 1 to this cycle, first meeting it at k, say (k may equal 1). Let i and j be the two neighbors of k in the cycle, with i the larger of the two. Delete the edge k<->j thereby forming a tree (in which j is a descendant of i) and take (i,j) as the marked inversion. To reverse this map, create a cycle by joining the smaller element of the marked inversion to the parent of the larger element. a(n) = binomial(n-1,2)*A129137(n). This is because, on the above marked trees, the marked inversion is uniformly distributed over 2-element subsets of {2,3,...,n} and so a(n)/binomial(n-1,2) is the number of trees on [n] (rooted at 1) for which (3,2) is an inversion. - David Callan, Mar 30 2007

Examples

			E.g., a(4)=15 because there are three different (labeled) 4-cycles and 12 different labeled graphs with a 3-cycle and an attached, external vertex.
		

References

  • F. Harary and E. M. Palmer, Graphical Enumeration, Academic Press, NY, 1973.
  • C. L. Mallows, Letter to N. J. A. Sloane, 1980.
  • R. J. Riddell, Contributions to the theory of condensation, Dissertation, Univ. of Michigan, Ann Arbor, 1951.

Crossrefs

A diagonal of A343088.
Cf. A000272 = labeled trees on n nodes; connected labeled graphs with n nodes and n+k edges for k=0..8: this sequence, A061540, A061541, A061542, A061543, A096117, A061544, A096150, A096224.
Cf. A001429 (unlabeled case), A052121.
For any number of edges we have A001187, unlabeled A001349.
This is the connected and covering case of A116508.
For #edges <= #nodes we have A129271, covering A367869.
For #edges > #nodes we have A140638, covering A367868.
This is the connected case of A367862 and A367863, unlabeled A006649.
The version with loops is A368951, unlabeled A368983.
This is the covering case of A370317.
Counting only covering vertices gives A370318.
A006125 counts graphs, A000088 unlabeled.
A006129 counts covering graphs, A002494 unlabeled.

Programs

  • Maple
    egf:= -1/2*ln(1+LambertW(-x)) +1/2*LambertW(-x) -1/4*LambertW(-x)^2:
    a:= n-> n!*coeff(series(egf, x, n+3), x, n):
    seq(a(n), n=1..25);  # Alois P. Heinz, Mar 27 2013
  • Mathematica
    nn=20; t=Sum[n^(n-1) x^n/n!, {n,1,nn}]; Drop[Range[0,nn]! CoefficientList[Series[Log[1/(1-t)]/2-t^2/4-t/2, {x,0,nn}], x], 1]  (* Geoffrey Critzer, Oct 07 2012 *)
    a[n_] := (n-1)!*n^n/2*Sum[1/(n^k*(n-k)!), {k, 3, n}]; Table[a[n], {n, 1, 20}] (* Jean-François Alcover, Jan 15 2014, after Vladeta Jovovic *)
    csm[s_]:=With[{c=Select[Subsets[Range[Length[s]],{2}],Length[Intersection@@s[[#]]]>0&]},If[c=={},s,csm[Sort[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
    Table[Length[Select[Subsets[Subsets[Range[n],{2}]],Union@@#==Range[n]&&Length[#]==n&&Length[csm[#]]<=1&]],{n,0,5}] (* Gus Wiseman, Feb 19 2024 *)
  • Sage
    # Warning: Floating point calculation. Adjust precision as needed!
    from mpmath import mp, chop, gammainc
    mp.dps = 200; mp.pretty = True
    for n in (1..100):
        print(chop((n^(n-2)*(1-3*n)+exp(n)*gammainc(n+1, n)/n)/2))
    # Peter Luschny, Jan 27 2016

Formula

The number of labeled connected graphs with n nodes and m edges is Sum_{k=1..n} (-1)^(k+1)/k*Sum_{n_1+n_2+..n_k=n, n_i>0} n!/(Product_{i=1..k} (n_i)!)* binomial(s, m), s=Sum_{i..k} binomial(n_i, 2). - Vladeta Jovovic, Apr 10 2001
E.g.f.: (1/2) Sum_{k>=3} T(x)^k/k, with T(x) = Sum_{n>=1} n^(n-1)/n! x^n. R. J. Riddell's thesis contains a closed-form expression for the number of connected graphs with m nodes and n edges. The present series applies to the special case m=n.
E.g.f.: -1/2*log(1+LambertW(-x))+1/2*LambertW(-x)-1/4*LambertW(-x)^2. - Vladeta Jovovic, Jul 09 2001
Asymptotic expansion (with xi=sqrt(2*Pi)): n^(n-1/2)*[xi/4-7/6*n^(-1/2)+xi/48* n^(-1)+131/270*n^(-3/2)+xi/1152*n^(-2)+4/2835*n^(-5/2)+O(n^(-3))]. - Keith Briggs, Aug 16 2004
Row sums of A098909: a(n) = (n-1)!*n^n/2*Sum_{k=3..n} 1/(n^k*(n-k)!). - Vladeta Jovovic, Oct 26 2004
a(n) = Sum_{k=0..C(n-1,2)} k*A052121(n,k). - Alois P. Heinz, Nov 29 2015
a(n) = (n^(n-2)*(1-3*n)+exp(n)*Gamma(n+1,n)/n)/2. - Peter Luschny, Jan 27 2016
a(n) = A062734(n,n+1) = A123527(n,n). - Gus Wiseman, Feb 19 2024

Extensions

More terms from Vladeta Jovovic, Jul 09 2001

A062734 Triangular array T(n,k) giving number of connected graphs with n labeled nodes and k edges (n >= 1, 0 <= k <= n(n-1)/2).

Original entry on oeis.org

1, 0, 1, 0, 0, 3, 1, 0, 0, 0, 16, 15, 6, 1, 0, 0, 0, 0, 125, 222, 205, 120, 45, 10, 1, 0, 0, 0, 0, 0, 1296, 3660, 5700, 6165, 4945, 2997, 1365, 455, 105, 15, 1, 0, 0, 0, 0, 0, 0, 16807, 68295, 156555, 258125, 331506, 343140, 290745, 202755, 116175, 54257, 20349
Offset: 1

Views

Author

Vladeta Jovovic, Jul 12 2001

Keywords

Comments

T(n,n-1) = n^(n-2) counts free labeled trees A000272.
T(n,n) counts labeled connected unicyclic graphs A057500. - Geoffrey Critzer, Oct 07 2012

Examples

			Triangle starts:
[1],
[0, 1],
[0, 0, 3,  1],
[0, 0, 0, 16,  15,   6,   1],
[0, 0, 0,  0, 125, 222, 205, 120, 45, 10, 1],
...
		

References

  • Cowan, D. D.; Mullin, R. C.; Stanton, R. G. Counting algorithms for connected labelled graphs. Proceedings of the Sixth Southeastern Conference on Combinatorics, Graph Theory, and Computing (Florida Atlantic Univ., Boca Raton, Fla., 1975), pp. 225-236. Congressus Numerantium, No. XIV, Utilitas Math., Winnipeg, Man., 1975. MR0414417 (54 #2519). - N. J. A. Sloane, Apr 06 2012
  • F. Harary and E. Palmer, Graphical Enumeration, Academic Press, 1973, Page 29, Exercise 1.5.

Crossrefs

Cf. A001187 (row sums), A054924 (unlabeled case), A061540 (a subdiagonal).
See A123527 for another version (without leading zeros in each row).

Programs

  • Mathematica
    nn=6;s=Sum[(1+y)^Binomial[n,2] x^n/n!,{n,0,nn}]; Range[0,nn]!CoefficientList[Series[Log[ s]+1,{x,0,nn}],{x,y}]//Grid  (* returns triangle indexed at n = 0, Geoffrey Critzer, Oct 07 2012 *)
    T[ n_, k_] := If[ n < 0, 0, Coefficient[ n! SeriesCoefficient[ Log[ Sum[ (1 + y)^Binomial[m, 2] x^m/m!, {m, 0, n}]], {x, 0, n}], y, k]]; (* Michael Somos, Aug 12 2017 *)
  • PARI
    {T(n, k) = if( n<0, 0, n! * polcoeff( polcoeff( log( sum(m=0, n, (1 + y)^(m * (m-1)/2) * x^m/m!)), n), k))}; /* Michael Somos, Aug 12 2017 */

Formula

G.f.: Sum_{n>=1, k>=0} T(n, k) * x^n/n! * y^k = log(Sum_{n>=0} (1 + y)^binomial(n, 2) * x^n/n!). - Ralf Stephan, Jan 18 2005

A123534 Triangular array T(n,k) giving number of 2-connected graphs with n labeled nodes and k edges (n >= 3, n <= k <= n(n-1)/2).

Original entry on oeis.org

1, 3, 6, 1, 12, 70, 100, 45, 10, 1, 60, 720, 2445, 3535, 2697, 1335, 455, 105, 15, 1, 360, 7560, 46830, 133581, 216951, 232820, 183540, 111765, 53627, 20307, 5985, 1330, 210, 21, 1, 2520, 84000, 835800, 3940440, 10908688, 20317528
Offset: 3

Views

Author

N. J. A. Sloane, Nov 13 2006

Keywords

Examples

			Triangle begins (n >= 3, k >= n):
  n
  3 | 1;
  4 | 3, 6, 1;
  5 | 12, 70, 100, 45, 10, 1;
  6 | 60, 720, 2445, 3535, 2697, 1335, 455, 105, 15, 1;
  ...
		

References

  • R. W. Robinson, Numerical implementation of graph counting algorithms, AGRC Grant, Math. Dept., Univ. Newcastle, Australia, 1977.

Crossrefs

Row sums give A013922.

Programs

  • Mathematica
    row[n_] := row[n] = Module[{s}, s = (n-1)!*Log[x/InverseSeries[#, x]& @ (x*D[#, x]& @ Log[Sum[(1+y)^Binomial[k, 2]*x^k/k!, {k, 0, n}] + O[x]^(n+1) ])]; CoefficientList[Coefficient[s, x, n-1]/y^n, y]];
    Table[row[n], {n, 3, 15}] // Flatten (* Jean-François Alcover, Aug 13 2019, after Andrew Howroyd *)
  • PARI
    row(n)={Vecrev((n-1)!*polcoef(log(x/serreverse(x*deriv(log(sum(k=0, n, (1 + y)^binomial(k, 2) * x^k / k!) + O(x*x^n))))), n-1)/y^n)}
    { for(n=3, 7, print(row(n))) } \\ Andrew Howroyd, Nov 30 2018

A123542 Triangular array T(n,k) giving number of 3-connected graphs with n labeled nodes and k edges (n >= 4, ceiling(3*n/2) <= k <= n(n-1)/2).

Original entry on oeis.org

1, 15, 10, 1, 70, 492, 690, 395, 105, 15, 1, 5040, 28595, 58905, 63990, 42392, 18732, 5880, 1330, 210, 21, 1, 16800, 442680, 2485920, 6629056, 10684723, 11716068, 9409806, 5824980, 2872317, 1147576, 373156, 98112, 20475, 3276
Offset: 4

Views

Author

N. J. A. Sloane, Nov 13 2006

Keywords

Examples

			Triangle begins:
n = 4
k = 6 : 1
Total( 4) = 1
n = 5
k = 8 : 15
k = 9 : 10
k = 10 : 1
Total( 5) = 26
n = 6
k = 9 : 70
k = 10 : 492
k = 11 : 690
k = 12 : 395
k = 13 : 105
k = 14 : 15
k = 15 : 1
Total( 6) = 1768
n = 7
k = 11 : 5040
k = 12 : 28595
k = 13 : 58905
k = 14 : 63990
k = 15 : 42392
k = 16 : 18732
k = 17 : 5880
k = 18 : 1330
k = 19 : 210
k = 20 : 21
k = 21 : 1
Total( 7) = 225096
		

References

  • R. W. Robinson, Numerical implementation of graph counting algorithms, AGRC Grant, Math. Dept., Univ. Newcastle, Australia, 1977.

Crossrefs

Row sums give A005644. Cf. A123527, A123534.

A123525 Arises in the normal ordering of functions of a*(a+)*a, where a and a+ are the boson annihilation and creation operators, respectively.

Original entry on oeis.org

2, 14, 102, 836, 7730, 79962, 916454, 11533832, 158149026, 2346622310, 37458934502, 640013453004, 11652216012242, 225169809833906, 4602407562991590, 99194703240441872
Offset: 1

Views

Author

Karol A. Penson, Oct 02 2006

Keywords

Crossrefs

Programs

  • Mathematica
    Rest[With[{nmax = 50}, CoefficientList[Series[(1/(1 - x)^2)*Exp[x/(1 - x)]*LaguerreL[1, 1/(x - 1)]*x, {x, 0, nmax}], x]*Range[0, nmax]!]] (* G. C. Greubel, Oct 14 2017 *)

Formula

E.g.f.: (1/(1-x)^2)*exp(x/(1-x))*LaguerreL(1,1/(x-1))*x.
From Vaclav Kotesovec, Nov 13 2017: (Start)
Recurrence: (n-2)*(n-1)*a(n) = 2*(n-2)*n^2*a(n-1) - (n-1)^3*n*a(n-2).
a(n) ~ exp(2*sqrt(n) - n - 1/2) * n^(n + 5/4) / sqrt(2) * (1 + 31/(48*sqrt(n))).
(End)
Showing 1-5 of 5 results.