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.

Previous Showing 21-30 of 153 results. Next

A369199 Irregular triangle read by rows where T(n,k) is the number of labeled loop-graphs covering n vertices with k edges.

Original entry on oeis.org

1, 0, 1, 0, 1, 3, 1, 0, 0, 6, 17, 15, 6, 1, 0, 0, 3, 46, 150, 228, 206, 120, 45, 10, 1, 0, 0, 0, 45, 465, 1803, 3965, 5835, 6210, 4955, 2998, 1365, 455, 105, 15, 1, 0, 0, 0, 15, 645, 5991, 27364, 79470, 165555, 264050, 334713, 344526, 291200, 202860, 116190, 54258, 20349, 5985, 1330, 210, 21, 1
Offset: 0

Views

Author

Gus Wiseman, Jan 18 2024

Keywords

Examples

			Triangle begins:
   1
   0   1
   0   1   3   1
   0   0   6  17  15   6   1
   0   0   3  46 150 228 206 120  45  10   1
Row n = 3 counts the following loop-graphs (loops shown as singletons):
  {1,23}   {1,2,3}     {1,2,3,12}    {1,2,3,12,13}   {1,2,3,12,13,23}
  {2,13}   {1,2,13}    {1,2,3,13}    {1,2,3,12,23}
  {3,12}   {1,2,23}    {1,2,3,23}    {1,2,3,13,23}
  {12,13}  {1,3,12}    {1,2,12,13}   {1,2,12,13,23}
  {12,23}  {1,3,23}    {1,2,12,23}   {1,3,12,13,23}
  {13,23}  {1,12,13}   {1,2,13,23}   {2,3,12,13,23}
           {1,12,23}   {1,3,12,13}
           {1,13,23}   {1,3,12,23}
           {2,3,12}    {1,3,13,23}
           {2,3,13}    {1,12,13,23}
           {2,12,13}   {2,3,12,13}
           {2,12,23}   {2,3,12,23}
           {2,13,23}   {2,3,13,23}
           {3,12,13}   {2,12,13,23}
           {3,12,23}   {3,12,13,23}
           {3,13,23}
           {12,13,23}
		

Crossrefs

The version without loops is A054548.
This is the covering case of A084546.
Column sums are A173219.
Row sums are A322661, unlabeled A322700.
The connected case is A369195, without loops A062734.
A000085, A100861, A111924 count set partitions into singletons or pairs.
A006125 counts simple graphs; also loop-graphs if shifted left.
A006129 counts covering graphs, unlabeled A002494.
A368927 counts choosable loop-graphs, covering A369140.
A369141 counts non-choosable loop-graphs, covering A369142.

Programs

  • Mathematica
    Table[Length[Select[Subsets[Subsets[Range[n],{1,2}],{k}],Length[Union@@#]==n&]],{n,0,5},{k,0,Binomial[n+1,2]}]
  • PARI
    T(n)={[Vecrev(p) | p<-Vec(serlaplace(exp(-x + O(x*x^n))*(sum(j=0, n, (1 + y)^binomial(j+1, 2)*x^j/j!)))) ]}
    { my(A=T(6)); for(i=1, #A, print(A[i])) } \\ Andrew Howroyd, Feb 02 2024

Formula

E.g.f.: exp(-x) * (Sum_{j >= 0} (1 + y)^binomial(j+1, 2)*x^j/j!). - Andrew Howroyd, Feb 02 2024

A005703 Number of n-node connected graphs with at most one cycle.

Original entry on oeis.org

1, 1, 1, 2, 4, 8, 19, 44, 112, 287, 763, 2041, 5577, 15300, 42419, 118122, 330785, 929469, 2621272, 7411706, 21010378, 59682057, 169859257, 484234165, 1382567947, 3952860475, 11315775161, 32430737380, 93044797486, 267211342954, 768096496093, 2209772802169
Offset: 0

Views

Author

Keywords

Comments

a(n) is the number of pseudotrees on n nodes. - Eric W. Weisstein, Jun 11 2012
Also unlabeled connected graphs covering n vertices with at most n edges. For this definition we have a(1) = 0 and possibly a(0) = 0. - Gus Wiseman, Feb 20 2024

Examples

			From _Gus Wiseman_, Feb 20 2024: (Start)
Representatives of the a(0) = 1 through a(5) = 8 graphs:
  {}  .  {12}  {12,13}     {12,13,14}     {12,13,14,15}
               {12,13,23}  {12,13,24}     {12,13,14,25}
                           {12,13,14,23}  {12,13,24,35}
                           {12,13,24,34}  {12,13,14,15,23}
                                          {12,13,14,23,25}
                                          {12,13,14,23,45}
                                          {12,13,14,25,35}
                                          {12,13,24,35,45}
(End)
		

References

  • J. Riordan, An Introduction to Combinatorial Analysis, Wiley, 1958, p. 150.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A000055, A000081, A001429 (labeled A057500), A134964 (number of pseudoforests, labeled A133686).
The labeled version is A129271.
The connected complement is A140636, labeled A140638.
Non-connected: A368834 (labeled A367869) or A370316 (labeled A369191).
A001187 counts connected graphs, unlabeled A001349.
A006125 counts simple graphs, unlabeled A000088.
A006129 counts covering graphs, unlabeled A002494.
A062734 counts connected graphs by number of edges.

Programs

  • Mathematica
    Needs["Combinatorica`"]; nn = 20; t[x_] := Sum[a[n] x^n, {n, 1, nn}];
    a[0] = 0;
    b = Drop[Flatten[
        sol = SolveAlways[
          0 == Series[
            t[x] - x Product[1/(1 - x^i)^a[i], {i, 1, nn}], {x, 0, nn}],
          x]; Table[a[n], {n, 0, nn}] /. sol], 1];
    r[x_] := Sum[b[[n]] x^n, {n, 1, nn}]; c =
    Drop[Table[
        CoefficientList[
         Series[CycleIndex[DihedralGroup[n], s] /.
           Table[s[i] -> r[x^i], {i, 1, n}], {x, 0, nn}], x], {n, 3,
         nn}] // Total, 1];
    d[x_] := Sum[c[[n]] x^n, {n, 1, nn}]; CoefficientList[
    Series[r[x] - (r[x]^2 - r[x^2])/2 + d[x] + 1, {x, 0, nn}], x] (* Geoffrey Critzer, Nov 17 2014 *)
  • PARI
    \\ TreeGf gives gf of A000081.
    TreeGf(N)={my(A=vector(N, j, 1)); for (n=1, N-1, A[n+1] = 1/n * sum(k=1, n, sumdiv(k, d, d*A[d]) * A[n-k+1] ) ); x*Ser(A)}
    seq(n)={my(t=TreeGf(n)); my(g(e)=subst(t + O(x*x^(n\e)), x, x^e) + O(x*x^n)); Vec(1 + g(1) + (g(2) - g(1)^2)/2 + sum(k=3, n, sumdiv(k, d, eulerphi(d)*g(d)^(k/d))/k + if(k%2, g(1)*g(2)^(k\2), (g(1)^2+g(2))*g(2)^(k/2-1)/2))/2)}; \\ Andrew Howroyd and Washington Bomfim, May 15 2021

Formula

a(n) = A000055(n) + A001429(n).

Extensions

More terms from Vladeta Jovovic, Apr 19 2000 and from Michael Somos, Apr 26 2000
a(27) corrected and a(28) and a(29) computed by Washington Bomfim, May 14 2008

A140638 Number of connected graphs on n labeled nodes that contain at least two cycles.

Original entry on oeis.org

0, 0, 0, 7, 381, 21748, 1781154, 249849880, 66257728763, 34495508486976, 35641629989151608, 73354595357480683904, 301272202621204113362497, 2471648811029413368450098688, 40527680937730440155535277704046, 1328578958335783199341353852258282496
Offset: 1

Views

Author

Washington Bomfim, May 21 2008

Keywords

Comments

These are the connected graphs that are neither trees nor unicyclic.
Also connected non-choosable graphs covering n vertices, where a graph is choosable iff it is possible to choose a different vertex from each edge. The unlabeled version is A140636. The complement is counted by A129271. - Gus Wiseman, Feb 20 2024

References

  • J. Riordan, An Introduction to Combinatorial Analysis, Dover, 2002, p. 2.

Crossrefs

The unlabeled version is A140636.
Cf. A000272 (trees), A001187 (connected graphs), A057500 (connected unicyclic graphs).
The complement is counted by A129271, unlabeled A005703.
The non-connected complement is A133686, covering A367869.
The non-connected version is A367867, unlabeled A140637.
The non-connected covering version is A367868.
A006125 counts graphs, A000088 unlabeled.
A006129 counts covering graphs, A002494 unlabeled.
A143543 counts simple labeled graphs by number of connected components.

Programs

  • Mathematica
    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[csm[#]]<=1&&Select[Tuples[#],UnsameQ@@#&]=={}&]],{n,0,5}] (* Gus Wiseman, Feb 19 2024 *)
  • PARI
    seq(n)={my(A=O(x*x^n), t=-lambertw(-x + A)); Vec(serlaplace( log(sum(k=0, n, 2^binomial(k, 2)*x^k/k!, A)) - log(1/(1-t))/2 - t/2 + 3*t^2/4), -n)} \\ Andrew Howroyd, Jan 15 2022

Formula

a(n) = A001187(n) - A129271(n).
a(n) = A001187(n) - A000272(n) - A057500(n).

Extensions

Definition clarified by Andrew Howroyd, Jan 15 2022

A137917 a(n) is the number of unlabeled graphs on n nodes whose components are unicyclic graphs.

Original entry on oeis.org

1, 0, 0, 1, 2, 5, 14, 35, 97, 264, 733, 2034, 5728, 16101, 45595, 129327, 368093, 1049520, 2999415, 8584857, 24612114, 70652441, 203075740, 584339171, 1683151508, 4852736072, 14003298194, 40441136815, 116880901512, 338040071375, 978314772989, 2833067885748, 8208952443400
Offset: 0

Views

Author

Washington Bomfim, Feb 24 2008

Keywords

Comments

a(n) is the number of simple unlabeled graphs on n nodes whose components have exactly one cycle. - Geoffrey Critzer, Oct 12 2012
Also the number of unlabeled simple graphs with n vertices and n edges such that it is possible to choose a different vertex from each edge. - Gus Wiseman, Jan 25 2024

Examples

			From _Gus Wiseman_, Jan 25 2024: (Start)
Representatives of the a(0) = 1 through a(5) = 5 simple graphs:
  {}  .  .  {12,13,23}  {12,13,14,23}  {12,13,14,15,23}
                        {12,13,24,34}  {12,13,14,23,25}
                                       {12,13,14,23,45}
                                       {12,13,14,25,35}
                                       {12,13,24,35,45}
(End)
		

Crossrefs

The connected case is A001429.
Without the choice condition we have A001434, covering A006649.
For any number of edges we have A134964, complement A140637.
The labeled version is A137916.
The version with loops is A369145, complement A368835.
The complement is counted by A369201, labeled A369143, covering A369144.
A006129 counts covering graphs, unlabeled A002494.
A054548 counts graphs covering n vertices with k edges, with loops A369199.
A129271 counts connected choosable simple graphs, unlabeled A005703.

Programs

  • Mathematica
    Needs["Combinatorica`"];
    nn=30;s[n_,k_]:=s[n,k]=a[n+1-k]+If[n<2k,0,s[n-k,k]];a[1]=1;a[n_]:=a[n]=Sum[a[i]s[n-1,i]i,{i,1,n-1}]/(n-1);rt=Table[a[i],{i,1,nn}];c=Drop[Apply[Plus,Table[Take[CoefficientList[CycleIndex[DihedralGroup[n],s]/.Table[s[j]->Table[Sum[rt[[i]]x^(k*i),{i,1,nn}],{k,1,nn}][[j]],{j,1,nn}],x],nn],{n,3,nn}]],1];CoefficientList[Series[Product[1/(1-x^i)^c[[i]],{i,1,nn-1}],{x,0,nn}],x]   (* Geoffrey Critzer, Oct 12 2012, after code given by Robert A. Russell in A000081 *)
    brute[m_]:=First[Sort[Table[Sort[Sort/@(m/.Rule@@@Table[{(Union@@m)[[i]],p[[i]]},{i,Length[p]}])],{p,Permutations[Range[Length[Union@@m]]]}]]];
    Table[Length[Union[brute/@Select[Subsets[Subsets[Range[n],{2}],{n}],Select[Tuples[#],UnsameQ@@#&]!={}&]]],{n,0,5}] (* Gus Wiseman, Jan 25 2024 *)

Formula

a(n) = Sum_{1*j_1 + 2*j_2 + ... = n} (Product_{i=3..n} binomial(A001429(i) + j_i -1, j_i)). [F. Ruskey p. 79, (4.27) with n replaced by n+1, and a_i replaced by A001429(i)].
Euler transform of A001429. - Geoffrey Critzer, Oct 12 2012

Extensions

Edited by Washington Bomfim, Jun 27 2012
Terms a(30) and beyond from Andrew Howroyd, May 05 2018
Offset changed to 0 by Gus Wiseman, Jan 27 2024

A004251 Number of graphical partitions (degree-vectors for simple graphs with n vertices, or possible ordered row-sum vectors for a symmetric 0-1 matrix with diagonal values 0).

Original entry on oeis.org

1, 1, 2, 4, 11, 31, 102, 342, 1213, 4361, 16016, 59348, 222117, 836315, 3166852, 12042620, 45967479, 176005709, 675759564, 2600672458, 10029832754, 38753710486, 149990133774, 581393603996, 2256710139346, 8770547818956, 34125389919850, 132919443189544, 518232001761434, 2022337118015338, 7898574056034636, 30873421455729728
Offset: 0

Views

Author

Keywords

Comments

In other words, a(n) is the number of graphic sequences of length n, where a graphic sequence is a sequence of numbers which can be the degree sequence of some graph.
In the article by A. Iványi, G. Gombos, L. Lucz, and T. Matuszka, "Parallel enumeration of degree sequences of simple graphs II", in Table 4 on page 260 the values a(30) = 7898574056034638 and a(31) = 30873429530206738 are incorrect due to the incorrect Gz(30) = 5876236938019300 and Gz(31) = 22974847474172100. - Wang Kai, Jun 05 2016

Examples

			For n = 3, there are 4 different graphic sequences possible: 0 0 0; 1 1 0; 2 1 1; 2 2 2. - Daan van Berkel (daan.v.berkel.1980(AT)gmail.com), Jun 25 2010
From _Gus Wiseman_, Dec 31 2020: (Start)
The a(0) = 1 through a(4) = 11 sorted degree sequences:
  ()  (0)  (0,0)  (0,0,0)  (0,0,0,0)
           (1,1)  (0,1,1)  (0,0,1,1)
                  (1,1,2)  (0,1,1,2)
                  (2,2,2)  (0,2,2,2)
                           (1,1,1,1)
                           (1,1,1,3)
                           (1,1,2,2)
                           (1,2,2,3)
                           (2,2,2,2)
                           (2,2,3,3)
                           (3,3,3,3)
For example, the graph {{2,3},{2,4}} has degrees (0,2,1,1), so (0,1,1,2) is counted under a(4).
(End)
		

References

  • R. A. Brualdi and H. J. Ryser, Combinatorial Matrix Theory, Cambridge Univ. Press, 1992.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • P. R. Stein, On the number of graphical partitions, pp. 671-684 of Proc. 9th S-E Conf. Combinatorics, Graph Theory, Computing, Congr. Numer. 21 (1978).

Crossrefs

Counting the positive partitions by sum gives A000569, ranked by A320922.
The version with half-loops is A029889, with covering case A339843.
The covering case (no zeros) is A095268.
Covering simple graphs are ranked by A309356 and A320458.
Non-graphical partitions are counted by A339617 and ranked by A339618.
The version with loops is A339844, with covering case A339845.
A006125 counts simple graphs, with covering case A006129.
A027187 counts partitions of even length, ranked by A028260.
A058696 counts partitions of even numbers, ranked by A300061.
A320921 counts connected graphical partitions.
A322353 counts factorizations into distinct semiprimes.
A339659 counts graphical partitions of 2n into k parts.
A339661 counts factorizations into distinct squarefree semiprimes.

Programs

  • Mathematica
    Table[Length[Union[Sort[Table[Count[Join@@#,i],{i,n}]]&/@Subsets[Subsets[Range[n],{2}]]]],{n,0,5}] (* Gus Wiseman, Dec 31 2020 *)

Formula

G.f. = 1 + x + 2*x^2 + 4*x^3 + 11*x^4 + 31*x^5 + 102*x^6 + 342*x^7 + 1213*x^8 + ...
a(n) ~ c * 4^n / n^(3/4) for some constant c > 0. Computational estimates suggest c ≈ 0.099094. - Tom Johnston, Jan 18 2023

Extensions

More terms from Torsten Sillke, torsten.sillke(AT)lhsystems.com, using Cor. 6.3.3, Th. 6.3.6, Cor. 6.2.5 of Brualdi-Ryser.
a(19) from Herman Jamke (hermanjamke(AT)fastmail.fm), May 19 2007
a(20)-a(23) from Nathann Cohen, Jul 09 2011
a(24)-a(29) from Antal Iványi, Nov 15 2011
a(30) and a(31) corrected by Wang Kai, Jun 05 2016

A123023 a(n) = (n-1)*a(n-2), a(0)=1, a(1)=0.

Original entry on oeis.org

1, 0, 1, 0, 3, 0, 15, 0, 105, 0, 945, 0, 10395, 0, 135135, 0, 2027025, 0, 34459425, 0, 654729075, 0, 13749310575, 0, 316234143225, 0, 7905853580625, 0, 213458046676875, 0, 6190283353629375, 0, 191898783962510625, 0, 6332659870762850625, 0, 221643095476699771875
Offset: 0

Views

Author

Roger L. Bagula, Sep 24 2006

Keywords

Comments

a(n) is the number of ways of separating n terms into pairs. - Stephen Crowley, Apr 07 2007
a(n) is the n-th moment of the standard normal distribution. - Hal M. Switkay, Nov 06 2019
a(n) is the number of fixed-point free involutions in the symmetric group of degree n. - Nick Krempel, Feb 26 2020

Examples

			From _Gus Wiseman_, Dec 23 2018: (Start)
The a(6) = 15 ways of partitioning {1,2,3,4,5,6} into disjoint pairs:
  {{12}{34}{56}},  {{12}{35}{46}},  {{12}{36}{45}},
  {{13}{24}{56}},  {{13}{25}{46}},  {{13}{26}{45}},
  {{14}{23}{56}},  {{14}{25}{36}},  {{14}{26}{35}},
  {{15}{23}{46}},  {{15}{24}{36}},  {{15}{26}{34}},
  {{16}{23}{45}},  {{16}{24}{35}},  {{16}{25}{34}}.
(End)
		

References

  • Richard Bronson, Schaum's Outline of Modern Introductory Differential Equations, MacGraw-Hill, New York, 1973, page 107, solved problem 19.18
  • Norbert Wiener, Nonlinear Problems in Random Theory, 1958, Equation 1.31

Crossrefs

Programs

  • Magma
    a:=[1,0]; [n le 2 select a[n] else (n-2)*Self(n-2): n in [1..30]]; // Marius A. Burtea, Nov 07 2019
  • Maple
    with(combstruct): ZL2 := [S, {S=Set(Cycle(Z, card=2))}, labeled]:
    seq(count(ZL2, size=n), n=0..36); # Zerinvary Lajos, Sep 24 2007
    a := n -> ifelse(irem(n, 2) = 1, 0, 2^(n/2) * pochhammer(1/2, n/2)):
    seq(a(n), n = 0..36); # Peter Luschny, Jan 11 2023
  • Mathematica
    RecurrenceTable[{a[0] == 1, a[1] == 0, a[n] == (n - 1) a[n - 2]}, a[n], {n, 0, 31}] (* Ray Chandler, Jul 30 2015 *)

Formula

a(n) = (1/2)*Gamma((1/2)*n + 1/2)*2^((1/2)*n)*(1 + (-1)^n)/sqrt(Pi). - Stephen Crowley, Apr 07 2007
E.g.f.: exp(x^2/2). - Geoffrey Critzer, Mar 15 2009
a(2n) = A001147(n). - R. J. Mathar, Oct 11 2011
From Sergei N. Gladkovskii, Nov 18 2012, Dec 05 2012, May 16 2013, May 24 2013, Jun 07 2013: (Start)
Continued fractions:
E.g.f.: E(0) where E(k) = 1 + x^2*(4*k+1)/((4*k+2)*(4*k+3) - x^2*(4*k+2)*(4*k+3)^2/(x^2*(4*k+3) + (4*k+4)*(4*k+5)/E(k+1))).
G.f.: 1/G(0) where G(k) = 1 - x^2*(k+1)/G(k+1).
G.f.: 1 + x^2/(1+x) + Q(0)*x^3/(1+x), where Q(k) = 1 + (2*k+3)*x/(1 - x/(x + 1/Q(k+1))).
G.f.: G(0)/2, where G(k) = 1 + 1/(1-x/(x+1/x/(2*k+1)/G(k+1))).
G.f.: (G(0) - 1)*x/(1+x) + 1, where G(k) = 1 + x*(2*k+1)/(1 - x/(x + 1/G(k+1))). (End)
For n even, a(n) = A001147(n/2) = A124794(3^(n/2)). a(n) is also the coefficient of x1*...*xn in Product_{1 <= i < j <= n} (1 + xi*xj). - Gus Wiseman, Dec 23 2018
a(n) = 2^(n/2)*Pochhammer(1/2, n/2)*(n+1 mod 2). - Peter Luschny, Jan 11 2023

Extensions

Edited by N. J. A. Sloane, Jan 06 2008
Better name by Sergei N. Gladkovskii, May 24 2013
Leading term 1 dropped, offset changed, and entry edited correspondingly by Andrey Zabolotskiy, Nov 07 2019

A368597 Number of n-element sets of singletons or pairs of distinct elements of {1..n} with union {1..n}, or loop-graphs covering n vertices with n edges.

Original entry on oeis.org

1, 1, 3, 17, 150, 1803, 27364, 501015, 10736010, 263461265, 7283725704, 223967628066, 7581128184175, 280103206674480, 11216492736563655, 483875783716549277, 22371631078155742764, 1103548801569848115255, 57849356643299101021960, 3211439288584038922502820
Offset: 0

Views

Author

Gus Wiseman, Jan 04 2024

Keywords

Comments

It doesn't matter for this sequence whether we use loops such as {x,x} or half-loops such as {x}.

Examples

			The a(0) = 1 through a(3) = 17 set-systems:
  {}  {{1}}  {{1},{2}}    {{1},{2},{3}}
             {{1},{1,2}}  {{1},{2},{1,3}}
             {{2},{1,2}}  {{1},{2},{2,3}}
                          {{1},{3},{1,2}}
                          {{1},{3},{2,3}}
                          {{2},{3},{1,2}}
                          {{2},{3},{1,3}}
                          {{1},{1,2},{1,3}}
                          {{1},{1,2},{2,3}}
                          {{1},{1,3},{2,3}}
                          {{2},{1,2},{1,3}}
                          {{2},{1,2},{2,3}}
                          {{2},{1,3},{2,3}}
                          {{3},{1,2},{1,3}}
                          {{3},{1,2},{2,3}}
                          {{3},{1,3},{2,3}}
                          {{1,2},{1,3},{2,3}}
		

Crossrefs

This is the covering case of A014068.
Allowing edges of any positive size gives A054780, covering case of A136556.
Allowing any number of edges gives A322661, connected A062740.
The case of just pairs is A367863, covering case of A116508.
The unlabeled version is A368599.
The version contradicting strict AOC is A368730.
The connected case is A368951.
A000085 counts set partitions into singletons or pairs.
A006129 counts covering graphs, connected A001187.
A058891 counts set-systems, unlabeled A000612.
A100861 counts set partitions into singletons or pairs by number of pairs.
A111924 counts set partitions into singletons or pairs by length.

Programs

  • Mathematica
    Table[Length[Select[Subsets[Subsets[Range[n],{1,2}], {n}],Union@@#==Range[n]&]],{n,0,5}]
  • PARI
    a(n) = sum(k=0, n, (-1)^(n-k) * binomial(n,k) * binomial(binomial(k+1,2), n)) \\ Andrew Howroyd, Jan 06 2024

Formula

a(n) = Sum_{k=0..n} (-1)^(n-k) * binomial(n,k) * binomial(binomial(k+1,2), n). - Andrew Howroyd, Jan 06 2024

Extensions

Terms a(7) and beyond from Andrew Howroyd, Jan 06 2024

A368927 Number of labeled loop-graphs covering a subset of {1..n} such that it is possible to choose a different vertex from each edge.

Original entry on oeis.org

1, 2, 7, 39, 314, 3374, 45630, 744917, 14245978, 312182262, 7708544246, 211688132465, 6397720048692, 210975024924386, 7537162523676076, 289952739051570639, 11949100971787370300, 525142845422124145682, 24515591201199758681892, 1211486045654016217202663
Offset: 0

Views

Author

Gus Wiseman, Jan 15 2024

Keywords

Comments

These are loop-graphs where every connected component has a number of edges less than or equal to the number of vertices. Also loop-graphs with at most one cycle (unicyclic) in each connected component.

Examples

			The a(0) = 1 through a(2) = 7 loop-graphs (loops shown as singletons):
  {}  {}     {}
      {{1}}  {{1}}
             {{2}}
             {{1,2}}
             {{1},{2}}
             {{1},{1,2}}
             {{2},{1,2}}
		

Crossrefs

Without the choice condition we have A006125.
The case of a unique choice is A088957, unlabeled A087803.
The case without loops is A133686, complement A367867, covering A367869.
For exactly n edges and no loops we have A137916, unlabeled A137917.
For exactly n edges we have A333331 (maybe), complement A368596.
For edges of any positive size we have A367902, complement A367903.
The covering case is A369140, complement A369142.
The complement is counted by A369141.
The complement for n edges and no loops is A369143, covering A369144.
The unlabeled version is A369145, complement A369146.
A000085, A100861, A111924 count set partitions into singletons or pairs.
A006129 counts covering graphs, unlabeled A002494.
A322661 counts labeled covering loop-graphs, connected A062740.

Programs

  • Mathematica
    Table[Length[Select[Subsets[Subsets[Range[n],{1,2}]], Length[Select[Tuples[#],UnsameQ@@#&]]!=0&]],{n,0,5}]
  • PARI
    seq(n)={my(t=-lambertw(-x + O(x*x^n))); Vec(serlaplace(exp(3*t/2 - 3*t^2/4)/sqrt(1-t) ))} \\ Andrew Howroyd, Feb 02 2024

Formula

Binomial transform of A369140.
Exponential transform of A369197 with A369197(1) = 2.
E.g.f.: exp(3*T(x)/2 - 3*T(x)^2/4)/sqrt(1-T(x)), where T(x) is the e.g.f. of A000169. - Andrew Howroyd, Feb 02 2024

Extensions

a(7) onwards from Andrew Howroyd, Feb 02 2024

A369141 Number of labeled loop-graphs covering a subset of {1..n} such that it is not possible to choose a different vertex from each edge (non-choosable).

Original entry on oeis.org

0, 0, 1, 25, 710, 29394, 2051522, 267690539, 68705230758, 35184059906570, 36028789310419722, 73786976083150073999, 302231454897259573627852, 2475880078570549574773324062, 40564819207303333310731978895956, 1329227995784915872613854321228773937
Offset: 0

Views

Author

Gus Wiseman, Jan 20 2024

Keywords

Comments

Also labeled loop-graphs having at least one connected component containing more edges than vertices.

Examples

			The a(0) = 0 through a(3) = 25 loop-graphs (loops shown as singletons):
  .  .  {{1},{2},{1,2}}  {{1},{2},{1,2}}
                         {{1},{3},{1,3}}
                         {{2},{3},{2,3}}
                         {{1},{2},{3},{1,2}}
                         {{1},{2},{3},{1,3}}
                         {{1},{2},{3},{2,3}}
                         {{1},{2},{1,2},{1,3}}
                         {{1},{2},{1,2},{2,3}}
                         {{1},{2},{1,3},{2,3}}
                         {{1},{3},{1,2},{1,3}}
                         {{1},{3},{1,2},{2,3}}
                         {{1},{3},{1,3},{2,3}}
                         {{2},{3},{1,2},{1,3}}
                         {{2},{3},{1,2},{2,3}}
                         {{2},{3},{1,3},{2,3}}
                         {{1},{1,2},{1,3},{2,3}}
                         {{2},{1,2},{1,3},{2,3}}
                         {{3},{1,2},{1,3},{2,3}}
                         {{1},{2},{3},{1,2},{1,3}}
                         {{1},{2},{3},{1,2},{2,3}}
                         {{1},{2},{3},{1,3},{2,3}}
                         {{1},{2},{1,2},{1,3},{2,3}}
                         {{1},{3},{1,2},{1,3},{2,3}}
                         {{2},{3},{1,2},{1,3},{2,3}}
                         {{1},{2},{3},{1,2},{1,3},{2,3}}
		

Crossrefs

Without the choice condition we have A006125, unlabeled A000088.
The case of a unique choice is A088957, unlabeled A087803.
The case without loops is A367867, covering A367868.
For edges of any positive size we have A367903, complement A367902.
For exactly n edges we have A368596, complement A333331 (maybe).
The complement is counted by A368927, covering A369140.
The covering case is A369142.
For n edges and no loops we have A369143, covering A369144.
The unlabeled version is A369146 (covering A369147), complement A369145.
A000085, A100861, A111924 count set partitions into singletons or pairs.
A006129 counts covering graphs, unlabeled A002494.
A054548 counts graphs covering n vertices with k edges, with loops A369199.
A133686 counts choosable graphs, covering A367869.
A322661 counts labeled covering loop-graphs, unlabeled A322700.

Programs

  • Mathematica
    Table[Length[Select[Subsets[Subsets[Range[n], {1,2}]],Length[Select[Tuples[#],UnsameQ@@#&]]==0&]],{n,0,5}]

Formula

Binomial transform of A369142.
a(n) = A006125(n + 1) - A368927(n). - Andrew Howroyd, Feb 02 2024

Extensions

a(6) onwards from Andrew Howroyd, Feb 02 2024

A000665 Number of 3-uniform hypergraphs on n unlabeled nodes, or equivalently number of relations with 3 arguments on n nodes.

Original entry on oeis.org

1, 1, 1, 2, 5, 34, 2136, 7013320, 1788782616656, 53304527811667897248, 366299663432194332594005123072, 1171638318502989084030402509596875836036608, 3517726593606526072882013063011594224625680712384971214848
Offset: 0

Views

Author

Keywords

Comments

The Qian reference has one incorrect term. The formula given in corollary 2.6 also contains a minor error. The second summation needs to be over p_i*p_j*p_h/lcm(p_i, p_j, p_h) rather than gcd(p_i, p_j, p_h)^2. - Andrew Howroyd, Dec 11 2018

Examples

			From _Gus Wiseman_, Dec 13 2018: (Start)
Non-isomorphic representatives of the a(5) = 34 hypergraphs:
  {}
  {{123}}
  {{125}{345}}
  {{134}{234}}
  {{123}{245}{345}}
  {{124}{134}{234}}
  {{135}{245}{345}}
  {{145}{245}{345}}
  {{123}{124}{134}{234}}
  {{123}{145}{245}{345}}
  {{124}{135}{245}{345}}
  {{125}{135}{245}{345}}
  {{134}{235}{245}{345}}
  {{145}{235}{245}{345}}
  {{123}{124}{135}{245}{345}}
  {{123}{145}{235}{245}{345}}
  {{124}{134}{235}{245}{345}}
  {{134}{145}{235}{245}{345}}
  {{135}{145}{235}{245}{345}}
  {{145}{234}{235}{245}{345}}
  {{123}{124}{134}{235}{245}{345}}
  {{123}{134}{145}{235}{245}{345}}
  {{123}{145}{234}{235}{245}{345}}
  {{124}{135}{145}{235}{245}{345}}
  {{125}{135}{145}{235}{245}{345}}
  {{135}{145}{234}{235}{245}{345}}
  {{123}{124}{135}{145}{235}{245}{345}}
  {{124}{135}{145}{234}{235}{245}{345}}
  {{125}{135}{145}{234}{235}{245}{345}}
  {{134}{135}{145}{234}{235}{245}{345}}
  {{123}{124}{135}{145}{234}{235}{245}{345}}
  {{125}{134}{135}{145}{234}{235}{245}{345}}
  {{124}{125}{134}{135}{145}{234}{235}{245}{345}}
  {{123}{124}{125}{134}{135}{145}{234}{235}{245}{345}}
(End)
		

References

  • F. Harary and E. M. Palmer, Graphical Enumeration, Academic Press, NY, 1973, p. 231.
  • 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

Row sums of A092337. Spanning 3-uniform hypergraphs are counted by A322451.
Column k=3 of A309858.

Programs

  • Mathematica
    (* about 85 seconds on a laptop computer *)
    Needs["Combinatorica`"];Table[A = Subsets[Range[n],{3}];CycleIndex[Replace[Map[Sort,System`PermutationReplace[A, SymmetricGroup[n]], {2}],Table[A[[i]] -> i, {i, 1, Length[A]}], 2], s] /. Table[s[i] -> 2, {i, 1, Binomial[n, 3]}], {n, 1, 8}] (* Geoffrey Critzer, Oct 28 2015 *)
    Table[Sum[2^PermutationCycles[Ordering[Map[Sort,Subsets[Range[n],{3}]/.Rule@@@Table[{i,prm[[i]]},{i,n}],{1}]],Length],{prm,Permutations[Range[n]]}]/n!,{n,8}] (* Gus Wiseman, Dec 13 2018 *)
    permcount[v_] := Module[{m = 1, s = 0, k = 0, t}, For[i = 1, i <= Length[v], i++, t = v[[i]]; k = If[i > 1 && t == v[[i - 1]], k + 1, 1]; m *= t*k; s += t]; s!/m];
    edges[p_] := Sum[Ceiling[(p[[i]] - 1)*((p[[i]] - 2)/6)], {i, 1, Length[p]}] + Sum[Sum[c = p[[i]]; d = p[[j]]; GCD[c, d]*(c + d - 2 + Mod[(c - d)/GCD[c, d], 2])/2 + Sum[c*d*p[[k]]/LCM[c, d, p[[k]]], {k, 1, j - 1}], {j, 1, i - 1}], {i, 2, Length[p]}];
    a[n_] := Module[{s = 0}, Do[s += permcount[p]*2^edges[p], {p, IntegerPartitions[n]}]; s/n!];
    a /@ Range[0, 12] (* Jean-François Alcover, Jan 08 2021, after Andrew Howroyd *)
  • PARI
    permcount(v) = {my(m=1,s=0,k=0,t); for(i=1,#v,t=v[i]; k=if(i>1&&t==v[i-1],k+1,1); m*=t*k;s+=t); s!/m}
    edges(p)={sum(i=1, #p, ceil((p[i]-1)*(p[i]-2)/6)) + sum(i=2, #p, sum(j=1, i-1, my(c=p[i], d=p[j]); gcd(c,d)*(c + d - 2 + (c-d)/gcd(c,d)%2)/2 + sum(k=1, j-1, c*d*p[k]/lcm(lcm(c,d), p[k]))))}
    a(n) = {my(s=0); forpart(p=n, s+=permcount(p)*2^edges(p)); s/n!} \\ Andrew Howroyd, Dec 11 2018

Extensions

Corrected and extended by Vladeta Jovovic
a(0)=1 prepended and a(12) from Andrew Howroyd, Dec 11 2018
Previous Showing 21-30 of 153 results. Next