%I#20 2018年3月28日10:28:32

%S 2047121341781217,25,9,91,9133,91,85,151687,15,9,25,9-2121,21，

%电话：169,25217,9121,9,15,49,15,25545,33,9,35,9,39133,39,21451,21,9，

%U 481,9,65,49,25,49、25,51、9,55,9、55,25、57、15481、15,9529,9、33

%N基N的最小强伪素数。

%C a（n）=9当且仅当n==1或8（mod 9）_罗伯特·伊斯雷尔，2018年3月27日

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

%H Eric Weisstein的数学世界，<a href=“http://mathworld.wolfram.com/StrongPseudoprime.html“>强伪素数</a>

%H<a href=“/index/Ps#伪素数”>与伪素数相关的序列的索引项</a>

%p过滤器：=proc（n，b）局部d，s，r；

%p如果isprime（n），则返回假fi；

%p s：=padic:-ordp（n-1,2）；

%p d：=（n-1）/2^s；

%p如果b&^d mod n=1，则返回真fi；

%p代表r从0到s-1 do

%p如果b&^（d*2^r）+1 mod n=0，则返回真fi

%p od；

%p假

%p端程序：

%p f：=proc（b）局部n；

%p代表n，从9乘2 do-if过滤器（n，b）然后返回n fiod

%p端程序：

%p映射（f，[2..100]美元）；#_罗伯特·伊斯雷尔，2018年3月27日

%t sppQ[n_？EvenQ，_]：=假；sppQ[n_？PrimeQ，_]：=假；sppQ[n_，b_]：=模块[{ans=False}，s=整数指数[n-1，2]；d=（n-1）/2^s；如果[PowerMod[b，d，n]==1，ans=True，Do[如果[PowerMod[b，d*2^r，n]==n-1，ans=True]，{r，0，s-1}]]；ans]；leastSPP[b_]：=模块[{k=3}，而[！sppQ[k，b]，k+=2]；k] ；表[leastSPP[n]，{n，2100}]（*在A020229*的Jean-François Alcover之后）

%o（PARI）is_a001262（n，a）={（bittest（n，0）&&！isprime（n）&&n>8）||return；my（s=估值（n-1，2））；if（1==a=Mod（a，n）^（n>>s），return（1））；while（a！=-1&&s--，a=a^2）；a=-1}\\after _M.F.Hasler_ in a001262

%o a（n）=复合材料（c=1，如果（is_a001262（c，n），返回（c）））\\_Felix Fröhlich_，2018年3月28日

%Y参见A001262、A020229、A020230、A020231、A020232、A020233、A020234、A020235。

%K非n

%氧2,1

%A _Amiram Eldar，2018年1月26日