login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A233422
Numbers n such that m - n^3 is a square, where m is the least square above n^3.
1
0, 2, 3, 6, 12, 20, 24, 30, 40, 42, 56, 60, 68, 75, 78, 84, 87, 120, 126, 160, 180, 248, 264, 270, 273, 308, 312, 318, 330, 336, 351, 360, 396, 564, 570, 588, 615, 620, 630, 635, 720, 738, 780, 840, 912, 1008, 1016, 1032, 1284, 1308, 1320, 1334, 1344, 1404, 1540, 1617
OFFSET
1,2
COMMENTS
Numbers n such that A070929(n) is a nonzero square.
The sequence of cubes a(n)^3 begins: 0, 8, 27, 216, 1728, 8000, 13824, 27000, 64000, 74088, 175616, 216000, 314432, ...
The sequence of m's begins: 1, 9, 36, 225, 1764, 8100, 13924, 27225, 64009, 74529, 176400, 216225, 314721, ...
The sequence of square roots of these m's begins: 1, 3, 6, 15, 42, 90, 118, 165, 253, 273, 420, 465, 561, 650, 689, 770, 812, ...
The sequence of squares m-n^3 begins: 1, 1, 9, 9, 36, 100, 100, 225, 9, 441, 784, 225, 289, 625, 169, 196, 841, ...
The sequence of their square roots begins: 1, 1, 3, 3, 6, 10, 10, 15, 3, 21, 28, 15, 17, 25, 13, 14, 29, 35, 43, 24, ... (note the first 12 terms are triangular numbers, A000217).
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
MATHEMATICA
fQ[n_]:=Module[{c=n^3, m}, m=(Floor[Sqrt[c]]+1)^2; IntegerQ[Sqrt[m-c]]]; Select[Range[0, 1650], fQ] (* Harvey P. Dale, Jan 03 2024 *)
PROG
(Python)
from math import isqrt
def isSquare(a):
sr = isqrt(a)
return (a==sr*sr)
for n in range(77777):
n3 = n*n*n
a = isqrt(n3)+1
if isSquare(a*a-n3): print(n, end=', ')
(PARI) is(n)=issquare((sqrtint(n=n^3)+1)^2-n) \\ Charles R Greathouse IV, Dec 09 2013
KEYWORD
nonn
AUTHOR
Alex Ratushnyak, Dec 09 2013
STATUS
approved