%I#8 2022年3月14日17:36:15
%S 1,0,1,0,1,2,0,1,9,6,0,0,17,36,12,0,-2,50325300,60,0,02824755250,
%电话:2700360,060,-882141617717580850264602520,0,0,-60854882392667,
%电话:499800311640705605040,0,-5046480,-5740473558835637213969000164052027216015120
%N按行读取的三角形。设R(n,k)=Y(n,k,B),其中Y是部分贝尔多项式,B是列表[Bernoulli(j,1),j=0..n]。T(n,k)是由第n行(A048803)中各项分母的lcm归一化的R(n,k)。
%F T(n,n)=产品{k=1..n}拉德(k)=产品_{k=1..n}A007947(k)=A048803(n)。
%e三角形开始:
%e[0]1;
%e[1]0,1;
%e[2]0,1,2;
%e[3]0,1,9,6;
%e[4]0,0,17,36,12;
%e[5]0,-2,50,325,300,60;
%e【6】0、0、28、2475、5250、2700、360;
%e【7】0、60、-882、14161、77175、80850、26460、2520;
%电子[8]0,0,-608,5488,239267,499800,311640,70560,5040;
%e、。
%e例如,第7行是2520*[R(7,k),k=0..7]=2520*[0,1/42,-7/20,2023/360,245/8,385/12,21/2,1],因为lcm(1,42,20,360,8,12,2,1)=A048803(7)=2520。相反,由于R(n,n)=1和T(n,n)=Product_{k=1..n}rad(k),R(n、k)可以通过将行n的项除以T(n、n)得到。
%tB[n_,k_]:=BellY[n,k,表[BernoulliB[j,1],{j,0,n}]];
%t P[n_]:=选择[Divisors[n],PrimeQ];
%tT[n_,k_]:=B[n,k]乘积[乘积[p,{p,p[j]}],{j,1,n}];
%t表[t[n,k],{n,0,9},{k,0,n}]//展平
%Y参见A027641/A027642、A048803、A007947、A264428(贝尔变换)。
%K符号,tabl
%0、6
%阿佩特·卢什尼,2022年3月14日
|