cp's OEIS Frontend

A376584 Fixed points in A375564.

%I A376584 #13 Oct 02 2024 06:43:05
%S A376584 1,2,3,4,15,51,63,363,437,729,1385,1388,1413,2807,5418,5422,5428,5435,
%T A376584 5451,5537,10790,10795,10803,10871,21505,21507,21899,22081,42969,
%U A376584 43073,43387,44069,171631,172889,342493,2747813,5464783,10908941,10911583,87259373,87327187,174579089
%N A376584 Fixed points in A375564.
%o A376584 (Python)
%o A376584 from itertools import count, islice
%o A376584 from math import gcd
%o A376584 from sympy import isprime
%o A376584 def A376584_gen(): # generator of terms
%o A376584     aset, a, b = {1,2}, 2, 3
%o A376584     yield from (1,2)
%o A376584     for i in count(3):
%o A376584         c = b
%o A376584         if isprime(a):
%o A376584             while c in aset or gcd(c,a)>1:
%o A376584                 c+=1
%o A376584         else:
%o A376584             while c in aset or gcd(c,a)==1:
%o A376584                 c+=1
%o A376584         aset.add(a:=c)
%o A376584         if i == c:
%o A376584             yield i
%o A376584         while b in aset:
%o A376584             b += 1
%o A376584 A376584_list = list(islice(A376584_gen(),20)) # _Chai Wah Wu_, Sep 30 2024
%Y A376584 Cf. A375564, A376189, A376191, A376192, A376193, A376195.
%K A376584 nonn
%O A376584 1,2
%A A376584 _Scott R. Shannon_, Sep 29 2024
%E A376584 a(40)-a(42) from _Scott R. Shannon_, Oct 02 2024