连接[{1，2}，表[CycleIndex[Join[PairGroup[SymmetricGroup[n]，Ordered]，排列[Range[n^2-n+1，n^2]，2]，s]/。表[s[i]->2，{i，1，n^2-n}]，{n，2，7}]](*杰弗里·克雷策2011年11月2日*）
permcount[v_]：=模[{m=1，s=0，k=0，t}，对于[i=1，i<=长度[v]，i++，t=v[i]]；k=如果[i>1&&t==v[[i-1]]，k+1，1]；m*=t*k；s+=t]；s/m] ；
edges[v_]：=Sum[2*GCD[v[[i]]，v[[j]]]，{i，2，Length[v]}，{j，1，i-1}]+总计[v]；
a[n_]：=（s=0；Do[s+=permcount[p]*2^edges[p]，{p，IntegerPartitions[n]}]；s/n！）；
dinorm[m_]：=If[m=={}，{}；If[Union@@m！=Range[Max@@Flatten[m]]，dinorm[m/.应用[Rule，Table[{（Union@@m）[[i]]，i}，｛i，Length[Union@m]}]，{1}]]，First[Sort[dinorm[m，1]]]]；
dinorm[m_，aft_]：=If[Length[Union@@m]<=aft，{m}，With[{mx=Table[Count[m，i，{2}]，{i，Select[Union@@m，#1>=aft&]}]}，Union@@（dinorm[#1，aft+1]&）/@Union[Table[Map[Sort，m/.{par+aft-1->aft，aft->par+aft-1}，{0}]，}，第一个/@位置[mx，Max[mx]}]]]]；
表[Length[Union[diform/@Subsets[Tuples[Range[n]，2]]]，{n，0，3}](*古斯·怀斯曼2019年6月17日*）
