|
（PARI）\\SetTypes函数被其他序列引用。
重量（v）={Vec（exp（x*Ser（dirmul（v，vector（#v，n，（-1）^（n-1）/n）））-1，-v}
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（n，w）={sumdiv（gcd（n，w），d，moebius（d）*二项式（n/d，w/d））/n}
S（n）={my（v=向量（n））；对于（w=0，n，fordiv（gcd（n，w），d，v[n/d]+=x^w*v（n/d，w/d））；v}
集合类型（ptyp，fx）={
my（lim=和（i=1，#ptyp，ptyp[i]），u=向量（lim，i，O（x*x^（lim\i）））；u[1]+=1；
对于（i=1，#ptyp，my（s=s（ptyp[i]），v=向量（#u））；对于（j=1，#s，my（g=lcm（j，k））；如果（g<=35;v，v[g]+=u[j]*s[k]*j*k/g））；u=v）；
u[1]=1；Vec（和（i=1，#u，subst（fx（u[i]），x，x^i））+O（x*x^lim），-lim）；}
a（n）={my（s=0）；forpart（p=n，s+=permcount（p）*权重（SetTypes（p，q->q））[n]）；序列号！} \\安德鲁·豪罗伊德2019年9月1日
| 