login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A308395 Numbers y such that x*(x+1) + y*(y+1) = z*(z+1) is solvable in positive integers x, z with x <= y. 2
2, 5, 6, 9, 10, 13, 14, 17, 18, 20, 21, 22, 24, 25, 26, 27, 29, 30, 33, 34, 35, 37, 38, 39, 41, 42, 44, 45, 46, 48, 49, 50, 51, 53, 54, 55, 56, 57, 58, 61, 62, 65, 66, 68, 69, 70, 73, 74, 75, 76, 77, 78, 80, 81, 82, 84, 85, 86, 89, 90, 91, 92, 93, 94, 95, 97 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
LINKS
EXAMPLE
14 is a term because 14*15 + 14*15 = 20*21.
MATHEMATICA
max = 220; lst = {}; For[x = 1, x < max, x++,
For[y = x, y < max, y++,
For[z = y, z < max, z++,
If[x (x + 1) + y (y + 1) == z (z + 1),
lst = AppendTo[lst, y]]]]]; Select[Union[lst], # < max/2 &]
PROG
(Python)
from sympy import integer_nthroot
A308395_list, y, w = [], 1, 0
while len(A308395_list) < 10000:
w += y
z = 0
for x in range(1, y+1):
z += x
if integer_nthroot(8*(w+z)+1, 2)[1]:
A308395_list.append(y)
break
y += 1 # Chai Wah Wu, Aug 02 2019
CROSSREFS
Cf. A012132.
Sequence in context: A340289 A133759 A188258 * A227149 A264120 A042963
KEYWORD
nonn
AUTHOR
Ralf Steiner, Jul 31 2019
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 04:13 EDT 2024. Contains 371235 sequences. (Running on oeis4.)