|
（PARI）
集合（sig，v，r，x）={forstep（r=r，1，-1，my（w=sig[r]+1）；v=向量（#v，k，和（j=1，k，二项式（#v-j，k-j）*v[j]*x^（k-j）*（1+x）^（w-#v+j-1））-v[k]）；v[#v]}
DomSetCount（sig，x）={my（v=[1]）；my（total=Collect（sig、v、#sig、x））；forstep（r=#sig，1，-1，my（w=sig[r]+1）；total+=Collect（1+x）^（j-1））；总计}
a（n）=DomSetCount（Vecrev（[1..n]），1）\\安德鲁·霍罗伊德2017年11月5日
| 