permcount(v)={my(m=1,s=0,k=0,t);对于(i=1,#v,t=v[i];k=if(i>1&&t==v[i-1],k+1,1);m*=t*k;s+=t);s!/m}
边(v)={和(i=2,#v,和(j=1,i-1,2*gcd(v[i],v[j]))+和(i=1,#v,v[i]-1)}
C(p)={((i,v)->if(i>#p,2^边(v),my(s=0);forpart(q=p[i],s+=置换数(q)*self()(i+1,concat(v,Vec(q)));s/p[i]!))(1,[])}
多重(v)={my(p=(#v)!,k=1);对于(i=2,#v,k=if(v[i]==v[i-1],k+1,p/=k!;1);p/k!}
wC(p)=乘积(p)*C(p
行(n)={[vecsum(应用(wC,vecsort([Vecrev(p)|p<-分区(n),#p==m],4))|m<-[1..n]]}
{对于(n=0,10,打印(行(n)))}