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 A341060 #11 May 16 2023 10:39:35 %S A341060 1,2,10,16,30,74,81,97,489,525,607,1861,4439,26051 %N A341060 Numbers k such that A340179(k) is a multiple of k. %e A341060 a(4) = 16 is a term because A340179(16) = 32 = 2*16. %p A341060 f:= proc(n) local C, s, c; %p A341060 C:=select(t -> igcd(t, n) = 1, [$1..n-1]); %p A341060 s:= convert(C, `+`); %p A341060 add(s mod c, c = C) %p A341060 end proc: %p A341060 select(t -> f(t) mod t = 0, [$1..5000]); %t A341060 A340179[n_] := Total@Mod[#2, #1]& @@ {#, Total@#}& @ Select[Range[n], GCD[#, n] == 1&]; %t A341060 Reap[For[k = 1, k <= 80000, k++, If[Divisible[A340179[k], k], Print[k]; Sow[k]]]][[2, 1]] (* _Jean-François Alcover_, May 16 2023, after _Michael De Vlieger_ in A340179 *) %Y A341060 Cf. A340179, A341059, A341061. %K A341060 nonn,more %O A341060 1,2 %A A341060 _J. M. Bergot_ and _Robert Israel_, Feb 04 2021