%I #10 Jun 13 2020 01:47:28
%S 370,407,1332,2331,2960,3256,4921,5957,8029,8288,9990,10656,10989,
%T 12691,12913,13357,13949,14023,14911,16021,16354,17353,18648,18907,
%U 19684,19721,20683,22681,23680,24605,24901,26048,27343,30007,30303,32893,34965,35964,36001,36556,37259,39331,39368,39627
%N Numbers which are the sum of two positive cubes and divisible by 37.
%H Robert Israel, <a href="/A102618/b102618.txt">Table of n, a(n) for n = 1..10000</a>
%p N:= 200000: # for terms <= N
%p G:= expand(add(x^(i^3),i=1..floor(N^(1/3)))^2):
%p select(t -> coeff(G,x,t) > 0, [seq(i,i=37..N,37)]); # _Robert Israel_, Jun 12 2020
%t upto[n_] := Block[{t}, Union@ Reap[ Do[If[ Mod[t = x^3 + y^3, 37] == 0, Sow@ t], {x, n^(1/3)}, {y, Min[x, (n - x^3)^(1/3) ]}]][[2, 1]]]; upto[40000] (* _Giovanni Resta_, Jun 12 2020 *)
%Y Cf. A003325. Other sequences of the form "sum of two positive cubes and divisible by ...": A224484, A224485, A101421, A101852, A094447, A099178, A102619, A101806, A224483, A102658.
%K nonn
%O 1,1
%A Jun Mizuki (suzuki32(AT)sanken.osaka-u.ac.jp), Jan 31 2005
%E Corrected by _Robert Israel_, Jun 12 2020
|