|
%I #22 Jul 26 2022 10:10:13
%S 1105,2465,10585,62745,278545,449065,825265,1050985,2531845,3224065,
%T 3664585,5632705,6054985,9582145,11119105,12945745,13187665,13992265,
%U 15403285,21584305,22665505,28787185,31692805,36121345,37354465,39353665,40280065,41298985,47006785,60112885,67371265,74165065,84417985
%N Carmichael numbers ending in 5.
%H Amiram Eldar, <a href="/A355305/b355305.txt">Table of n, a(n) for n = 1..10000</a>
%H <a href="/index/Ca#Carmichael">Index entries for sequences related to Carmichael numbers</a>
%t Select[10*Range[0, 10^7] + 5, CompositeQ[#] && Divisible[# - 1, CarmichaelLambda[#]] &] (* _Amiram Eldar_, Jul 07 2022 *)
%o (PARI) Korselt(n) = my(f=factor(n)); for(i=1, #f[, 1], if(f[i, 2]>1||(n-1)%(f[i, 1]-1), return(0))); 1;
%o isok(n) = ((n%10)==5) && !isprime(n) && Korselt(n) && n>1; \\ _Michel Marcus_, Jul 07 2022; after A002997
%o (Python)
%o from itertools import islice
%o from sympy import factorint, nextprime
%o def A355305_gen(): # generator of terms
%o p, q = 3, 5
%o while True:
%o for n in range(p+2+(-p+3)%10, q, 10):
%o f = factorint(n)
%o if max(f.values()) == 1 and not any((n-1) % (p-1) for p in f):
%o yield n
%o p, q = q, nextprime(q)
%o A355305_list = list(islice(A355305_gen(),10)) # _Chai Wah Wu_, Jul 24 2022
%Y Intersection of A002997 and A017329.
%Y Cf. A352970, A354609, A355307, A355309.
%K nonn,base
%O 1,1
%A _Omar E. Pol_, Jul 03 2022
|
