The OEIS is supported by the many generous donors to the OEIS Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A357241 a(n) is the number of j in the range 1 <= j <= n such that j / rad(j) = n / rad(n). 1
 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, 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, 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 (list; graph; refs; listen; history; text; internal format)
 OFFSET 1,2 LINKS Michel Marcus, Table of n, a(n) for n = 1..10000 FORMULA a(n) = |{j <= n : A003557(j) = A003557(n)}|. MATHEMATICA Table[Length[Select[Range[n], #/Last[Select[Divisors[#], SquareFreeQ]] == n/Last[Select[Divisors[n], SquareFreeQ]] &]], {n, 1, 80}] PROG (PARI) f(n) = n/factorback(factor(n)[, 1]); \\ A003557 a(n) = my(x=f(n)); sum(j=1, n, f(j) == x); \\ Michel Marcus, Sep 20 2022 CROSSREFS Cf. A001694 (positions of 1's), A003557, A008479. Sequence in context: A253887 A026370 A078446 * A055447 A055448 A055449 Adjacent sequences: A357238 A357239 A357240 * A357242 A357243 A357244 KEYWORD nonn AUTHOR Ilya Gutkovskiy, Sep 19 2022 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified July 19 15:03 EDT 2024. Contains 374410 sequences. (Running on oeis4.)