# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/ Search: id:a058957 Showing 1-1 of 1 %I A058957 #18 Jul 15 2018 08:03:00 %S A058957 9,15,16,21,24,25,27,32,33,35,36,39,40,45,48,49,51,55,56,57,60,63,64, %T A058957 65,69,72,75,77,80,81,84,85,87,88,91,93,95,96,99,100,104,105,108,111, %U A058957 112,115,117,119,120,121,123,125,128,129,132,133,135,136,140,141,143 %N A058957 Numbers having at least two representations as b^2 - c^2 with b > c >= 0. %C A058957 This is the union of the squares > 4 and A306102: numbers that are the difference of two positive squares in at least two ways (where c=0 is excluded). - _M. F. Hasler_, Jul 12 2018 %F A058957 A058957 = { n | A034178(n) >= 2 }. - _M. F. Hasler_, Jul 10 2018 %e A058957 9 is in the sequence since 9 = 3^2 - 0^2 = 5^2 - 4^2; %e A058957 45 is in the sequence since 45 = 7^2 - 2^2 = 9^3 - 6^2 = 23^2 - 22^2. %t A058957 Select[Range@143, Length@ FindInstance[x^2 - y^2 == # && x>y>= 0, {x, y}, Integers, 2] == 2 &] (* _Giovanni Resta_, Jul 10 2018 *) %o A058957 (PARI) select( is(n)=fordiv(n, d, d^2 > n && return; (n-d^2)%(2*d) || c++<2 || return(1)), [1..200]) \\ _M. F. Hasler_, Jul 10 2018 %Y A058957 Cf. A034178, A016825. %Y A058957 Cf. A306102 (subsequence and variant using c > 0). %K A058957 nonn %O A058957 1,1 %A A058957 _Henry Bottomley_, Jan 13 2001 %E A058957 Name edited by _M. F. Hasler_, Jul 10 2018 # Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE