OFFSET
1,1
COMMENTS
If k is the average of two positive cubes, then k^4 is also the average of two positive cubes. So this sequence focuses on the solutions that are not trivial.
LINKS
Robert Israel, Table of n, a(n) for n = 1..502
EXAMPLE
329 is a term because 329 is not the average of two positive cubes while 329^4 = (1833^3 + 2585^3)/2.
MAPLE
Q:= proc(x) local t;
for t in select(t -> t^3<=x and 4*t^3 > x and x/t - t^2 mod 3 = 0,
numtheory:-divisors(x)) do
if issqr((x/t - t^2)/3) then return true fi
od:
false
end proc:
select(x -> not(Q(x)) and Q(x^4), [$1..10000]); # Robert Israel, May 26 2016
MATHEMATICA
Q[x_] := Module[{s, t}, s = Select[Divisors[x], #^3 <= x && 4*#^3 > x && Mod[x/# - #^2, 3] == 0 &]; For[t = 1, t <= Length[s], t++, If[IntegerQ@Sqrt[(x/s[[t]] - s[[t]]^2)/3], Return[True]]]; False];
Reap[For[x = 1, x <= 10000, x++, If[!Q[x] && Q[x^4], Print[x]; Sow[x]]]][[2, 1]] (* Jean-François Alcover, May 18 2023, after Robert Israel *)
PROG
(PARI) isA003325(n) = for(k=1, sqrtnint(n\2, 3), ispower(n-k^3, 3) && return(1));
lista(nn) = for(n=1, nn, if(isA003325(2*n^4) && !isA003325(2*n), print1(n, ", ")));
(PARI) T=thueinit('z^3+1);
isA003325(n)=#select(v->min(v[1], v[2])>0, thue(T, n))>0
is(n)=isA003325(2*n^4) && !isA003325(2*n) \\ Charles R Greathouse IV, May 27 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
Altug Alkan, May 26 2016
STATUS
approved