cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A371075 Fixed points of A374941.

This page as a plain text file.
%I A371075 #28 Aug 25 2024 22:23:56
%S A371075 12,30,80,324,448,888,11264,53248,1114112,3684352,4980736,18055168,
%T A371075 96468992
%N A371075 Fixed points of A374941.
%F A371075 Equals the ordered set {k: A374941(k) = k}.
%e A371075 12 is a term because A374941(12) = 12.
%o A371075 (PARI) f(n) = my(d=divisors(n)); sum(i=2, #d-1, if (isprime(d[i]), d[i], f(d[i])));
%o A371075 isok(k) = f(k) == k; \\ _Michel Marcus_, Mar 10 2024
%o A371075 (Python)
%o A371075 from sympy import divisors, isprime
%o A371075 from functools import cache
%o A371075 @cache
%o A371075 def f(n): return sum(di if isprime(di) else f(di) for di in divisors(n)[1:-1])
%o A371075 def ok(n): return n == f(n)
%o A371075 print([k for k in range(1, 1000) if ok(k)]) # _Michael S. Branicky_, Mar 31 2024 after _Michel Marcus_
%Y A371075 Cf. A000396, A005101, A299795, A374941.
%K A371075 nonn,more
%O A371075 1,1
%A A371075 _Tanmaya Mohanty_, Mar 10 2024
%E A371075 a(13) from _Michael S. Branicky_, Mar 31 2024