%I#20 2020年7月14日10:47:56

%S 2，5，6，7，9，10，11，13，19，21，22，25，27，30，33，37，42，43，55，57，58，61，67，70，

%电话：73,75,78,82,85,91,93102105115121123130133145147157163165，

%电话：177187190193205210211213217235253267273277283310

%N和时（SP）数：不能被4整除的数k，使得（a+b）/gcd（k+1，4）对于每个因式分解k=a*b是素数。

%C低于10^n的SP编号：5、31、123、532、2728、15402、98294。。。

%H Amiram Eldar，n的表，a（n）表示n=1..10000</a>

%H清水贤一，<a href=“http://ousar.lib.okayama-u.ac.jp/56015“>具有互补除数素数和的正整数的算术，冈山大学数学杂志，第60卷（2018年），第155-164页。

%H Shimizu Kenichi和后藤和夫，<a href=“https://repository.lib.tottori-u.ac.jp/4656（https://repository.lib.tottori-u.ac.jp/4656）“>与二次域相关的整数属性，J.Fac.Educ.Tottori Univ.（Nat.Sci.），第47卷（1998），第5-12页。

%e9是一个SP数，因为9=1*9=3*3和（1+9）/2=5以及（3+3）/2=3都是素数。

%e42是一个SP数，因为42=1*42=2*21=3*14=6*7，而1+42=43、2+21=23、3+14=17和6+7=13都是素数。

%t div[n_]：=开关[Mod[n，4]，0,0,1,2,2,1,3,4]；pQ[n_，div_]：=模[{}，d=除数[n]；m=长度[d]；mm=如果[OddQ[m]，（m+1）/2，m/2]；dd=取[d，mm]；ret=真；做[p=（dd[[k]]+n/dd[[k]）/div；如果[！PrimeQ[p]，ret=False；中断[]]，{k，1，mm}]；ret]；spQ[n_]：=模块[{div1=div[n]，ret}，如果[div1==0，ret=False，ret=pQ[n，div1]]；ret]；s={}；Do[If[spQ[n]，AppendTo[s，n]]，{n，1100}]；秒

%o（PARI）isok（n）={if（n%4，my（d=除数（n））；对于（i=1，ceil（#d/2），如果（！isprime（（d[i]+n/d[i]）/gcd（n+1，4）），返回（0）；）；返回（1）；）

%Y A042968的子序列。

%K非n

%O 1,1号机组

%2018年4月24日，阿尔达米拉