%我#18 2012年3月31日10:28:26

%S 3,5,11,13,19,37,43,53,61,67,83101107131139149163173，

%电话：18121127728330731733449373397421461467491499509，

%电话：523541547557563571587613619653659661691709733739

%N除A003499（（p-1）/2）+6而不除A00349（N）+6的素数p，其中N<（p-1”/2。

%C已检查了所有小于3000000的奇数，似乎形式8N+/-3的大多数素数（例如3,5,11,13,19,37…）都符合测试，但有一些例外（例如29），形式8N+/-1的素数或复合素数都不符合测试。

%t f=编译[{{Co，_Integer}，{S0，_Integer}，}，[S1，_Iinteger}，[Caa，_Integer}}，模块[{xCo=Co，xS0=S0，xS1=S1，Temp}，而[Temp=Mod[6 xS1-xS0-6，Caa]；xCo>0&&温度>0，xS0=xS1；xS1=温度；xCo-]；xCo]]；Caa=5；收获[While[Caa<1000，Co=（Caa-3）/2；S0=2；S1=3；If[Co，S0，S1，Caa]==1，Sow[Caa]]；Caa+=2]]（为了简化代码，跳过了*prime 3。Bill Simpson提供的上述代码比我的原始代码快20倍。请注意，通过不搜索8n+/-1*形式的数字，似乎也可以将速度提高10倍。）

%Y参考A003499、A001541（等于A003499.的1/2）。

%K nonn公司

%O 1,1号机组

%A _Kenneth J Ramsey，2012年2月16日