OFFSET
0,4
COMMENTS
This sequence is in reply to an extension request made in A100450.
Note that gcd(0,m) = m for any m.
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..1000 (terms n = 0..200 from Robert Price)
FORMULA
a(n) = (A213208(n) + 1)/2.
MAPLE
h:= proc() true end:
b:= proc(n) local c, i, j, p;
c:=0;
for i to iquo(n, 3) do
for j from i to iquo(n-i, 2) do
if igcd(i, j, n-i-j)=1 then p:= i*j*(n-i-j);
if h(p) then h(p):= false; c:=c+1 fi
fi
od
od; c
end:
a:= proc(n) a(n):= `if`(n=0, 1, a(n-1) +b(n)) end:
seq(a(n), n=0..60); # Alois P. Heinz, Mar 02 2013
MATHEMATICA
f[n_] := Length[ Union[ Flatten[ Table[ If[ i+j+k <= n&& GCD[i, j, k] <= 1, i*j*k, 0], {i, 0, n}, {j, 0, n}, {k, 0, n}], 2]]]; Table[ f[n], {n, 0, 200}]
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Robert Price, Mar 02 2013
STATUS
approved