%I#38 2024年8月2日18:15:13

%S 3111129,3777331311112931311193190141848731079317731129，

%电话：311064731038194250793103975113310793109313731036131259331，

%电话：36389191176757654383310310333752332494131393127139337631034567124791323697195623333711317773493103934412909

%N A316941中应用于奇数复合数（A071904）的迭代过程（如果没有达到素数，则a（N）=-1）。

%Ca（n）=A316941（A071904）。

%e第一个条目是从9=3*3开始的。33=3*11，311是素数。

%e较长的10步进展是从45开始的A（9）。明确地，45=3*15连接到315=3*105连接到3105=3*1035连接到31035=3*10345连接到310345=5*62069连接到562069=41*13709连接到4113709=19*216511连接到19216511=17*1130383=3*57043461连接到357043461=3*119014487连接到3119014487，即素数。a（9）则为3119014487。

%t映射[NestWhile[Function[n，FromDigits@Flatten@IntegerDigits@{#，n/#}&[FactorInteger[n][[1，1]]]，#！PrimeQ@#&]&，选择[Range[9，157，2]，CompositeQ]]（*_Michael De Vlieger_，2016年9月13日*）

%o（PARI）genit（iend）={i5=9；while（i5<=iend，n=i5；where（i素数（n），n+=2）；i5=n；无止境=0；whiles（无止境<99999，dun=0；z=除法（n）；

%o a=z[2]；b=不适用；k=长度（数字（b））；q=a*10^k+b；如果（i素数（q），dun=1；断裂）；无止境+=1；n＝q）；如果（dun>0，打印1（q，“，”））；i5+=2）；}

%o（Python）

%o来自sympy import primepi，primefactors，factorint

%o定义A276662（n）：

%o如果n==1：返回311

%o m，k=n，素数（n）+n+（n>>1）

%o而m！=克：

%o m，k=k，素数（k）+n+（k>>1）

%o当总和（（f:=因子（m））.values（））>1时：

%o m=int（str（p:=最小值（f））+str（m//p））

%o返回m#_Chai Wah Wu_，2024年8月2日

%Y参见A037274、A071904、A316941。

%K nonn，基础

%O 1,1号机组

%A Bill McEachen，2016年9月11日

%E编辑：N.J.A.Sloane，2016年10月2日