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