|
%I #29 Jun 13 2021 08:46:22
%S 2701,18721,31621,49141,83333,83665,88561,90751,93961,104653,107185,
%T 176149,204001,226801,228241,276013,282133,534061,563473,574561,
%U 622909,653333,665281
%N Composite numbers k such that 2^k-2 and 3^k-3 are both divisible by k and k is not a Carmichael number (A002997).
%H Amiram Eldar, <a href="/A153513/b153513.txt">Table of n, a(n) for n = 1..10000</a> (terms 1..200 from Robert Israel)
%p filter:= proc(n) local p;
%p if isprime(n) or (2 &^n - 2 mod n <> 0) or (3 &^n - 3 mod n <> 0) then return false fi;
%p if n::even then return true fi;
%p if not numtheory:-issqrfree(n) then return true fi;
%p for p in numtheory:-factorset(n) do
%p if n-1 mod (p-1) <> 0 then return true fi
%p od;
%p false
%p end proc:
%p select(filter, [$2..10^6]); # _Robert Israel_, Jan 29 2017
%t Reap[Do[If[CompositeQ[n] && Divisible[2^n-2, n] && Divisible[3^n-3, n] && Mod[n, CarmichaelLambda[n]] != 1, Print[n]; Sow[n]], {n, 2, 10^6}]][[2, 1]] (* _Jean-François Alcover_, Mar 25 2019 *)
%Y Intersection of A153514 and A153508 (excluding the number 1).
%Y Cf. A002997, A122780.
%K nonn
%O 1,1
%A _Artur Jasinski_, Dec 28 2008
|
