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.

A341060 Numbers k such that A340179(k) is a multiple of k.

This page as a plain text file.
%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