A380034 Records in A380032.
1, 2, 3, 4, 5, 6, 7, 9, 12, 13, 14, 15, 17, 21, 28, 31, 35, 41, 42, 45, 51, 54, 60, 61, 67, 75, 89, 90, 111, 118, 133, 147, 155, 163, 176, 185, 186, 195, 205, 231, 246, 277, 307, 323, 343, 368, 369, 377, 383, 411, 429, 455, 471, 479, 490, 502, 545, 565, 567, 627
Offset: 1
Keywords
Links
- Michael De Vlieger, Table of n, a(n) for n = 1..66
Programs
-
Mathematica
(* Load function f at A025487 *) r = 0; s = Select[Union@ Flatten@ f[8][[3 ;; -1]], Not @* SquareFreeQ]; nn = Length[s]; Print[nn]; Reap[Monitor[ Do[k = s[[i]]; If[# > r, r = #; Sow[r]] &@ Count[Transpose@ {#, k/#} &@ #[[2 ;; Ceiling[Length[#]/2]]] &@ Divisors[k], _?((m = GCD @@ {##}; And[! MemberQ[{1, #2}, m], m == #1, ! Divisible[#1, rad[#2]]]) & @@ # &)], {i, nn}], i] ][[-1, 1]]
Comments