|
|
|
|
9, 91, 1729, 4187, 6533, 8149, 8401, 10001, 11111, 19201, 21931, 50851, 79003, 83119, 94139, 100001, 102173, 118301, 118957, 134863, 139231, 148417, 158497, 166499, 188191, 196651, 201917, 216001, 226273, 231337, 237169, 251251, 287809, 302177
(列表;图表;参考;听;历史;文本;内部格式)
|
|
|
抵消
|
1,1
|
|
链接
|
|
|
例子
|
9是以10为基数的强伪素数。仅检查10^8=1 mod 9是不够的。由于8=1*2^3,我们还需要验证10=1mod 9和10^2=1mod 9。因为它们都等于1,所以我们看到9实际上是以10为基数的强伪素数。
91也是以10为基数的强伪素数。除了检查10^90=1 mod 91,因为90=45*2,我们还检查10^45=-1 mod 91;-1足以满足强伪素数的定义。
99是以10为基数的费马伪素数(参见A005939号)但它不是以10为基数的强伪素数。虽然10^98=1 mod 99,但由于98=49*2,我们必须检查10^49 mod 99。在这里,我们发现不是-1,也不是1,而是10。因此99不在这个序列中。(结束)
|
|
数学
|
strong伪素数Q[b_,n_]:=模块[{rems=表[PowerMod[b,(n-1)/2^expo,n],{expo,0,IntegerExponent[n-1,2]},(rems[[-1]]==1||MemberQ[rems,n-1])&&PowerMod[b,n-1,n]==1];最大值=5000;选择[Complement[Range[2,max],Prime[Range[PrimePi[max]]],strongPseudo-primeQ[10,#]&](*阿隆索·德尔·阿特2018年8月10日*)
|
|
交叉参考
|
|
|
关键字
|
非n
|
|
作者
|
|
|
状态
|
经核准的
|
|
|
|