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] ；
边[p_]：=总和[天花板[（p[i]]-1）*（p[[i]-2）/6）]，{i，1，长度[p]}]+总和[c=p[i]；d=p[[j]]；GCD[c，d]*（c+d-2+Mod[（c-d）/GCD[c，d]，2]）/2+总和[c*d*p[[k]]/LCM[c，d，p[k]]，{k，1，j-1}]，{j，1，i-1}]，}i，2，长度[p]}]；
a[n_]：=模块[{s=0}，Do[s+=permcount[p]*2^edges[p]，{p，IntegerPartitions[n]}]；s/n！]；
a/@范围[0，12](*Jean-François Alcover公司2021年1月8日之后安德鲁·霍罗伊德*)
