login
Even pseudoprimes to base 9.
14

%I #16 Sep 19 2024 01:50:00

%S 4,8,28,52,286,364,532,616,946,1036,1288,2806,2926,3052,4376,4636,

%T 5356,6364,8744,8866,11476,12124,15964,17446,19096,19684,21196,21736,

%U 24046,24388,26596,31876

%N Even pseudoprimes to base 9.

%H Amiram Eldar, <a href="/A090083/b090083.txt">Table of n, a(n) for n = 1..10000</a> (terms 1..1000 from Charles R Greathouse IV)

%H <a href="/index/Ps#pseudoprimes">Index entries for sequences related to pseudoprimes</a>.

%t Do[s=Mod[ -1+9^(n-1), n]; If[Equal[s, 0]&&!PrimeQ[n]&&EvenQ[n], Print[n]], {n, 1, 1000000}]

%o (PARI) is(n)=Mod(9, n)^(n-1)==1&&!isprime(n)&&n%2==0 \\ _Charles R Greathouse IV_, Apr 12 2012

%o (PARI) p=2; forprime(q=3, 1e8, forstep(n=p+1, q-1, 2, if(Mod(9, n)^(n-1)==1, print1(n", "))); p=q) \\ _Charles R Greathouse IV_, Apr 12 2012

%Y Cf. A020138.

%Y Cf. A006935, A130433, A090082, A090084, A090085, A130434, A130435, A130436, A130438, A130439, A130440, A130441, A130442, A130443.

%K nonn,changed

%O 1,1

%A _Labos Elemer_, Nov 25 2003