%I#27 2025年4月22日03:49:06

%S 2,2,3,2,2,3,1,2,5,2,3,5,2,2,5,3,2,7,3,3,2,2,5,1,3,5,1,5,3,1,7,2,5，

%温度2,3,5,2,7,11,2,3,3,7,2,3,2,11,5,5,2,3,3,11,7,5，

%U 2,11,3,2,5,7,2,3,13,2,3,3,3,11,2,7,2,5,13,2,2,7,17,3,5,3,2

%N素数因子<=第N个半素数的其他素数因子。

%C A001358（n）的素因子的较小值。-宋嘉宁，2022年8月5日

%H Zak Seidov，n的表格，n=1..1000的a（n）</a>

%H Eric Weisstein的数学世界，<a href=“https://mathworld.wolfram.com/Semiprime.html“>半素数</a>

%F a（n）=A020639（A001358（n））。

%F a（n）=A001358（n）/A006530（A001359（n））。【由米歇尔·马库斯修订，2020年7月18日】

%F a（n）=A001358（n）/A084127（n）。

%t系数整数[#][[1,1]]&/@Select[Range[500]，PrimeOmega[#]==2&]（*哈维·P·戴尔，2018年6月25日*）

%o（哈斯克尔）

%o a084126=a020639。a001358---Reinhard Zumkeller_，2012年11月25日

%o（Python）

%o来自sympy import primepi，primerange，primefactors

%o定义A084126（n）：

%o定义平分（f，kmin=0，kmax=1）：

%o而f（kmax）>kmax:kmax<<=1

%o kmin=kmax>>1

%o当kmax-kmin>1时：

%o kmid=kmax+kmin>>1

%o如果f（kmid）<=kmid：

%o kmax=kmid

%o其他：

%o kmin=kmid

%o返回kmax

%o定义f（x）：返回int（n+x+（（t:=primepi（s:=isqrt（x）））*（t-1）>>1）-和

%o回归最小值（素数（二等分（f，n，n））#_Chai Wah Wu_，2025年4月3日

%Y参见A001358（半素数），A084127（半素数中素数因子的较大者）。

%Y参见A068318、A087718、A087 794、A089994、A089 995、A096 916、A106550、A106554、A108541、A131284、A138510、A138511。

%K nonn公司

%O 1，1

%A _Reinhard Zumkeller，2003年5月15日