(PARI)
lista(nn)={local(Cache=Map());
my(excl=向量(nn,n,sumdiv(n,d,2^(n-d)));
my(a(n,m=n,b=0)=
如果(n==0,1,
而(m>n||位测试(b,0),m--;b> >=1);
我的(hk=[n,m,b],z);
如果(!mapisdefined(缓存、hk和z),
z=如果(m,self()(n,m-1,b>>1)+self()(n-m,m,bitor(b,excl[m])),0);
mapput(缓存,hk,z));z) );
对于(n=1,nn,打印1(a(n),“,”)
|