OFFSET
1,2
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..100
FORMULA
a(n) = (1/n) * Sum_{d | n} totient(n/d) * A002426(d). - Andrew Howroyd, Mar 02 2017
EXAMPLE
All solutions for n=5:
.-1...-1...-1...-1...-1...-1...-1...-1....0...-1...-1
.-1....1....0...-1....0....0....0....0....0....1...-1
..1...-1....0....1....1...-1....1....0....0....0....0
..1....1....1....0...-1....1....0....0....0....0....1
..0....0....0....1....1....1....0....1....0....0....1
MATHEMATICA
comps[r_, m_, k_] := Sum[(-1)^i*Binomial[r - 1 - i*m, k - 1]*Binomial[k, i], {i, 0, Floor[(r - k)/m]}]; a[n_Integer, k_] := DivisorSum[n, EulerPhi[n/#] comps[#*(k + 1), 2 k + 1, #] &]/n; a[n_] = a[n, 1]; Array[a, 31] (* Jean-François Alcover, Nov 01 2017, after Andrew Howroyd *)
CROSSREFS
KEYWORD
nonn
AUTHOR
R. H. Hardin, Feb 29 2012
STATUS
approved