%I#41 2019年8月14日15:28:55
%S 1,1,1,1,1,1,1,1,2,1,1,1,2,1,1,2,1,1,1,3,2,2,1,2,1,3,2,1,0,0,1,0,
%T 1,1,4,2,4,1,5,1,1,0,1,0,0,1,0,1,4,3,4,1,7,1,1,2,0,2,0,1,0,0,0,0,
%U 0,1,1,5,3,6,2,9,1,2,1,3,0,4,0,1,0,0,0,1,1,0,00,0,1,15,36,2
%N不规则表T(N,k):N的分区数,其中所有部分的最小公共倍数等于k。
%H Alois P.Heinz,行n=0..30,扁平</a>
%e n=4的5个分区是1+1+1+1(lcm=1)、1+1+2(lcm=2)、2+2(lcm=2),1+3(lcm=3)和4(lcm=4)。所以k=1、3和4出现一次,k=2出现两次。
%e三角形开始:
%e 1;
%e 1;
%e 11;
%e 11 1;
%e 1 2 1 1;
%e 1 2 1 1 1 1;
%e 1 3 2 2 1 2;
%e 1 3 2 1 3 1 0 0 1 0 1;
%e。。。
%p A256067:=程序(n,k)
%p局部a,p;
%p a:=0;
%组合[划分](n)do中p的p
%p ilcm(op(p));
%p如果%=k,则
%pa:=a+1;
%p end if;
%p端do:
%p a;
%p端程序:
%p#第二个Maple程序:
%p b:=proc(n,i)选项记忆`如果`(n=0或i=1,x,
%p b(n,i-1)+(p->加(系数(p,x,t)*x^ilcm(t,i),
%p t=1.度(p))(添加(b(n-i*j,i-1),j=1..n/i))
%p端:
%pT:=n->(p->seq(系数(p,x,i),i=1..度(p)))(b(n$2)):
%p序列(T(n),n=0..12);#_阿洛伊斯·海因茨(Alois P.Heinz),2015年3月27日
%tb[n_,i_]:=b[n,i]=如果[n==0||i==1,x,b[n、i-1]+函数[{p},总和[系数[p,x,t]*x^LCM[t,i],{t,1,指数[p,x]}][总和[b[n-i*j,i-1],{j,1,n/i}]];T[n_]:=函数[{p},表[系数[p,x,i],{i,1,指数[p,x]}][b[n,n]];表[T[n],{n,0,12}]//扁平(*_Jean-François Alcover_,2015年6月22日,以_Alois P.Heinz_*命名)
%Y参见A000041(行总和)、A000793(行长度)、A213952、A074761(对角线)、A074752(第6列)、A008642(第4列)、C002266(第5列)、P002264(第3列)、A132270(第7列)、AO08643(第8列)、A308649(第9列)、A258470(第10列)。
%Y参考A009490(行的非零项数)、A074064(行的最后一个元素)、A168532(gcd相同)、A181844(总和k*T(n,k))。
%K nonn,标签
%O 0.9
%A R.J.Mathar_,2015年3月18日
%E T(0,1)=1由_Alois P.Heinz_编写,2015年3月27日
|