OFFSET
1,1
COMMENTS
Primitive means that the 4 summands are coprime.
Not every term is the sum of two coprime cubes.
a(1) = A047696(2).
LINKS
Robert Israel, Table of n, a(n) for n = 1..1000 (first 352 terms from Rosalie Fay)
EXAMPLE
189 = 4^3 + 5^3 = 6^3 + (-3)^3 and 4, 5, 6, -3 are coprime, so 189 is in the sequence.
35208 = 34^3 + (-16)^3 = 33^3 + (-9)^3 and 34, -16, 33, -9 are coprime, so 35208 is in the sequence.
MAPLE
g:= proc(s, n) local x;
x:= s/2 + sqrt(12*n/s-3*s^2)/6;
if not x::integer then return NULL fi;
[x, s - x];
end proc:
filter:= proc(n)
local pairs, i, j;
pairs:= map(g, numtheory:-divisors(n), n);
for i from 2 to nops(pairs) do
for j from 1 to i-1 do
if igcd(op(pairs[i]), op(pairs[j]))=1 then return true fi
od od;
false
end proc:
select(filter, [seq(seq(9*i+j, j=[1, 2, 7, 8, 9]), i=0..4000)]); # Robert Israel, Oct 22 2017
MATHEMATICA
g[s_, n_] := Module[{x}, x = s/2 + Sqrt[12*n/s - 3*s^2]/6; If[!IntegerQ[x], Return[Nothing]]; {x, s - x}];
filter[n_] := Module[{pairs, i, j}, pairs = g[#, n]& /@ Divisors[n]; For[i = 2, i <= Length[pairs], i++, For[j = 1, j <= i - 1, j++, If[GCD @@ Join[pairs[[i]], pairs[[j]]] == 1, Return[True]]]]; False];
Select[Flatten[Table[Table[9*i + j, {j, {1, 2, 7, 8, 9}}], {i, 0, 4000}]], filter] (* Jean-François Alcover, May 28 2023, after Robert Israel *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Rosalie Fay, Oct 16 2017
STATUS
approved