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.

A122082 Number of unlabeled bicolored graphs on 2n nodes which are invariant when the two color classes are interchanged.

Original entry on oeis.org

1, 2, 5, 16, 67, 404, 3904, 64840, 1930842, 104698904, 10401039400, 1900637187280, 641429385018832, 401454435464761376, 467919402404052870944, 1019758699013228238271040, 4171161230867751509749228304
Offset: 0

Views

Author

N. J. A. Sloane, Oct 18 2006

Keywords

References

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

Crossrefs

Row sums of A123548.

Programs

  • Mathematica
    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_] := Sum[GCD[v[[i]], v[[j]]], {i, 2, Length[v]}, {j, 1, i - 1}] + Total @ Quotient[v + 1, 2]
    a[n_] := (s=0; Do[s += permcount[p]*2^edges[p], {p, IntegerPartitions[n]}]; s/n!);
    Table[a[n], {n, 0, 20}] (* Jean-François Alcover, Jul 06 2018, 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) = {sum(i=2, #v, sum(j=1, i-1, gcd(v[i],v[j]))) + sum(i=1, #v, (v[i]+1)\2)}
    a(n) = {my(s=0); forpart(p=n, s+=permcount(p)*2^edges(p)); s/n!} \\ Andrew Howroyd, Oct 23 2017

Formula

a(n) = 2*A007139(n) - A002724(n). - Vladeta Jovovic, Feb 27 2007

Extensions

More terms from Vladeta Jovovic, Feb 27 2007