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.
%I A351087 #9 Feb 05 2022 17:24:18 %S A351087 0,1,7,8,16,25,80,91,161,175,203,343,392,440,539,588,616,968,1100, %T A351087 1225,1331,1333,1337,1375,1400,2556,3025,3773,5112,5875,6655,7280, %U A351087 7668,8281,8575,8600,9604,9800,10224,10820,10868,11011,11492,12595,12980,13013,15925,17303,17576,17875,20449,22308,23677,24067,24167 %N A351087 After initial 0, numbers k such that A327860(k) is a multiple of k. %t A351087 Select[Range[25000], Block[{i, m, n = #, p}, m = i = 1; While[n > 0, p = Prime[i]; m *= p^Mod[n, p]; n = Quotient[n, p]; i++]; # == GCD[#, If[m < 2, 0, m Total[#2/#1 & @@@ FactorInteger[m]]]]] &] (* _Michael De Vlieger_, Feb 04 2022 *) %o A351087 (PARI) %o A351087 A327860(n) = { my(s=0, m=1, p=2, e); while(n, e = (n%p); m *= (p^e); s += (e/p); n = n\p; p = nextprime(1+p)); (s*m); }; %o A351087 isA351087(n) = (!n || !(A327860(n)%n)); %Y A351087 Fixed points of A351083. %Y A351087 Cf. A003415, A276086, A327860. %Y A351087 Cf. A328110 (a subsequence). %K A351087 nonn %O A351087 1,3 %A A351087 _Antti Karttunen_, Feb 03 2022