%I#39 2023年9月5日21:23:57
%S 1,1,1,1,3,2,1,6,11,6,1,10,31,34,12,1,15,81189182,60,12,1141393,
%电话:494282,60,1828813803245374020344201,364562716829313268,
%电话:113384908840、1,45726557822207473515759840602155642520
%N行读取的三角形:T(N,k)=c(N,k)}lcm(c)中的Sum_{c,其中c(N、k)是{1,2,…,N}的所有k-子集的集合。
%C C(n,k)也称为n与大小k的组合(参见A181842)。
%C主对角线给出:A003418。下对角线给出:A094308。第k=1列给出:A000217.-_Alois P.Heinz,2013年7月29日
%H Alois P.Heinz,行n=0..46,扁平</a>
%电子[0]1
%e【1】1 1
%e[2]1 3 2
%e[3]1 6 11 6
%电子[4]1 10 31 34 12
%电子[5]1 15 81 189 182 60
%电子[6]1 21 141 393 494 282 60
%p与(combstruct):
%p a181853_row:=进程(n)局部k,L,L,R,comb;
%p R:=空;
%p代表k从0到n do
%p L:=0;
%p梳:=iterstructs(组合(n),大小=k):
%p未完成时(梳)do
%p l:=nextstruct(梳);
%p L:=L+ilcm(op(L));
%p od;
%p R:=R,L;
%p od;
%p R端:
%p#第二个Maple程序:
%p b:=proc(n,k)选项记忆`如果`(k=0,[1],
%p[`如果`(k<n,b(n-1,k),[])[],seq(ilcm(c,n),c=b(n-1,k-1))])
%p端:
%p T:=(n,k)->添加(c,c=b(n,k)):
%p-seq(seq(T(n,k),k=0..n),n=0..10);#_Alois P.Heinz,2013年7月29日
%p#第三个Maple程序:
%p b:=proc(n,m)选项记忆;展开(`if`(n=0,m,
%p b(n-1,ilcm(m,n))*x+b(n-1,m))
%p端:
%pT:=n->(p->seq(系数(p,x,i),i=0..n))(b(n,1)):
%p序列(T(n),n=0..10);#_阿洛伊斯·海因茨,2023年9月5日
%tt[_,0]=1;t[n_,k_]:=和[LCM@@c,{c,子集[Range[n],{k}]}];表[t[n,k],{n,0,8},{k,0,n}]//Flatten(*_Jean-François Alcover_,2013年7月29日*)
%o(圣人)#(以阿洛伊斯·海因茨(Alois P.Heinz)命名)
%o@CachedFunction
%o定义b(n,k):
%o如果k==0:返回[1]
%o w=b(n-1,k),如果k<n else[0]
%o返回w+[b(n-1,k-1)中c的lcm(c,n)]
%o def T(n,k):返回add(b(n,k))
%o压扁([[T(n,k)代表k in(0..n)]代表n in(0..10)])
%o#_Peter Luschny_,2013年7月29日
%Y行总和为A226037。
%Y参考A065567、A096179、A181854。
%K nonn,表
%0、5
%A _彼得·卢什尼,2010年12月6日
|