The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation. Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A154778 Numbers of the form a^2 + 5b^2 with positive integers a,b. 21

%I

%S 6,9,14,21,24,29,30,36,41,45,46,49,54,56,61,69,70,81,84,86,89,94,96,

%T 101,105,109,116,120,126,129,134,141,144,145,149,150,161,164,166,174,

%U 180,181,184,189,196,201,205,206,214,216,224,225,229,230,241,244,245,246

%N Numbers of the form a^2 + 5b^2 with positive integers a,b.

%C Subsequence of A020669 (which allows for a=0 and/or b=0). See there for further references. See A155560 ff for intersection of sequences of type (a^2 + k b^2).

%C Also, subsequence of A000408 (with 5b^2 = b^2 + (2b)^2).

%e a(1) = 6 = 1^2 + 5*1^2 is the least number that can be written as A+5B where A,B are positive squares.

%e a(2) = 9 = 2^2 + 5*1^2 is the second smallest number that can be written in this way.

%t formQ[n_] := Reduce[a > 0 && b > 0 && n == a^2 + 5 b^2, {a, b}, Integers] =!= False; Select[ Range, formQ] (* _Jean-François Alcover_, Sep 20 2011 *)

%t Timing[mx = 300; limx = Sqrt[mx]; limy = Sqrt[mx/5]; Select[Union[Flatten[Table[x^2 + 5 y^2, {x, limx}, {y, limy}]]], # <= mx &]] (* _T. D. Noe_, Sep 20 2011 *)

%o (PARI) isA154778(n,/* use optional 2nd arg to get other analogous sequences */c=5) = { for( b=1,sqrtint((n-1)\c), issquare(n-c*b^2) & return(1))}

%o for( n=1,300, isA154778(n) & print1(n","))

%Y Cf. A033205 (subsequence of primes). [From _R. J. Mathar_, Jan 26 2009]

%K easy,nonn

%O 1,1

%A _M. F. Hasler_, Jan 24 2009

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

Last modified May 31 00:06 EDT 2020. Contains 334747 sequences. (Running on oeis4.)