%I#22 2020年4月11日14:58:01

%S 1,3,2,5,1,7,2,3,1,11,13,1,2,17,19,1,1,1,23,1,5,1,3,1,29,1,31,2，

%T 1,1,1,1,37,1,1,1,41,1,43,1,1,47,1,7,1,1,1,1,53,1,1,11,1,59,1,61,1,1，

%U 2,1,1,67,1,1,1,71,1,73,1,1,1,11,1,1,1,79,1,3,1,83,1,1,1_1,1,89,1,1,1.1,1,1,97,1

%三角数的N LCM变换。

%H A.Nowicki，<A href=“http://arxiv.org/abs/1310.2416“>强可分性和LCM序列</a>，arXiv:1310.2416[math.NT]，2013。

%H A.Nowicki，<A href=“http://www.jstor.org/stable/10.4169/amer.math.monthly.122.10.958“>强可分性和LCM序列</a>，《美国数学》第122卷（2015年），第958-966页。

%F From _Andrey Zabolotskiy_，2020年4月11日：（开始）

%当n>1时，F a（n）=A051543（n-1）。

%当n>1时，F a（n）=A014963（n+1）。（结束）

%p LCMXfm:=程序（a）局部L，i，n，g，b；

%p L：=nops（a）；

%p g：=数组（1..L，0）；b： =数组（1..L，0）；

%p b[1]：=a[1]；g[1]：=a[1]；

%n从2到L的p do g[n]：=ilcm（g[n-1]，a[n]）；b[n]：=g[n]/g[n-1]；od；

%p L打印（[seq（b[i]，i=1..L）]）；

%p端；

%p t1:=[序列（n*（n+1）/2，n=1..100）]；

%p LCMXfm（t1）；

%t LCMXfm[a_List]：=模块[{L=长度[a]，b，g}，b[1]=g[1]=a[[1]；b[_]=0；g[_]=0；Do[g[n]=LCM[g[n-1]，a[[n]]]；b[n]=g[n]/g[n-1]，{n，2，L}]；数组[b，L]]；

%t LCMXfm[表[n*（n+1）/2，{n，1100}]]（*_Jean-François Alcover_，2017年12月5日，来自Maple*）

%Y参考A051543、A020500、A014963。

%K nonn公司

%O 1,2号机组

%A _N.J.A.Sloane，2016年1月2日