%I M3592 #29 Nov 07 2021 15:28:17
%S 0,1,1,4,22,178,2278,46380,1578060,92765486,9676866173,1821391854302,
%T 625710416245358,395761853562201960,464128290507379386872,
%U 1015085639712281997464676,4160440039279630394986003604,32088534920274236421098827156776
%N Number of unlabeled rooted nonseparable graphs with n nodes.
%D R. W. Robinson, personal communication.
%D R. W. Robinson, Numerical implementation of graph counting algorithms, AGRC Grant, Math. Dept., Univ. Newcastle, Australia, 1976.
%D N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
%H Andrew Howroyd, <a href="/A004115/b004115.txt">Table of n, a(n) for n = 1..40</a> (terms 1..26 from R. W. Robinson)
%H R. W. Robinson, <a href="/A004115/a004115.pdf">Rooted non-separable graphs - computer printout</a>
%o (PARI) \\ See links in A339645 for combinatorial species functions.
%o edges(v) = {sum(i=2, #v, sum(j=1, i-1, gcd(v[i], v[j]))) + sum(i=1, #v, v[i]\2)}
%o graphsCycleIndex(n)={my(s=0); forpart(p=n, s+=permcount(p) * 2^edges(p) * sMonomial(p)); s/n!}
%o graphsSeries(n)={sum(k=0, n, graphsCycleIndex(k)*x^k) + O(x*x^n)}
%o cycleIndexSeries(n)={my(g=graphsSeries(n), gcr=sPoint(g)/g); x*sSolve( sLog( gcr/(x*sv(1)) ), gcr )}
%o { my(N=15); Vec(OgfSeries(cycleIndexSeries(N)), -N) } \\ _Andrew Howroyd_, Dec 25 2020
%Y Cf. A002218, A007145, A013922.
%K nonn,nice
%O 1,4
%A _N. J. A. Sloane_
|