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 A357241 #10 Nov 01 2022 13:41:03
%S A357241 1,2,3,1,4,5,6,1,1,7,8,2,9,10,11,1,12,2,13,3,14,15,16,2,1,17,1,4,18,
%T A357241 19,20,1,21,22,23,1,24,25,26,3,27,28,29,5,3,30,31,2,1,2,32,6,33,2,34,
%U A357241 4,35,36,37,7,38,39,4,1,40,41,42,8,43,44,45,1,46,47,3,9,48,49,50,3
%N A357241 a(n) is the number of j in the range 1 <= j <= n such that j / rad(j) = n / rad(n).
%H A357241 Michel Marcus, <a href="/A357241/b357241.txt">Table of n, a(n) for n = 1..10000</a>
%F A357241 a(n) = |{j <= n : A003557(j) = A003557(n)}|.
%t A357241 Table[Length[Select[Range[n], #/Last[Select[Divisors[#], SquareFreeQ]] == n/Last[Select[Divisors[n], SquareFreeQ]] &]], {n, 1, 80}]
%o A357241 (PARI) f(n) = n/factorback(factor(n)[, 1]); \\ A003557
%o A357241 a(n) = my(x=f(n)); sum(j=1, n, f(j) == x); \\ _Michel Marcus_, Sep 20 2022
%Y A357241 Cf. A001694 (positions of 1's), A003557, A008479.
%K A357241 nonn
%O A357241 1,2
%A A357241 _Ilya Gutkovskiy_, Sep 19 2022