%I#33 2023年2月2日15:28:42

%S 2,5,7,17,19,29,31,37,41,43,53,61,67,73,79,89,97101103113127137，

%电话：139149151157163173193197199211223233241257269271281，

%电话：2832933073173313373493533636367379389397401409439

%N个素数除以形式为3^k+1的数。

%C此序列可用于计算形式3^k+2的素数候选的P-1值，以帮助素性测试。

%Ca（1）=2除以形式3^k+1的每个数字。这是唯一具有此属性的术语。

%C对于k>2，A000040（k）是一个成员，当且仅当A062117（k）是偶数。-_罗伯特·伊斯雷尔，2018年5月23日

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

%e每个3^k+1的值都是偶数，所以2在序列中。

%对于任何整数k，3^k+1的值都不是3的倍数，因此3不在序列中。

%e3^2+1=10，是5的倍数，所以5在序列中。

%pf:=p->numtheory:-顺序（3，p）：：偶数：

%p f（2）：=真：

%p选择（i素数和f，[2，seq（p，p=5..1000,2）]）；#_罗伯特·伊斯雷尔，2018年5月23日

%t连接[{2}，选择[Range[5，1000，2]，PrimeQ[#]&&EvenQ@MultiplicativeOrder[3，#]&]]（*_Jean-François Alcover_，2023年2月2日*）

%o（PARI）isok（p）=如果（p！=3，m=Mod（3，p））；nb=zn阶（m）；对于（k=1，nb，如果（m^k==Mod（-1，p），返回（1）；）；）；返回（0）；\\_米歇尔·马库斯，2018年5月18日

%o（PARI）列表（lim）=我的（v=列表（[2]），t）；对于因子（n=4，lim\1+1，如果（n[2][，2]==[1]~，my（p=n[1]，m=Mod（3，p）））；对于（k=2，znorder（m，t），m*=3；如果（m==-1，则列表输入（v，p）；断裂）；t=n）；Vec（v）\\_Charles R Greathouse IV_，2018年5月23日

%o（PARI）isok（p）=isprime（p）&&if（p<4，p==2，znorder（Mod（3，p））%2==0）Jeppe Stig Nielsen，2020年6月27日

%o（PARI）isok（p）=！isprime（p）&&return（0）；p<4&&返回（p==2）；s=估价（p-1,2）；模式（3，p）^（（p-1）>>s）=2020年6月27日，Jeppe Stig Nielsen

%Y参考A000040、A034472、A062117、A301917、A320481。

%K非n

%O 1,1号机组

%A _Luke W.Richards_，2018年3月28日