OFFSET
1,5
LINKS
FORMULA
a(n) = Sum_{k=1..floor(n/3)} Sum_{i=k..floor((n-k)/2)} floor(1/gcd(k,i)) * floor(1/gcd(k,n-i-k)).
MAPLE
f:= proc(n) local k, t;
t:= 0:
for k from 1 to floor(n/3) do
t:= t + nops(select(s -> igcd(s, k) = 1 and igcd(n-s, k) = 1, [$k..(n-k)/2]))
od;
t
end proc:
map(f, [$1..100]); # Robert Israel, Nov 25 2024
MATHEMATICA
Table[Sum[Sum[Floor[1/GCD[k, i]]*Floor[1/GCD[k, n - i - k]], {i, k, Floor[(n - k)/2]}], {k, Floor[n/3]}], {n, 80}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Wesley Ivan Hurt, Jan 02 2021
STATUS
approved