%I#21 2021年1月12日09:06:52
%S 1,1,2,2,4,5,6,8,13,12,19,23,29,35,45,52,68,80,98111141163198230,
%电话283320376443517585719799932108512541417166818612138,
%电话:2449280431663666408346625277596066765651849496351080312157
%N避免等距三项算术级数的N的分区数。
%H Fausto A.C.Cariboni,n表,n=0..300的A(n)(术语0..150来自Joerg Arndt和Alois P.Heinz)
%e a(8)=13此类分区为:
%e 01:[1 1 2 4]
%e 02:[1 1 3 3]
%e 03:[1 1 6]
%e 04:[1 2 2 3]
%e 05:[1 2 5]
%e 06:[1 3 4]
%e 07:[17]
%e 08:[2 2 4]
%e 09:[2 3 3]
%e 10:[2 6]
%e 11:[3 5]
%e 12:[4 4]
%e 13:[8]
%e注意,第四个分区具有算术级数1,2,3,但不在等距位置。
%p b:=程序(n,i,l)局部j;
%从2到iquo的j的p(nops(l)+1,2)do
%p如果l[1]-l[j]=l[j]-l[2*j-1],则返回0 fiod;
%p`if`(n=0,1,`if`)(i<1,0,b(n,i-1,l)+
%p`if`(i>n,0,b(n-i,i,[i,l[]]))
%p端:
%pa:=n->b(n,n,[]):
%p序列(a(n),n=0..40);
%tb[n_,i_,l_]:=b[n,i,l]=模[{j},对于[j=2,j<=商[l]+1,2],j++,如果[l[[1]]-l[j]]==l[[j]]-l[[2*j-1]],返回[0]]];如果[n==0,1,如果[i<1,0,b[n,i-1,l]+如果[i>n,0,b[n-i,i,Prepend[l,i]]]];
%ta[n]:=b[n,n,{}];
%t表[a[n],{n,0,40}](*_Jean-François Alcover_,2018年5月21日,翻译自Maple*)
%Y参考A238432(成分相同)。
%Y参见A238571(避免任何三项算术级数的分区)。
%Y参见A238424(算术级数中避免三个连续部分的分区)。
%Y参考A238687。
%K nonn公司
%0、3
%A _Joerg Arndt_和_Alois P.Heinz,2014年3月1日
|