%I#54 2023年1月6日15:23:50

%S 1,2,3,5,6,7,9,10,11,13,14,15,17,18,19,21,22,23,25,26,29,30,31,33,34，

%电话：35,37,38,39,41,42,43,45,46,47,49,50,51,53,55,57,58,59,61,62,63,65,66，

%U 67,69,70,71,73,73,74,75,77,78,79,83,85,86,87,89,90,91,91,93,94,95,97,98

%对于任何素数p，N个数都不能被p^p整除。

%C如果a（n）=乘积p_i^e_i，则p_i>e_i对于所有i。

%A100716的C补语；A129251（a（n））=0.-_Reinhard Zumkeller_，2007年4月7日

%C密度为0.72199023441955…=Product_{p>=2}（1-p^-p），其中p覆盖素数。[Charles R Greathouse IV_，2012年1月25日]

%C A027748（a（n），k）<=A124010（a（n]，k），1<=k<=A001221（a（n-））。[_Reinhard Zumkeller_，2012年4月28日]

%H Reinhard Zumkeller，n表，n=1..10000的a（n）</a>

%F a（n）~kn with k=1/Product_{p>=2}（1-p^-p）=Product_{p>=2}（1+1/（p^p-1））=1.3850602852…，其中乘积是所有素数的乘积

%e6=2^1*3^1可以，但12=2^2*3^ 1不行。

%e 625=5^4存在，因为它不能被5^5整除。

%t{1}~加入~选择[Range@120，Times@@Boole@Map[First@#>Last@#&，FactorInteger@#]>0&]（*_Michael De Vlieger_，2016年8月19日*）

%o（哈斯克尔）

%o a048103 n=a048103_列表！！（n-1）

%o a048103_list=过滤器（\x->和$

%o zipWith（>）（a027748_row x）（映射到Integer$a12410_row x））[1..]

%o——_ Inhard Zumkeller_，2012年4月28日

%o（方案，与Antti Karttunen的IntSeq-library合作）

%o（定义A048103（零位置1 1 A129251））

%o_Antti Karttunen，2016年8月18日

%o（PARI）isok（n）=我的（f=系数（n））；对于（i=1，#f~，如果（f[i，1]<=f[i、2]，返回（0））；返回（1）；\\_Michel Marcus，2020年11月13日

%o（Python）

%o从itertools导入计数，islice

%o来自sympy进口保理商

%o定义A048103_gen（startvalue=1）：#术语生成器>=startvalue

%o返回过滤器（lambda n:all（map（lambdad:d[1]<d[0]，factorint（n）.items（））），count（max（startvalue，1））

%o A048103_list=list（岛屿（A048103_ gen（），30））#_Chai Wah Wu_，2023年1月5日

%Y补码：A100716。

%A129251中0的Y位置，A327936和A342007中1的位置。

%Y参考A048102、A048104、A054743、A05474、A359550（特性函数）。

%Y参见A276086（该序列的排列）。

%Y参见A276092（子序列）。

%Y参考A051674（p^p）。

%Y第一次与子序列A276078不同，n=451，其中a（451）=625，而A276079中缺少该值。

%K nonn，简单

%O 1,2号机组

%A _N.J.A.斯隆_

%E来自James A.Sellers_的更多条款，2000年4月22日