OFFSET
1,3
COMMENTS
Former name was "Number of partitions of 1 into rational parts a/b satisfying 0<a<b<=n", but that was wrong, that is A269926, a different sequence. - N. J. A. Sloane, Dec 20 2024
FORMULA
a(n) = Sum_{k=2..n} A000837(k).
MATHEMATICA
p[n_] := p[n] = IntegerPartitions[n];
l[n_] := Length[p[n]];
p[n_, j_] := p[n, j] = Part[p[n], j]
g[n_, j_] := g[n, j] = Apply[GCD, p[n, j]]
h[n_] := h[n] = Table[g[n, j], {j, 1, l[n]}]
c[n_] := c[n] = Count[h[n], 1]
Table[c[n], {n, 0, 45}] (* A000837 *)
s[n_] := Sum[c[k], {k, 1, n}]
Table[s[n] - 1, {n, 1, 45}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Clark Kimberling, Mar 09 2012
EXTENSIONS
New name and edits made by Clark Kimberling, Dec 18 2024
STATUS
approved