%I#22 2019年11月17日09:33:19

%S 9,9117294187653381498401101111119201219315085179003，

%电话：831199413910000110217311801118957134863139231148417158497，

%电话：166499188191196651201917216001226273231337237169251287809302177

%N以10为基数的强伪素数。

%H Amiram Eldar，n表，n=1..100000的a（n）（R.J.Mathar的术语1..203）

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

%e From _Alonso del Arte，2018年8月10日：（开始）

%e9是碱基10的强伪素数。仅检查10^8=1 mod 9是不够的。由于8=1*2^3，我们还需要验证10=1mod 9和10^2=1mod 9。因为它们都等于1，所以我们看到9实际上是以10为基数的强伪素数。

%e91也是以10为基数的强伪素数。除了检查10^90=1 mod 91，因为90=45*2，我们还检查10^45=-1 mod 91；-1足以满足强伪素数的定义。

%e99对于基数10是费马伪素数（参见A005939），但对于基数10不是强伪素数。虽然10^98=1 mod 99，但由于98=49*2，我们必须检查10^49 mod 99。在这里，我们发现不是-1，也不是1，而是10。因此99不在这个序列中。（结束）

%t 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，#]&]（*_Alonso del Arte_2018年8月10日*）

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

%K nonn公司

%O 1,1号机组

%大卫·W·威尔逊_