%I#21 2021年8月20日04:23:08

%S 1,2,3,4,6,10,12,14,15,20,21,24,30,35,36,39,40,42,44,52,55,56,60,65，

%电话66,70,72,76,77,78,84,85,90,91,95,99102105110114115117119120，

%电话：126130132133136138140143152153154155156161165170

%N除以A000793（k）的数字k。

%C在1之后，是A225649的子集。

%C对于所有n，A225650（a（n））=a（n。

%H Alois P.Heinz，n的表格，n=1..10000的a（n）</a>

%p b:=proc（n，i）选项记住；局部p；

%p p:=`if`（i<1，1，ithprime（i））；

%p`if`（n=0或i<1，1，最大值（b（n，i-1），

%p序列（p^j*b（n-p^j，i-1），j=1.ilog[p]（n））

%p端：

%p g：=n->b（n，`if`（n<8，3，数字理论[pi]（细胞（1.328*isqrt（n*ilog（n））））：

%p a:=proc（n）选项记忆；局部k；

%p表示从1+`if`（n=1,0，a（n-1））开始的k

%p而不是irem（g（k），k）=0做od；k个

%p端：

%p序列（a（n），n=1..70）；#_Alois P.Heinz，2013年5月22日

%t Reap[For[n=1，n<=40，n++，If[Divisible[Max[LCM@@@IntegerPartitions[n]]，n]，Sow[n]]][[2，1]]

%t（*或，对于大量术语：*）

%tb[n_，i_]：=b[n，i]=模[{p}，p=如果[i<1，1，素数[i]]；如果[n==0||i<1，1，Max[b[n，i-1]，表[p^j*b[n-p^j，i-1]{j，1，Log[p，n]//Floor}]]]；g[n_]：=b[n，如果[n<8，3，PrimePi[Ceiling[1.328*Sqrt[n*Log[n]//Floor]]]]；收获[对于[k=1，k<=1000，k++，如果[可分割[g[k]，k]，母猪[k]][[2，1]]（*Jean-François Alcover_，2016年2月28日，在_Alois P.Heinz_*之后）

%o（与Antti Karttunen_的IntSeq-library合作的方案）：

%o（定义A225651（匹配POS 1 1（λ（i）（零？（模（A000793 i）i））））

%Y参见A225648、A225649、A225650、A225653、A225655-A225657。

%K nonn公司

%O 1,2号机组

%2013年5月16日，安蒂·卡图内