%I#52 2023年7月9日02:43:42

%第1,2,5,18,68283116148431998583074349670146586107257页，

%电话255478351061156554403961131833079809764292461231705433101，

%电话：13171160795654628372949322574622982349339325821411385937083186901596000664153613715159245162736805917843710

%N k≤P（N）的数量，即rad（k）|P（N。

%H Bert Dobbelaere，Python程序</a>

%F a（n）=A010846（A002110（n））。

%F a（n）>=2^n。

%e a（0）=1，因为P（0）=1和1|1。

%e a（1）=2，因为P（1）=2，以及1|2和2|2。

%e a（2）=5，因为P（2）=6，rad（m）|6，因为m={1，2，3，4，6}。

%e a（3）=18，因为P（3）=30，rad（m）|30，对于m={1，2，3，4，5，6，8，9，10，12，15，16，18，20，24，25，27，30}等。

%e关于a（3），我们看到在素数幂范围2、3和5的张量积中有18项不超过30：

%电子邮箱5^0X|2^0 2^1 2^2 2^3 2^4 5^1X|2 ^0 2 ^1 2 ^2 5^2X|2

%e-----------------------------------------------

%电子邮箱3^0|1 2 4 8 16 3^0|5 10 20 3^0| 25

%电子邮箱3^1|3 6 12 24 3^1|15 30

%e 3^2|9 18

%e 3^3|27

%e因此，a（3）=18。事实证明，这种方法对于较大的n很方便。

%t f[1]=1；f[n_]：=函数[w，

%t ToExpression@StringJoin[“块[{n=”，ToString@n，

%t“，k=0}，展平@表[k++，”，

%t大多数@Flatten@Map[{#，“，”}&，#]，“]；k]”]&@

%t映射索引[

%t函数[p，StringJoin[“{”，ToString@Last@p，“，0，Log[”，

%t ToString@First@p，“，n/（”，

%t ToString@InputForm[Times@@Map[Power@@#&，Take[w，First@#2-1]]，

%t“）]}”]]@w[[First@#2]]&，w]]@

%t映射[{#，ToExpression[“p”<>ToString@PrimePi@#]}&，

%t系数整数[n][[All，1]]]；

%t地图[f，FoldList[Times，1，Prime@Range@9]]

%Y参考A002110、A007947、A010846。

%K nonn，硬

%0、2

%2023年6月16日，A _Michael De Vlieger_

%E修正了2023年6月27日_Bert Dobbelaere中的a（15）并增加了a（16）-a（23）

%E a（24）-a（26）来自Martin Ehrenstein，2023年7月8日