<<数学世界`Graphs`
长度/@(gp=选择[#,GraphicalPartitionQ]&/@
图形/@范围[9])
nn=20;g=总和[NumberOfGraphs[n]x^n,{n,0,nn}];系数列表[系列[g(1-x),{x,0,nn}],x](*Geoffrey Critzer,2012年4月14日*)
sysnorm[m_]:=如果[Union@@m!=范围[Max@@Flatten[m]],sysnorm[m/.Rule@@@表[{(Union@@m)[[i]],i},{i,长度[Union@m]}]],第一个[Sort[sysnormal[m,1]]];
sysnorm[m_,aft_]:=If[Length[Union@@m]<=aft,{m},With[{mx=Table[Count[m,i,{2}],{i,Select[Union@@m,#>=aft&]}]},Union@@(sysnorm[#,aft+1]&/@Union[Table[Map[Sort,m/.{par+aft-1->aft,aft->par+aft_1},{0,1}],},[par,First/@Position[mx,Max[mx]]}])]])];
表[Length[Union[sysnorm/@Select[Subsets[Select[Subsets[Range[n]],Length[#]==2&]],Union@@#=Range[n]&]]],{n,6}](*古斯·怀斯曼2018年8月2日*)
b[n_,i_,l_]:=如果[n==0||i==1,1/n!*2^(函数[p,Sum[Ceiling[(p[j]-1)/2]+Sum[GCD[p[[k]],p[[j]],{k,1,j-1}],{j,1,长度[p]}]][Join[l,表[1,{n}]]),总和[b[n-i*j,i-1,连接[l,表格[i,{j}]]]]/j/i^j,{j,0,n/i}]];
a[n]:=b[n,n,{}]-如果[n>0,b[n-1,n-1,{}],0];