|
%I #15 Sep 25 2020 05:39:26
%S 49,91,121,133,169,175,217,231,247,259,301,325,341,343,361,385,403,
%T 427,435,451,469,475,481,511,529,553,559,561,589,595,637,645,651,671,
%U 679,703,715,721,763,775,781,793,805,817,841,847,861,871,889,891,925,931,949,961,973,1001,1015,1027,1035,1045
%N Odd composite numbers q such that there exists a, 2<=a<=q-2, such that a^d == 1 mod q where d = A000265(q-1). Thus q is a strong pseudoprime in base a.
%C Odd composite numbers q such that gcd(A000010(q), A000265(q-1)) > 1. - _Robert Israel_, Dec 20 2017
%H Robert Israel, <a href="/A178705/b178705.txt">Table of n, a(n) for n = 1..10000</a>
%F a^d == 1 mod q
%e 18^3 == 1 mod 49
%p filter:= proc(n)
%p if isprime(n) then return false fi;
%p igcd((n-1)/2^padic:-ordp(n-1,2), numtheory:-phi(n)) > 1
%p end proc:
%p select(filter, [seq(i,i=9..2000,2)]); # _Robert Israel_, Dec 20 2017
%t filterQ[n_] := If[PrimeQ[n], False, GCD[(n-1)/2^IntegerExponent[n-1, 2], EulerPhi[n]] > 1];
%t Select[Range[9, 2000, 2], filterQ] (* _Jean-François Alcover_, Sep 25 2020, after _Robert Israel_ *)
%Y Cf. A000010, A000265.
%K nonn
%O 1,1
%A _Karsten Meyer_, Dec 26 2010
%E Corrected by _Robert Israel_, Dec 20 2017
|
