cp's OEIS Frontend

A185635 Fixed points of A075075.

%I A185635 #26 Oct 26 2024 21:17:38
%S A185635 1,2,8,36,39,49,370,626,632,1030,13155,32317,61358,86704,2535431,
%T A185635 13360009
%N A185635 Fixed points of A075075.
%C A185635 A075075(a(n)) = A160516(a(n)) = a(n);
%C A185635 A075075(a(n)-1) * A075075(a(n)+1) mod a(n) = 0.
%o A185635 (Haskell)
%o A185635 a185635 n = a185635_list !! (n-1)
%o A185635 a185635_list = filter (> 0) $
%o A185635    zipWith (\x y -> if x == y then y else 0) [1..] a075075_list
%o A185635 -- _Reinhard Zumkeller_, Dec 19 2012
%o A185635 (Python)
%o A185635 from math import gcd
%o A185635 A185635_list, l1, l2, m, b = [1, 2], 2, 1, 2, {1, 2}
%o A185635 for n in range(3, 10**4):
%o A185635     i = m
%o A185635     while True:
%o A185635         if not i in b:
%o A185635             if n == i:
%o A185635                 A185635_list.append(i)
%o A185635             l1, l2, m = i, l1, i//gcd(l1, i)
%o A185635             b.add(i)
%o A185635             break
%o A185635         i += m # _Chai Wah Wu_, Dec 09 2014
%K A185635 nonn,more
%O A185635 1,2
%A A185635 _Reinhard Zumkeller_, Dec 19 2012
%E A185635 a(15) from _Chai Wah Wu_, Dec 10 2014
%E A185635 a(16) from _Chai Wah Wu_, Sep 22 2019