OFFSET
1,5
COMMENTS
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
FORMULA
EXAMPLE
Let s = A005117.
Table of n, s(n), a(n), a(n)+1 for select n:
n s(n) a(n) List of k.
-----------------------------------------------------------------
1 1 0 {}
2 2 1 {1}
5 6 = 2 * 3 4 {1, 2, 3, 8}
7 10 = 2 * 5 2 {1, 4}
19 30 = 2 * 3 * 5 10 {1, 2, 3, 4, 5, 8, 9, 15, 24, 80}
28 42 = 2 * 3 * 7 9 {1, 2, 3, 6, 7, 8, 27, 48, 63}
41 66 = 2 * 3 * 11 7 {1, 2, 3, 8, 11, 32, 242}
49 78 = 2 * 3 * 13 6 {1, 2, 3, 8, 12, 26}
69 110 = 2 * 5 * 11 3 {1, 4, 10}
85 138 = 2 * 3 * 23 5 {1, 2, 3, 8, 23}
521 858 = 2 * 3 * 11 * 13 11 {1, 2, 3, 8, 11, 12, 26, 32, 143, 242, 351}
529 870 = 2 * 3 * 5 * 29 12 {1, 2, 3, 4, 5, 8, 9, 15, 24, 29, 80, 144}
1163 1914 = 2 * 3 * 11 * 29 8 {1, 2, 3, 8, 11, 32, 87, 242}
MATHEMATICA
s = Select[Range[150], SquareFreeQ]; rad[x_] := Times @@ FactorInteger[x][[All, 1]]; f[x_] := Select[Range[x^2, 1, -1], Divisible[x, rad[#]] &]; Table[Count[Partition[f[s[[n]]], 2, 1], _?(#1 - #2 == 1 & @@ # &)], {n, Length[s]}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Michael De Vlieger, May 05 2026
STATUS
approved
