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-10 of 16 results. Next

A191646 Triangle read by rows: T(n,k) = number of connected multigraphs with n >= 0 edges and 1 <= k <= n+1 vertices, with no loops allowed.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 1, 2, 2, 0, 1, 3, 5, 3, 0, 1, 4, 11, 11, 6, 0, 1, 6, 22, 34, 29, 11, 0, 1, 7, 37, 85, 110, 70, 23, 0, 1, 9, 61, 193, 348, 339, 185, 47, 0, 1, 11, 95, 396, 969, 1318, 1067, 479, 106, 0, 1, 13, 141, 771, 2445, 4457, 4940, 3294, 1279, 235
Offset: 0

Views

Author

Alberto Tacchella, Jul 04 2011

Keywords

Examples

			Triangle T(n,k) (with rows n >= 0 and columns k >= 1) begins as follows:
  1;
  0, 1;
  0, 1, 1;
  0, 1, 2,  2;
  0, 1, 3,  5,  3;
  0, 1, 4, 11, 11,  6;
  0, 1, 6, 22, 34, 29, 11;
  ...
		

Crossrefs

Row sums give A076864. Diagonal is A000055.
Cf. A034253, A054923, A192517, A253186 (column k=3), A290778 (column k=4).

Programs

  • PARI
    EulerT(v)={my(p=exp(x*Ser(dirmul(v,vector(#v,n,1/n))))-1); Vec(p/x,-#v)}
    InvEulerMT(u)={my(n=#u, p=log(1+x*Ser(u)), vars=variables(p)); Vec(serchop( sum(i=1, n, moebius(i)*substvec(p + O(x*x^(n\i)), vars, apply(v->v^i,vars))/i), 1))}
    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(v,x)={sum(i=2, #v, sum(j=1, i-1, my(g=gcd(v[i],v[j])); g*x^(v[i]*v[j]/g))) + sum(i=1, #v, my(t=v[i]); ((t-1)\2)*x^t + if(t%2,0,x^(t/2)))}
    G(n,m)={my(s=0); forpart(p=n, s+=permcount(p)*EulerT(Vec(edges(p,x) + O(x*x^m), -m))); s/n!}
    R(n)={Mat(apply(p->Col(p+O(y^n),-n), InvEulerMT(vector(n, k, 1 + y*Ser(G(k,n-1), y)))))}
    { my(A=R(10)); for(n=1, #A, for(k=1, n, print1(A[n,k], ", "));print) } \\ Andrew Howroyd, May 14 2018

Formula

T(n,k=3) = A253186(n) = A034253(n,k=2) for n >= 1. - Petros Hadjicostas, Oct 02 2019

A054923 Triangle read by rows: number of connected graphs with k >= 0 edges and n nodes (1<=n<=k+1).

Original entry on oeis.org

1, 0, 1, 0, 0, 1, 0, 0, 1, 2, 0, 0, 0, 2, 3, 0, 0, 0, 1, 5, 6, 0, 0, 0, 1, 5, 13, 11, 0, 0, 0, 0, 4, 19, 33, 23, 0, 0, 0, 0, 2, 22, 67, 89, 47, 0, 0, 0, 0, 1, 20, 107, 236, 240, 106, 0, 0, 0, 0, 1, 14, 132, 486, 797, 657, 235, 0, 0, 0, 0, 0, 9, 138, 814, 2075, 2678, 1806, 551, 0, 0, 0, 0, 0, 5, 126, 1169, 4495, 8548, 8833, 5026, 1301
Offset: 0

Views

Author

Keywords

Comments

The diagonal n = k+1 is A000055(n). - Jonathan Vos Post, Aug 10 2008

Examples

			Triangle begins:
  1;
  0, 1;
  0, 0, 1;
  0, 0, 1, 2;
  0, 0, 0, 2, 3;
  0, 0, 0, 1, 5   6;
  0, 0, 0, 1, 5, 13,  11;
  0, 0, 0, 0, 4, 19,  33,  23;
  0, 0, 0, 0, 2, 22,  67,  89,  47;
  0, 0, 0, 0, 1, 20, 107, 236, 240, 106;
  ... (so with 5 edges there's 1 graph with 4 nodes, 5 with 5 nodes and 6 with 6 nodes). [Typo corrected by Anders Haglund, Jul 08 2008]
		

References

  • F. Harary and E. M. Palmer, Graphical Enumeration, Academic Press, NY, 1973, p. 93, Table 4.2.2; p. 241, Table A2.

Crossrefs

Main diagonal is A000055.
Subsequent diagonals give the number of connected unlabeled graphs with n nodes and n+k edges for k=0..2: A001429, A001435, A001436.
Cf. A002905 (row sums), A001349 (column sums), A008406, A046751 (transpose), A054924 (transpose), A046742 (w/o left column), A343088 (labeled).

Programs

  • PARI
    InvEulerMT(u)={my(n=#u, p=log(1+x*Ser(u)), vars=variables(p)); Vec(serchop( sum(i=1, n, moebius(i)*substvec(p + O(x*x^(n\i)), vars, apply(v->v^i,vars))/i), 1))}
    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(v,t) = {prod(i=2, #v, prod(j=1, i-1, my(g=gcd(v[i],v[j])); t(v[i]*v[j]/g)^g )) * prod(i=1, #v, my(c=v[i]); t(c)^((c-1)\2)*if(c%2, 1, t(c/2)))}
    G(n, x)={my(s=0); forpart(p=n, s+=permcount(p)*edges(p,i->1+x^i)); s/n!}
    T(n)={Mat([Col(p+O(y^n), -n) | p<-InvEulerMT(vector(n, k, G(k, y + O(y^n))))])}
    {my(A=T(10)); for(n=1, #A, print(A[n,1..n]))} \\ Andrew Howroyd, Oct 23 2019

Extensions

a(83)-a(89) corrected by Andrew Howroyd, Oct 24 2019

A076864 Number of connected loopless multigraphs with n edges.

Original entry on oeis.org

1, 1, 2, 5, 12, 33, 103, 333, 1183, 4442, 17576, 72810, 314595, 1410139, 6541959, 31322474, 154468852, 783240943, 4077445511, 21765312779, 118999764062, 665739100725, 3807640240209, 22246105114743, 132672322938379, 807126762251748
Offset: 0

Views

Author

N. J. A. Sloane, Nov 23 2002

Keywords

Comments

Inverse Euler transform of A050535.

Crossrefs

Programs

  • Mathematica
    A050535 = Cases[Import["https://oeis.org/A050535/b050535.txt", "Table"], {, }][[All, 2]];
    (* EulerInvTransform is defined in A022562 *)
    Join[{1}, EulerInvTransform[A050535 // Rest]] (* Jean-François Alcover, Feb 11 2020, updated Mar 17 2020 *)

Extensions

More terms from Sean A. Irvine, Oct 02 2011
Name and comment swapped by Gus Wiseman, Nov 28 2018
a(0)=1 prepended by Andrew Howroyd, Oct 23 2019

A191970 Number of connected graphs with n edges with loops allowed.

Original entry on oeis.org

1, 2, 2, 6, 12, 33, 93, 287, 940, 3309, 12183, 47133, 190061, 796405, 3456405, 15501183, 71681170, 341209173, 1669411182, 8384579797, 43180474608, 227797465130, 1229915324579, 6790642656907, 38311482445514, 220712337683628, 1297542216770482, 7779452884747298
Offset: 0

Views

Author

Alberto Tacchella, Jun 20 2011

Keywords

Comments

Inverse Euler transform of A053419.
From R. J. Mathar, Jul 25 2017: (Start)
The Multiset Transform gives the number of graphs with n edges (loops allowed) and k components (0<=k<=n):
1
0 2
0 2 3
0 6 4 4
0 12 15 6 5
0 33 36 24 8 6
0 93 111 64 33 10 7
0 287 324 207 92 42 12 8
0 940 1036 633 308 120 51 14 9
0 3309 3408 2084 966 409 148 60 16 10
0 12183 11897 6959 3243 1305 510 176 69 18 11
0 47133 43137 24415 10970 4432 1644 611 204 78 20 12
0 190061 163608 88402 38763 15125 5628 1983 712 232 87 22 13
0 796405 644905 332979 140671 53732 19316 6824 2322 813 260 96 24 14
0 3456405 2639871 1299054 529179 195517 68878 23515 8020 2661 914 288 105 26 15 (End)

Examples

			a(1)=2: Either one node with the edge equal to a loop, or two nodes connected by the edge. a(2)=2: Either three nodes on a chain connected by the two edges, or two nodes connected by an edge, one node with a loop. Apparently multi-loops are not allowed (?). - _R. J. Mathar_, Jul 25 2017
		

Crossrefs

Programs

  • PARI
    \\ See A322114 for InvEulerMT, G.
    seq(n)={vecsum([Vec(p+O(y^n), -n) | p<-InvEulerMT(vector(n, k, G(k, y + O(y^n))))])} \\ Andrew Howroyd, Oct 22 2019

Extensions

Terms a(25) and beyond from Andrew Howroyd, Oct 22 2019

A322115 Triangle read by rows where T(n,k) is the number of unlabeled connected multigraphs with loops with n edges and k vertices.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 4, 4, 2, 1, 6, 11, 9, 3, 1, 9, 25, 34, 20, 6, 1, 12, 52, 104, 99, 49, 11, 1, 16, 94, 274, 387, 298, 118, 23, 1, 20, 162, 645, 1295, 1428, 881, 300, 47, 1, 25, 263, 1399, 3809, 5803, 5088, 2643, 765, 106, 1, 30, 407, 2823, 10187, 20645, 24606, 17872, 7878, 1998, 235
Offset: 0

Views

Author

Gus Wiseman, Nov 26 2018

Keywords

Examples

			Triangle begins:
  1
  1   1
  1   2   1
  1   4   4   2
  1   6  11   9   3
  1   9  25  34  20   6
  1  12  52 104  99  49  11
		

Crossrefs

Row sums are A007719. Diagonal k = n-1 is A000055.

Programs

  • PARI
    EulerT(v)={my(p=exp(x*Ser(dirmul(v,vector(#v,n,1/n))))-1); Vec(p/x,-#v)}
    InvEulerMT(u)={my(n=#u, p=log(1+x*Ser(u)), vars=variables(p)); Vec(serchop( sum(i=1, n, moebius(i)*substvec(p + O(x*x^(n\i)), vars, apply(v->v^i,vars))/i), 1))}
    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(v,x)={sum(i=2, #v, sum(j=1, i-1, my(g=gcd(v[i],v[j])); g*x^(v[i]*v[j]/g))) + sum(i=1, #v, my(t=v[i]); ((t+1)\2)*x^t + if(t%2, 0, x^(t/2)))}
    G(n,m)={my(s=0); forpart(p=n, s+=permcount(p)*EulerT(Vec(edges(p,x) + O(x*x^m), -m))); s/n!}
    R(n)={Mat(apply(p->Col(p+O(y^n), -n), InvEulerMT(vector(n, k, 1 + y*Ser(G(k,n-1), y)))))}
    { my(T=R(10)); for(n=1, #T, print(T[n, 1..n])) } \\ Andrew Howroyd, Nov 30 2018

Extensions

Terms a(28) and beyond from Andrew Howroyd, Nov 30 2018

A070166 Irregular triangle read by rows giving T(n,k) = number of rooted graphs on n + 1 nodes with k edges (n >= 0, 0 <= k <= n(n-1)/2).

Original entry on oeis.org

1, 1, 1, 1, 2, 2, 1, 1, 2, 4, 6, 4, 2, 1, 1, 2, 5, 11, 17, 18, 17, 11, 5, 2, 1, 1, 2, 5, 13, 29, 52, 76, 94, 94, 76, 52, 29, 13, 5, 2, 1, 1, 2, 5, 14, 35, 83, 173, 308, 487, 666, 774, 774, 666, 487, 308, 173, 83, 35, 14, 5, 2, 1, 1, 2, 5, 14, 37, 98, 252, 585, 1239, 2396, 4135, 6340
Offset: 0

Views

Author

Keywords

Comments

T(n,k) is also the number of graphs with n nodes and k edges which may contain loops. (Delete the root node and change every edge leading to it into a loop.)
T(n,k) is also the number of symmetric relations with n points and k relations.

Examples

			Triangle begins:
1;
1, 1;
1, 2, 2, 1;
1, 2, 4, 6, 4, 2, 1;
1, 2, 5, 11, 17, 18, 17, 11, 5, 2, 1; <- gives either the numbers of rooted graphs on 5 nodes, or the numbers of graphs with loops on 4 nodes; with from 0 to 10 edges
1, 2, 5, 13, 29, 52, 76, 94, 94, 76, 52, 29, 13, 5, 2, 1;
...
		

References

  • E. Palmer and F. Harary, Graphical Enumeration, Academic Press, 1973.

Crossrefs

Programs

  • Mathematica
    Join[{{1},{1,1}},CoefficientList[Table[CycleIndex[Join[PairGroup[SymmetricGroup[n]],Permutations[Range[Binomial[n,2]+1,Binomial[n,2]+n]],2],s]/.Table[s[i]->1+x^i,{i,1,n^2-n}],{n,2,7}],x]]//Grid  (* Geoffrey Critzer, Oct 01 2012 *)
    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[v_, t_] := Product[Product[g = GCD[v[[i]], v[[j]]]; t[v[[i]]*v[[j]]/g]^g, {j, 1, i - 1}], {i, 2, Length[v]}]*Product[c = v[[i]]; t[c]^Quotient[c + 1, 2]*If[OddQ[c], 1, t[c/2]], {i, 1, Length[v]}];
    row[n_] := Module[{s = 0}, Do[s += permcount[p]*edges[p, 1 + x^# &], {p, IntegerPartitions[n]}]; s/n!] // Expand // CoefficientList[#, x] &
    row /@ Range[0, 7] // Flatten (* Jean-François Alcover, Jan 07 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(v, t) = {prod(i=2, #v, prod(j=1, i-1, my(g=gcd(v[i], v[j])); t(v[i]*v[j]/g)^g )) * prod(i=1, #v, my(c=v[i]); t(c)^((c+1)\2)*if(c%2, 1, t(c/2)))}
    G(n, A=0) = {my(s=0); forpart(p=n, s+=permcount(p)*edges(p, i->1+x^i+A)); s/n!}
    { for(n=0, 7, print(Vecrev(G(n)))) } \\ Andrew Howroyd, Oct 23 2019, updated Jan 09 2024

Extensions

Offset changed by Andrew Howroyd, Oct 23 2019

A322137 Number of labeled connected graphs with n edges (the vertices are {1,2,...,k} for some k).

Original entry on oeis.org

1, 1, 3, 17, 140, 1524, 20673, 336259, 6382302, 138525780, 3384988809, 91976158434, 2751122721402, 89833276321440, 3179852538140115, 121287919647418118, 4959343701136929850, 216406753768138678671, 10037782414506891597734, 493175891246093032826160
Offset: 0

Views

Author

Gus Wiseman, Nov 27 2018

Keywords

Crossrefs

Programs

  • Mathematica
    csm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[OrderedQ[#],UnsameQ@@#,Length[Intersection@@s[[#]]]>0]&]},If[c=={},s,csm[Union[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
    Table[Length[Select[Subsets[Subsets[Range[n+1],{2}],{n}],And[Union@@#==Range[Max@@Union@@#],Length[csm[#]]==1]&]],{n,6}]
  • PARI
    Connected(v)={my(u=vector(#v));for(n=1, #u, u[n]=v[n] - sum(k=1, n-1, binomial(n-1,k)*v[k]*u[n-k])); u}
    seq(n)={Vec(vecsum(Connected(vector(2*n, j, (1 + x + O(x*x^n))^binomial(j,2)))))} \\ Andrew Howroyd, Nov 28 2018

Extensions

Terms a(8) and beyond from Andrew Howroyd, Nov 28 2018

A368983 Number of connected graphs with loops (symmetric relations) on n unlabeled vertices with n edges.

Original entry on oeis.org

1, 1, 1, 3, 6, 14, 33, 81, 204, 526, 1376, 3648, 9792, 26485, 72233, 198192, 546846, 1515687, 4218564, 11782427, 33013541, 92759384, 261290682, 737688946, 2086993034, 5915398230, 16795618221, 47763406249, 136028420723, 387928330677, 1107692471888, 3166613486137
Offset: 0

Views

Author

Andrew Howroyd, Jan 11 2024

Keywords

Comments

The graphs considered here can have loops but not parallel edges.

Examples

			Representatives of the a(3) = 3 graphs are:
   {{1,2}, {1,3}, {2,3}},
   {{1}, {1,2}, {1,3}},
   {{1}, {1,2}, {2,3}}.
		

Crossrefs

A diagonal of A322114.
The labeled version is A368951.
Cf. A000081, A001429, A027852, A068051, A283755, A368984 (not necessarily connected).

Programs

  • 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 + (sum(d=1, n, eulerphi(d)/d*log(1/(1-g(d)))) + ((1+g(1))^2/(1-g(2))-1)/2 - (g(1)^2 + g(2)))/2)}

Formula

a(n) = A000081(n) + A001429(n) = A068051(n) - A027852(n) for n > 0.

A322147 Regular triangle read by rows where T(n,k) is the number of labeled connected graphs with loops with n edges and k vertices, 1 <= k <= n+1.

Original entry on oeis.org

1, 1, 1, 0, 2, 3, 0, 1, 10, 16, 0, 0, 12, 79, 125, 0, 0, 6, 162, 847, 1296, 0, 0, 1, 179, 2565, 11436, 16807, 0, 0, 0, 116, 4615, 47100, 185944, 262144, 0, 0, 0, 45, 5540, 121185, 987567, 3533720, 4782969, 0, 0, 0, 10, 4720, 220075, 3376450, 23315936, 76826061, 100000000
Offset: 0

Views

Author

Gus Wiseman, Nov 28 2018

Keywords

Examples

			Triangle begins:
  1
  1     1
  0     2     3
  0     1    10    16
  0     0    12    79   125
  0     0     6   162   847  1296
  0     0     1   179  2565 11436 16807
		

Crossrefs

Row sums are A322151. Last column is A000272.
Column sums are A062740.

Programs

  • Mathematica
    multsubs[set_,k_]:=If[k==0,{{}},Join@@Table[Prepend[#,set[[i]]]&/@multsubs[Drop[set,i-1],k-1],{i,Length[set]}]];
    csm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[OrderedQ[#],UnsameQ@@#,Length[Intersection@@s[[#]]]>0]&]},If[c=={},s,csm[Union[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
    Table[If[n==0,1,Length[Select[Subsets[multsubs[Range[k],2],{n}],And[Union@@#==Range[k],Length[csm[#]]==1]&]]],{n,0,6},{k,1,n+1}]
  • PARI
    Connected(v)={my(u=vector(#v)); for(n=1, #u, u[n]=v[n] - sum(k=1, n-1, binomial(n-1, k)*v[k]*u[n-k])); u}
    M(n)={Mat([Col(p, -(n+1)) | p<-Connected(vector(2*n, j, (1 + x + O(x*x^n) )^binomial(j+1,2)))[1..n+1]])}
    { my(T=M(10)); for(n=1, #T, print(T[n,][1..n])) } \\ Andrew Howroyd, Nov 29 2018

Extensions

Terms a(28) and beyond from Andrew Howroyd, Nov 29 2018

A322151 Number of labeled connected graphs with loops with n edges (the vertices are {1,2,...,k} for some k).

Original entry on oeis.org

1, 2, 5, 27, 216, 2311, 30988, 499919, 9431026, 203743252, 4960335470, 134382267082, 4009794148101, 130668970606412, 4617468180528235, 175867725701333896, 7182126650899080024, 313063334893103361130, 14507460736615554141354, 712192629608088061633746
Offset: 0

Views

Author

Gus Wiseman, Nov 28 2018

Keywords

Crossrefs

Row sums of A322147. The unlabeled version is A191970.

Programs

  • Mathematica
    multsubs[set_,k_]:=If[k==0,{{}},Join@@Table[Prepend[#,set[[i]]]&/@multsubs[Drop[set,i-1],k-1],{i,Length[set]}]];
    csm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[OrderedQ[#],UnsameQ@@#,Length[Intersection@@s[[#]]]>0]&]},If[c=={},s,csm[Union[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
    Table[Length[Select[Subsets[multsubs[Range[n+1],2],{n}],And[Union@@#==Range[Max@@Union@@#],Length[csm[#]]==1]&]],{n,5}]
  • PARI
    Connected(v)={my(u=vector(#v)); for(n=1, #u, u[n]=v[n] - sum(k=1, n-1, binomial(n-1, k)*v[k]*u[n-k])); u}
    seq(n)={Vec(vecsum(Connected(vector(2*n, j, (1 + x + O(x*x^n))^binomial(j+1,2)))))} \\ Andrew Howroyd, Nov 28 2018

Extensions

Terms a(7) and beyond from Andrew Howroyd, Nov 28 2018
Showing 1-10 of 16 results. Next