|
%I #23 Apr 29 2023 08:10:34
%S 1,2,44,3572,1005584,1060875152,4382913876704,71987098738435232,
%T 4721068803628864289024,1237845578934919489219757312,
%U 1298046978912816702510086132201984,5444486716626952189940499391640815580672,91343710775311761525117954724021374685703481344
%N Number of labeled connected digraphs with n nodes such that complement is also connected.
%H Alois P. Heinz, <a href="/A054914/b054914.txt">Table of n, a(n) for n = 1..50</a>
%H V. A. Liskovets, <a href="http://www.cs.uwaterloo.ca/journals/JIS/VOL3/LISK/Derseq.html">Some easily derivable sequences</a>, J. Integer Sequences, 3 (2000), #00.2.2.
%F a(n) = 2*A003027(n) - A053763(n).
%p b:= n-> 2^(n^2-n):
%p g:= proc(n) option remember; local k; `if`(n=0, 1,
%p b(n)- add(k*binomial(n,k) *b(n-k)*g(k), k=1..n-1)/n)
%p end:
%p a:= n-> 2*g(n)-b(n):
%p seq (a(n), n=1..20); # _Alois P. Heinz_, Oct 21 2012
%t nn=20; g=Sum[2^(2Binomial[n,2])x^n/n!,{n,0,nn}];
%t Drop[Range[0,nn]!CoefficientList[Series[2(Log[g]+1)-g,{x,0,nn}],x],1] (* _Geoffrey Critzer_, Oct 21 2012 *)
%o (Magma)
%o m:=30;
%o f:= func< x | (&+[2^(n*(n-1))*x^n/Factorial(n): n in [0..m+3]]) >;
%o R<x>:=PowerSeriesRing(Rationals(), m);
%o Coefficients(R!(Laplace( 1 + 2*Log(f(x)) - f(x) ))); // _G. C. Greubel_, Apr 28 2023
%o (SageMath)
%o m=30
%o def f(x): return sum(2^(n*(n-1))*x^n/factorial(n) for n in range(m+4))
%o def A054914_list(prec):
%o P.<x> = PowerSeriesRing(QQ, prec)
%o return P( 2 + 2*log(f(x)) - f(x) ).egf_to_ogf().list()
%o a=A054914_list(40); a[1:] # _G. C. Greubel_, Apr 28 2023
%Y Cf. A003027, A053763.
%K nonn,easy
%O 1,2
%A _N. J. A. Sloane_, May 23 2000
%E More terms from _Vladeta Jovovic_, Jul 17 2000
| 