%I#18 2017年11月19日08:02:11
%S 1,1,2,2,4,4,6,8,11,12,18,22,28,34,44,54,69,82102125154185226,
%电话:2713273934745626737979471124132915631846216425412974,
%电话:3480406247385508640437432861499661153013307153451767020337
%N按频率小于或等于其在总和列表中位置的部分划分。
%C设分区的和为s(1)<s(2)<…<s(k)和s(i)的频率为f(i)。然后我们计算f(i)<=i的分区。
%H Alois P.Heinz,n的表格,n=0..2000的a(n)</a>
%e分区1+1不计算在内,因为它的最小部分1出现了两次。
%e计算分区3+2+2+1是因为它的最小部分1出现一次;它的第二个最小部分2出现两次(2<=2),第三个部分3出现一次(1<=3)。
%pb:=proc(n,i,t)选项记忆`如果`(n=0,1,`if`(n<i,0,
%p加(b(n-i*j,i+1,t+`如果`(j=0,0,1)),j=0..分钟(t,n/i))
%p端:
%pa:=n->b(n,1美元2):
%p序列(a(n),n=0..60);#_Alois P.Heinz,2017年11月18日
%t<<组合数学`;
%t nend=30;
%t对于[n=1,n<=nend,n++,计数[n]=0;
%t部分=分区[n];
%t对于[i=1,i<=长度[part],i++,
%t t=理货[部分[i]]];
%t条件=真;
%t对于[j=1,j<=长度[t],j++,
%t如果[t[[-j,2]]>j,条件=假]];
%t如果[条件,计数[n]+]]];
%t打印[表[count[i],{i,1,nend}]]
%Y参考A244395。
%K nonn公司
%0、4
%A _David S.Newman,2017年11月18日
%E来自_Alois P.Heinz的更多条款,2017年11月18日
|