%I#20 2019年3月5日08:24:43

%S 2,27，911031631992032247305371377397421551559757779，

%电话：883991102110791087112311891199123713511521154315671609，

%电话：165117531769179918071873188319192009209207122612539

%N不除A000213的任何项的非平凡整数（即没有具有此属性的非平凡因子）。

%H Robert Israel，n的表，n=1..1275的a（n）</a>

%p nd:=proc（p）局部a，b，c，r，r；

%pa:=1；b： =1；c： =1；R[1,1,1]：=真；

%p do公司

%p r：=a+b+c mod p；

%p如果r＝0，则返回false fi；

%p a：=b；b： =c；c： =r；

%p如果分配（R[a，b，c]）或nops（{a，b、c}）=1

%p然后返回true

%p其他R[a，b，c]：=真

%p fi；

%日期

%p端程序：

%p N:=10^4:#获取所有术语<=N

%p V：=矢量（N）：结果：=空：

%p代表n从1到n do

%p如果V[n]=0，则

%p如果nd（n），则Res:=Res，n；V[[seq（k*n，k=2..楼层（n/n））]:=1；fi（菲涅耳）

%p fi；

%日期：

%p资源；#_罗伯特·伊斯雷尔，2017年2月26日

%t非除数[n_]：=模[{a=1，b=1，c=1，t}，对于[i=1，i<=n^2，i++，t=Mod[a+b+c，n]；如果[t！=0，a=b；b=c；c=t，则返回[False]]；如果[c==1&b==1&&a==1，则返回[True]]]；

%tokQ[n_]：=Do[If[nondivisor[d]，Return[n==d]]，{d，Divisors[n]}]；

%t选择[Range[3000]，okQ]（*Jean-François Alcover_，2019年3月5日，来自PARI*）

%o（PARI）非除数（n）=my（a=1，b=1，c=1，t）；对于（i=1，n^2，t=（a+b+c）%n；如果（t，a=b；b=c；c=t，返回（0））；如果（c==1&b==1&&a==1，返回（1））

%o是（n）=fordiv（n，d，if（nondivisor（d），return（n==d）））；2012年8月29日，查尔斯·格里特豪斯四世

%K nonn公司

%O 1，1

%A·热心的W·威尔逊_

%亨利·阿尤拉（Henry.Ayoola（AT）googlemail.com）修正的E定义，2009年2月3日

%E a（1）由_Charles R Greathouse IV_于2012年8月29日添加