The OEIS mourns the passing of Jim Simons and is grateful to the Simons Foundation for its support of research in many branches of science, including the OEIS.
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!)
A176835 Number of positive integers k for which k^2 - n*k is a square. 4
1, 1, 2, 1, 2, 2, 2, 2, 3, 2, 2, 2, 2, 2, 5, 3, 2, 3, 2, 2, 5, 2, 2, 5, 3, 2, 4, 2, 2, 5, 2, 4, 5, 2, 5, 3, 2, 2, 5, 5, 2, 5, 2, 2, 8, 2, 2, 8, 3, 3, 5, 2, 2, 4, 5, 5, 5, 2, 2, 5, 2, 2, 8, 5, 5, 5, 2, 2, 5, 5, 2, 8, 2, 2, 8, 2, 5, 5, 2, 8, 5, 2, 2, 5, 5, 2, 5, 5, 2, 8, 5, 2, 5, 2, 5, 11, 2, 3, 8, 3 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,3
COMMENTS
The number of solutions to the Diophantine equation x^2 - 4*y^2 = n^2 for x,y >= 0, where k = (1+x)/2 and k^2 - n*k = y^2. - R. J. Mathar, May 02 2010
LINKS
FORMULA
a(n) = A176836(n) - A176837(n).
EXAMPLE
A176835(9) = 3 because the three numbers 9^2 - 9*9 = 0, 12^2 - 9*12 = 36 and 25^2 - 9*25 = 400 are squares.
MAPLE
A176835 := proc(n) N := n^2 ; a := 0 ; for d in numtheory[divisors](N) do d2 := N/d ; x := (d+d2)/2 ; y := (d-d2)/4 ; if type(x, 'integer') and type(y, 'integer') and x>=0 and y >= 0 then a := a+1 ; end if; end do: a ; end proc: seq(A176835(n), n=1..100) ; # R. J. Mathar, May 02 2010
MATHEMATICA
a[n_] := Solve[k >= n && x >= 0 && k^2 - n k == x^2, {k, x}, Integers] // Length;
Table[a[n], {n, 1, 100}] (* Jean-François Alcover, Oct 20 2023 *)
PROG
(PARI) A176835(n) = { my(n2 = n^2, s=0); fordiv(n2, d, my(d2=n2/d, x = (d+d2)/2, y = (d-d2)/4); if((1==denominator(x))&&(1==denominator(y))&&(x>=0)&&(y>=0), s++)); (s); }; \\ Antti Karttunen, Sep 28 2018, after R. J. Mathar's Maple-program
CROSSREFS
Sequence in context: A326034 A332207 A339929 * A353865 A237110 A078704
KEYWORD
nonn
AUTHOR
Jens Voß, Apr 27 2010
EXTENSIONS
More terms from R. J. Mathar, May 02 2010
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 June 15 10:58 EDT 2024. Contains 373407 sequences. (Running on oeis4.)